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Abstract 


The solar magnetic field, including the interplanetary field is discussed. It is pointed 
out that in the absence of a theory of the Zeeman-effect in a turbulent atmosphere it is at 
present impossible to derive a value of the sun’s general magnetic field from spectroscopic 


measurements. 


From the magneto-hydrodynamic theory of sunspots a value of the sun’s general magne- 
tic field can be derived. Ananthakrisnan’s diagrams of the prominence activity can be inter- 
preted as due to the same disturbance wave as is causing the sunspots. This gives further 
confirmation to the theory, which indicates magnetic field of H, = 12 gauss. 


The magnetic field outside the sun (interplanetary field) is certainly far from a dipole 
field. An approximate model of the field is tentatively suggested. This model is in agree- 
ment with information derived from cosmic ray and magnetic storm data and is also 


reconcilable with the coronal ray structure. 


I. Introduction 


At present the discussion about the general 
magnetic field of the sun is very confused. 
This seems to be due to the almost general 
acceptance of two altogether unfounded as- 
sumptions about the field: 

1. It is usually assumed that outside the solar 
surface the field is necessarily a dipole field. 
This would be true if in this region there were 
no electric current. The conspicuous violent 
disturbances in the corona, the occurrence of 
magnetic storms and aurorae, and the variations 
in the cosmic radiation are clear indications of 
changes in the electromagnetic state of inter- 
planetary space which must be associated with 
electric currents and drastic deviations from 
a dipole field. 

2. Also quite without any justification it is 
assumed that Zeeman-effect measurements give 
an accurate value of the solar magnetic field. 
The experimental technique has advanced in an 
admirable way from the time when Hale’s 
measurements indicated the existence of a 
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general field of + so gauss to the present days 
when Bascock (1955 a, b) claims a value of 
— 0.5 gauss. There is probably no objection to 
the present measurements of the Zeeman-effect 
displacement, but still there does not exist any 
theory of the Zeeman effect in a turbulent atmos- 
phere. The granulation shows that the photo- 
sphere is in a turbulent state and this must affect 
also the magnetic field, so that it becomes 
turbulent. This means that the temperature T, 
the density @ and the magnetic field H are 
coupled and it is not at all likely that the 
average Zeeman effect gives an average of the 
magnetic field. 


The turbulence converts kinetic energy into 
magnetic energy and back again. Under the 
assumption of equipartition the turbulent 
magnetic field was estimated to a few hundred 
gauss (H = 160 gauss—s500 gauss) (ALFVEN, 
1950, p. 146). BATCHELOR (1950) has raised 
some doubts about the equipartition for large 
turbulence elements, but this objection seems 
not to be applicable to the granulation. Com- 
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pare also CHANDRASEKHAR (1950, 1951) and 
KuLsrup (1955). COWLING (1953) has discussed 
the problem from different points of view and 
his results essentially confirms the existence of 
at least an approximate equipartition. 

Thus there are good reasons to suppose that 
in the photosphere there is a turbulent magnetic 
field of the order of, say H = 300 gauss. This 
cannot be detected by the present methods 
because of the lack of resolution. Babcock’s 
measurements refer to surface elements which 
contain many thousand granulae. Suppose that 
in a small surface element do of the photo- 
sphere the magnetic field component in the 
direction of the line of sight is Hz. It seems to 
be generally supposed that the measured field 
is the simple average 


H. = [ H.do | [do (1) 


This would be true in a uniform atmosphere 
without turbulence, but in the actual case we 
must introduce a weight function f (0, T) because 
the absorption of a spectral line is a function of 
the density and the temperature. Hence what 
is measured is 


Hi =[H.f(e, T)dolf fle, T)do (x) 


Let us assume that there is a general field H, 
of the order H, = 10 gauss. This is super- 
imposed by the turbulent field H,, so that in 
some parts of a granula the resulting field is 
H, + H;=+ 10 + 200 = + 310 gauss, whereas 
in other parts of it the field is H,-H,= 
+10-300=-290 gauss. The average from 
(1°) depends on the function f (0, T) and the 
way in which o and T are coupled with A. 
If f is systematically smaller for positive values 


of H, than for negative values, H' may very 


well be negative even if the real field H is 
positive, and its absolute value may differ by 
orders of magnitude).t Hence no conclusion 
can be drawn about the general magnetic field 
from the Zeeman-effect measurements and the 
results of these measurements are not relevant 
for our problem. 


* An example of this was discussed recently (ALFVÉN 
1952). It was shown that with a certain (rather specialized) 
coupling between H, o, and T a reversal of the field takes 
place. It is obvious that the reversal can be produced also 
under more general conditions. 


“= 


2. Field inside the sun 


The solar magnetic field is obviously gener- 
ated by some current system in the sun’s 
interior, probably in the solar core. It is very 
likely that the currents are produced by some 
magneto-hydrodynamic mechanism, but al- 
though much work has been done in the field 
by CowziNG (1934, 1945), ELSASSER (1946, 
1947, 1950, 1954), BULLARD (1948, 1949, 1954), 
PARKER (1954 a, b), FERRARO (1937), BIERMANN 
(1950), TAKEUCHI and Suimazu (1952) we 
are still far from a complete understanding of 
the origin of the sun’s—as well as the earth’s— 
general magnetic field. 

We can get information about the solar 
field by studying the way in which magneto- 
hydrodynamic disturbances are propagated. 

According to the magneto-hydrodynamic 
theory of sunspots (ALFVÉN, 1950, p. 116) 
magneto-hydrodynamic waves are generated 
in the central parts of the sun and travel out- 
wards along the magnetic lines of force. When 
they reach the solar surface they produce sun- 
spots. The theory has been worked out under 
the assumption that inside the sun the magnetic 
field can be approximated as a dipole field in a 
spherical shell between a central distance R, 
and the solar surface Ry and as a homogeneous 
magnetic field inside Ry. Recent results (ALE 
VEN, 1953) indicate that the field inside the sun 
also has a tangential component, but as the 
m.-h. velocity is proportional to the magnetic 
field, this component affects only the longi- 
tudinal displacement, and in the present paper 
only latitude displacements will be discussed. 
LEHNERT (1954, 1955) has shown that the 
Coriolis force modifies the m.-h. velocitv but 
for the relatively high frequencies associated 
with sunspot whirl rings the change is not very 
large. : 

When a wave produced by a disturbance in 
the central parts of the sun proceeds outwards, 
it reaches the solar surface at first at high 
latitudes and proceeds towards the equator. 
In this way the progression of the sunspot 
zones from high latitudes towards the equator 
is explained. 

It at a certain time the wave front coincides 
with the equatorial plane near the centre of the 
sun, the time T(x) when it reaches the latitude 
a of the solar surface is given by 


T(x) = AF(«) + BG(&) + const. 
Tellus VIII (1956), 1 
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where the first term derives from the motion 
in the dipole field and the second term from 
the homogeneous field (ALFVEN, 1945). A and 
B are constants, F(x) is a function which is 
cd paper (ALFVÉN, 1945, 
p. 11) and G(«) can be approximated as 


G(a) = sin? « 


If B =o, corresponding to-R, =0 (only 
dipole field), ES mallee alas for 
& = 90° so that the wave front reaches the 
surface at first at the pole and from there 
proceeds to the equator. However, if B is 
sufficiently large, corresponding to a sufficiently 
large value of Ry, T may have a minimum 
value at a latitude ei, somewhere between the 
pole and the equator. This means that the 
wave front reaches the solar surface at first at 
latitude amin from where it proceeds both 
towards the equator and towards the pole. 

In the cited paper (ALFVEN, 19451) it was 
found that the observational sunspot progres- 
sion curve was intermediate between two 
theoretical curves with 


= 2200 DB — 1087 ER. 70 
and 
A070 : B=.25 70, Ry = 05 


and perhaps more similar to the first one than 
to the second one. The two curves are shown in 


07010 20° 30° 40 50° 60° 70° 60° 90° 
latitude 


Fig. 1. Theoretical curve for the progression of the m.-h. 
wave front. (From Cosm. EI. p. 123.) 

10,2% B=108-10°; R,= 7.0 1010 

B=) 2 CHOR, > 0.5.1020 


De IT A222 


A= 0.41-Io -38. 


1 Atthe bottom of p. 15 in this paper there is a numeri- 
cal error. The value of F(o) should be 22.8 - 104? instead 
of 28.4 + 104°. This changes the values of A, B, a and Rj. 
The corrected values are used in the present paper. 
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Fig. ı together with the sunspot data. The 
first one has a maximum at about « = 33°, 
the second one at about « = 50°. The figure 
is the same as in the cited paper (with the 
difference that a small numerical error has 
been corrected and the high-latitude part has 
been included). 

In the earlier paper the high-latitude parts 
of the theoretical curves were not drawn 
because they were not supposed to be of any 
importance as no sunspots are observed above 
30° to 40°. It is of interest to note, however, 
that high latitude solar activity shows a tenden- 
cy to move polewards periodically. The pole- 
ward motion of solar prominences is well-known 
since long through the works by ABETTI (1938), 
d'AZAMBUJA (1941) and others. Recently a 
diagram published by ANANTHAKRISNAN (1952) 
shows very clearly that the prominence activity 
at high latitudes moves towards the pole in a 
way which is similar to the progression of the 
high-latitude branch of the m.-h. wave front 
according to Fig. 1. 

Hence there is a possibility that not only 
the sunspots but also the high-latitude prominence 
activity is caused by the same train of m.-h. wave. 
When the waves reach the solar surface they 
produce electromagnetic disturbances, which 
for some reason manifest themselves as sun- 
spots at low latitudes but as prominence 
activity at high latitudes. We are going to 
discuss this possibility. 

Hence we want to investigate whether both 
the observational sunspot progression curve 
and the maximum prominence activity in 
Ananthakrisnan’s diagram could be represented 
by the same dheoredcal curve and then to 
calculate the parameters A and B of this curve. 
The observational sunspot progression curves? 
have been calculated for the cycles 11—16 
and Ananthakrisnan’s diagram covers the cycles 
14—18. Hence the cycles 14, 15 and 16 are 
suited for our analysis. 

From the sunspot progression curves we 
take the time T,, when the curves pass the 
latitude « = 30° and the time T, when it 
passes the equator. On Re entbalrisnen’ s dia- 
Bins the best defined part is the time of arrival 

To of the activity wave to the pole. There are 
also some fairly well defined maxima at 55°, 


1 See Cosmical Electrodynamics p. 120—122. Fig. 5.6 
and 5.7: 
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Table 1 4 
Fe. WESEN A ON re En EEE a re En m —— en nn 
To T 30 55 Too A ARE 
Cycle ee PC RTE 
N S N S N S N S N S) N S 
14 08.7 09.3 04.1 7, 06.1 07.2 2.6 Zell 
15 20.3 20.7 15.7 16.1 14.9 15.0 17.6 18.0 5.4 5.7 am, DR 
16 35,2 31.0 26.6 26.4 23.5 24.1 20778270) TT 6.9 4.5 4.0 


but these are sometimes double. We design the 
time of these maxima by T,,. Table 1 gives 
the values of T, and T3) from the diagrams 
in Cosmical Electrodynamics p. 120—121 and 
T:s and Ty, from Ananthakrisnan’s paper. 
The table also gives the differences 


IN — Tas 
Re 0” Too 


The dilerence A, = T = Te, is always 
N3°= 4.6 years because the observational 
points are fitted to the average observational 
-curve (fig. 5.7 in Cosmical Electrodynamics). 
‘The maximum T;, for cycle 14 is not included 
in Ananthakrisnan’s diagram. 


The meanvof AR is 
NS OS VEAL. 


"We can now determine the parameters A 
cand B from the equations 


EN Oa Brei area: 
+ B [sin? 0° - sin? 30°] 
EL (0)= Too) = AF (0) Foo) + 
+ B [sin? 0° -sin? 90°] 
or 
Ay = 4 
A = 4 


[F (0) —F(30)] -0.25 B 
[F (0) -F (90)] -B 


In the earlier paper two alternative hypoth- 
ses were worked out: either are the sun- 
spots produced when the m.-h. waves reach 
the solar surface Ry = 7.0 101% or when they 
reach the surface Ry = 6.5 + 101°, In the first 
case we have 7°F(0) = 12.2: 1040; 7-°F(30) = 
= —22.0° 104; 7°F (90) = —65.8 : ve which 
gives AO TON ED ST > a 
= 2.2 ° 10% gauss : Cm°: H,= 12 oat here 
ais the dipole moment and H, the field at the 
pole. 


Mean 3.1 years 


In the second case we have 55F(o) = 22.2: 
+ 1040: &5F(30) = 28.6. 1029 7&5F(90) =77.0- 

104, which gives A = 0.46: 107%; B= 
= 37-107; a = 3.8 - 10°° gauss : cm®; H, = 22 
gauss. 

The curves corresponding to these two 
alternatives are shown in Fig. 2. They are not 
very different and both agree with the observa- 
tional sunspot progression curve within an 
error of a few months. 

The reason why the curve for Rg = 6.5 : 101° 
was introduced into the sunspot theory was 
that it is possible that the m.-h. waves need 
only reach a certain distance below the solar 
surface in order to produce a sunspot, which 
may have its “roots” at a certain depth. It is 
more likely that prominences which are 
observed above the photosphere should be 
caused when the disturbance reaches or passes 
the photosphere itself. Hence it is reasonable 
to choose the curve for Ry = 7.0: 101°, 

Fig. 3 consists of Ananthaleneagn: s diagram 
into which the observational dots of the sun- 
spot progression curves are introduced. They 
are taken from fig. 5.6 in Cosmical Electro- 
dynamics. The oe curve R,=7.0° 101°; 
A=0.83 + 10-38; B—55 - 107, is de through 


Atimd | | | LA | 


| > 

0° 10° 20° 30° 40° 50° 60° 70° 80° 90° 
latitude 

Fig. 2. Theoretical curves for the progression of the m.-h. 


waves 
A = 0.83 :10.°% B= 55-107; Ry = 7.0- 1010 
Sis =. A 0462105"; Bi 3728, rot Rn CS roi 
Dots are taken from Cosm. EI. Fig. 5,7. 
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Fig. 3. Ananthakrisnan’s diagram of solar prominences. The sunspot data from Cosm. EI. Fig. 5, 6 
are introduced. The curves are the full curve of Fig. 2, fitted to the sunspot dots. 


these dots. It is seen that at high latitudes it 
represents the maximum prominence activity 
rather well, especially for the cycle 15. For the 
earlier part of cycle 14 the prominence data are 
lacking. In cycle 16 the maximum activity at 
so’—s5° comes about one year too early. 
Since this analysis was made, the sunspot 
data for cycle 17 have become available. The 
sunspot progression curve has been constructed 
in the same way as earlier and plotted in the 
diagram. It is evident from Fig. 3 that the curve 
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represents the prominence activity in a satis- 


factory way. 
In the cited paper (ALFVÉN, 1945) the solar 
magnetic field was determined to be between 


the limits 
37 gauss > H, > 9 gauss. 


The value which corresponds to the curves in 
Fig. 3 is 
H, = 12 gauss. 
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3. Interplanetary field 


In the absence of currents outside the sun, 
the field would be a dipole field. However, this 
field will be violently disturbed by the ejection 
of ionized gases from the sun. From solar 
observations we know that often there is an 
ejection of matter from prominences. Schuster, 
Chapman, and others have given good reasons 
for supposing that magnetic storms are pro- 
duced by beams of ionized gas, which are 
ejected from the sun with velocities of the 
order of 108 cm/sec. and which go out to the 
earth’s orbit and beyond it. Such beams cannot 
pass through a magnetic field without disturb- 
ing it violently, and if the solar field at a certain 
moment were a dipole field the first beam shot 
out from the sun would change its character 
drastically. In the region where the beam passes, 
the lines of force will be transported outwards 
with the beam. At some distance from the 
sun the field may derive to a larger extent from 
currents in space than from currents inside the 
sun. Hence the term “the interplanetary magnetic 
field” seems to be more appropriate than “the 


solar field”. 


The outward transport of field lines cannot 
go on indefinitely without a corresponding 
inward transport. This can be assumed to take 
place in the sections of the sun’s field where for 
the moment no beam is ejected. In average over 
a long time the inward transport of magnetic 
flux must equal the outward transport. 


The ejection of beams takes place in the 
sunspot zones on both sides of the equator. It 
is essentially limited to a belt between + 30° 
and — 30° latitude. Near the poles the solar 
field is much less disturbed and is probably not 
very far from a dipole field. 


It is obviously difficult to construct a model 
of a field which is so complicated and so 
rapidly changing as the interplanetary field. 
Still it is very important to have a picture which 
gives some of the essential properties. The 
model, which we are going to propose, re- 
presents an ideal average field, a stationary 
state which always is considerably disturbed 
and probably has some sort of turbulence 
superimposed on itself. The model is a kine- 
matic one. The problem of the dynamics of 
interplanetary space is very complicated and 
cannot be analysed in a satisfactory way at 
present. 


u 


In order to make the model as simple as 
possible we assume that the field is always 
meridional, ie. the tangential component 
jel, = ©, : 

Further we assume that the field consists of 

four different parts limited by latitude angles 
Po and 9, and longitude angles A, and A, +/\A, 
the numerical values of which are discussed 
later. 
1. In the polar regions p > po and P<-% 
the field is not very much disturbed. We 
approximate it to a dipole field. Hence if the 
dipole moment is a, we have 


Rs 3 . 
Hr = H,(—) sing (1) 
R 
3 
He = 4, (+) cos ® (2) 
with 
I a 


The numerical values are 
I 
= H, = H.=6 gauss 


= x 33 3 
a = 2.1" 10% gaus cm (4) 


The total flux is: 


n'2 
®= [ Hr22Rcos gdp =2 R* H, cos? yy (5) 
Po 


or 
D = 1.9+ 1075 cos? p, gauss - cm? (6) 


2. In the equatorial region +9, >P >-9ı 
the field is disturbed by beams emitted from 
the sun. These beams occupy one or several 
sectors with longitudes À, < 1< 4, + /\2 
where /\1 is about 30°. Each part of a beam is 
supposed to move outwards radially with con- 
stant velocity. Because of the high conductivity 
the magnetic field components perpendicular 
to the velocity are “frozen in”. As according to 
our assumption every volume element of the 
beam moves outwards radially at a constant 
velocity, a surface element in the equatorial 
plane expands as R-1 This means that inside 
the beams the field is 

db 


Ha, (7) 


where ay is a constant which we shall determine 
later. Besides H, there is a radial field Hr which 
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we shall discuss later. As the field is meridional 
the total flux should equal the flux 


D="RYH,AA (8) 

which passes the photosphere. Hence we obtain 
Rı 

= R?H, AA= | Hp AARAR (9) 
+ Ro 


where R, is the outer limit to which the beam 
moves with its lines of force still going to the 
sun. As R, > Ry we obtain approximately 


from (7), (8) and (9) 


and 


i 
PE RoR R iat) 


The numerical value of R, should be de- 
termined in such a way that the field in the 
neigbourhood of the earth’s orbit gets a value 
which is in agreement with observations. From 
the theory of the storm effect of cosmic 


radiation we find values from (H;);=7 u 
gauss to 70 u gauss for small and for big storm 
effects. These values are acceptable in the 
theory of magnetic storms and aurorae (ALF- 
VEN, 1954; 1955). For our model we choose 
an intermediate value of 


(HE); = 30 u gauss (12) 


This gives 
Ra = 6.510!" cm 
and 
RE TOR SION 4;S 2100 
ur R R 
R 
= 0 I 
Gun (13) 


The value of R, is about 4 times the earth’s 
orbital radius. This means that out to this 
limit the lines of force in the beam go back to 
the sun. The conditions outside R, shall not be 
discussed here. 
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3. In the equatorial region + 9, >P>-9, 
outside the beam there is a field which falls as 


aq 


(14) 


where a, is a constant. As in this region the 
field which once has been transported out- 
wards by a beam is tending to go back to a 
dipole field, the exponent n should be some- 
where between the dipole exponent 3 and the 
beam exponent 1: 


Ten ei (15) 
In reality the exponent n may vary with r, but 
in our model we assume it to be constant. 
BRUNBERG and DATTNER (1954) have shown 
that the absence of an increase in cosmic 
radiation before a magnetic storm requires 
that near the earth the field is at least 10 u gauss. 
They derive an average value of n = 2.6, but 
this is under the assumption that the value of a 
is higher than we have assumed. With our 
value n goes down somewhat, and we are 
still in a reasonable agreement with their results 
if for the case of simplicity we put 


n=2 


(15?) 


The same flux condition as above gives 
Rı 
Su (2-4) = | Hi(en- AA) RAR 
Ro 


where rather arbitrarily we have integrated 
between the same limits as for the beam. This 
gives 

a 
7 R, ln Ril Rs 


Ag 
and 


R 2 
Hy = 0.88 (+) gauss 


4. In order to satisfy the flux condition there 
must be regions p > y > pand-ps <p< 
<-9, between the dipole regions and the 
equatorial region. In these regions the field is 
essentially radial. As flux from the interior 
dipole field is penetrating the solar surface, 
there must also be a radial field in the equatorial 
regions. 
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a) — solar surface. The field outside the equatorial 
| = regions is of no importance in the theory of 
ZZ EE magnetic storms and aurorae. In the cosmic- 
ray theory the discussion was essentially con- 
fined to the conditions near the equatorial 
plane. The introduction of the present model 
makes it possible to extend the analysis also 
to regions far from the equatorial plane. This 
may call for modifications of some details in 

the cosmic-ray paper. 
The beam is associated with an electric field 


Ey == H 


= 


itt Fe 
RE RAN: in the beam, so that there is a voltage difference 


V, over the beam 


N ha 
RoRı 


Vite RAAE,=— 


= const. 


Outside the beam (in region 3) there is a field 


associated with a voltage difference 


va (2%-A Aa 1 


V, = (2-AA) RE, 


Sie Ren 


Fig. 5. The solar magnetic field outside a beam. 


Fig. 4 and Fig. s illustrate the model. The 
limiting angles have been chosen as gy = 30°; 
9, = 20°. In the pictures the difference between 
the regions 2 and 4 is probably exaggerated. 
The beams go out with great velocity, and the 
flux they take with them has to be connected 
to the sun by the region 4. On the other hand, 
region 3 is characterized by the resetting of the 
field and the transition to 4 should be less 
pronounced. 

In the equatorial regions (2 and 3) the model 
is the same as used in the theory of cosmic 
storm variation and local generation of cosmic 
radiation and in the theory of magnetic storms 
and aurorae (ALFVEN, 1954; 1955). The model 
field in this region is shown in Fig. 6. 

2 ie on = Fig. 6. Field model in the equatorial plane. (Tellus: 
g g near the Origin of Cosm. Radiation, 6, p. 238). 
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Table 2 
cn ee es 4 à À © 


Magnetic field Hy 


Region 


Polar Region | Dipole field 3 (6 gauss) 
Equatorial Inside beam I 0.065 gauss 
un Outside beam 2 0.88 gauss 
In our stationary model we have 
Va > V, 
À In(R,/R 
v=v TAN ( 1/ 0) R 


TE R, 


The compression of the gas in the region 3 
may in part be compensated by a flow parallel 
to the lines of force. 

The electric conditions in region 4 are more 
complicated. In fact they cannot be described 
in a satisfactory way in a stationary model, 
which shows one of the limitations of our 
model. In reality the conditions in region 4 are 
non-stationary and depend on the life time of 
the beam etc. 

The driving mechanism of the ejection of 
beams is probably the intense heating of some 
parts of the solar atmosphere by solar flares 
and other kinds of solar activity. The heating 
makes the gas expand and this weakens the 
magnetic field. The very weak magnetic field 
in the beam near the solar surface is reconcilable 
with this process. 

Table 2 gives the field in different parts of our 
model. 

The values of H, refer to a point near the 
limit of the transition region 4. In the equatorial 
plane we have H, = 0. 

The shape of the lines of force is given in 
the following way. A line of force which 
passes the solar surface at a latitude 9, intersects 
the equatorial plane at a point at the distance 
R. The flux ©, between the equator and the 
latitude is 


10, = D, sin? y 


where 
D = ake, 


is the total flux through the solar surface. 
Inside the beam the flux 8% between the 
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Radial field Hy 


| N | Sun’s equator | Earth’s orbit | Sun’s surface | Earth’s orbit 


(0.6 u gauss) 10 gauss — 
30 u gauss 2 gauss 50 u gauss 
19 u gauss 2 gauss 15 u gauss 


equator and a circle in the equatorial plane 
with the radius R, is 


R-R, 


Dh = @ 


which gives for the two ends of the line of 
force ; 


R 1 ; : 
—=1I+(-1})sin®®=1I+ 929sin? 


The same relations for the region outside the 
beam are 


Selo IR; 
FREE 
and 
In (R/R,) = (à) sin? @ = 4.5 sin? y 
0 


The following table gives the intersections 
with the equatorial plane of the lines of force 
passing the solar surface at the latitude 9. 


Table 3 
Distance in solar radii 
Latitude 
(0) inside the outside the 
beam beam 
o I I 
5 8.1 1.05 
10 29 1.23 
15 63 1.58 
20 IIo 2.23 
25 167 3.38 
30 233 5.50 
39 9.3 
40 17.0 


4. Consequences of the solar model field 


Near the equatorial plane the present model 
is essentially the same as used in the papers 
about the cosmic radiation and the magnetic 
storms and aurorae. This means that the model 
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gives electric fields which are needed in order 
to explain the magnetic storms and aurorae, 
and it also explains the variation of cosmic 
radiation. 

Further the model gives the progression of 
the sunspot zones ads the equator and the 
poleward motion of the prominences as a 
simple consequence of the motion of m.-h. 
waves in the sun. 

The model seems to be in striking confliet 
with the Zeeman-effect measurements if these 
are interpreted in the naive way as really 
meaning an average field. It is possible that a 
magneto-hydrodynamic treatment of the tur- 
bulence in the photosphere will reveal that 
no real conflict exists. 

The cut-off in the cosmic-ray spectrum has 
from time to time been quoted as an argument 
for or against a solar magnetic field. The 
present model is so far from a dipole field that 
the usual arguments are invalid. It may be 
that the cut-off is related to the solar magnetic 
field but in the way which has been discussed 
in one of the cited papers. 


An argument which has been used against 
the assumption of a solar field is that ch 
flares often are accompanied by increases in 
cosmic radiation which follow within an hour 
after the flare (FORBUSH, STINCHCOMB and 
SCHEIN, 1950). This phenomenon has been 
studied, especially by Simpson’s group (Simp- 
son, Frror and TREIMAN, 1954), with the neu- 
tron counter technique, which is sensitive in 
the rigidity region of some times 107 gauss : cm. 
If the solar field were ‘an ideal dipole field, 
cosmic rays produced near the solar flare could 
not reach the earth, unless the dipole moment 
were very small (corresponding to a surface 
field of 10-4 gauss). The time lag is so small 
that particles generated near the sun cannot 
have time to diffuse along a very complicated 
orbit out to the earth. 

The present model removes this difficulty 
because the magnetic lines of force run out 
radially from the equatorial region. Table 3 
shows that in the beam there are lines reaching 
the neighbourhood of the earth (R/R, = 200) 
coming from @ = 25°—30°. Only relatively 


Fig. 7. Photograph of the solar corona. (Wallenquist 1954.) 
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small disturbances are needed in order to 
shift this point. Outside the beam the condi- 
tions are not so favourable, but even there a 
reasonable drift perpendicular to the lines of 
force, or a diffusion, will allow cosmic-ray 
particles to travel almost directly from the 
flare out to the earth. Usually they should 
reach the earth from the sun’s direction, but 
due to the complicated structure of the magnet- 
ic field they may also come in from other 
directions. This explains that many of the 
flares produce increases in Firor’s zones (Frror, 
1954) but that there also are many noteworthy 
exceptions, when flare effects have been observ- 
ed far away from the expected zones. 


5. Relation to the structure of the corona 


If we compare Fig. 4 or especially Fig. 5 
with a photograph of the corona during a 
sunspot minimum, we find some resemblance: 
The fine structure of the corona is similar to 
the magnetic lines of force of the model. The 
reasons for this may be the following. 

The proposed model gives a simplified 
picture of the real conditions. There may be 
a superimposed turbulence, and especially near 
the sun’s surface the model field is usually 
much disturbed by local fields from sunspots, 
etc. It is well known that the solar corona is 
very much influenced by sunspots and prom- 
inences (see e.g. van de Hulst). The corona 
should also be controlled to some extent by 
the general field. We whould expect that the 
general field is most important for the structure 
of the corona at times of sunspot minimum, 
when the influence of local solar activity is 
small, and at large distance from the centres of 
disturbance. Hence a comparison between the 


field model and the corona photographs should 
preferably be made for a sunspot minimum and 
refer especially to the outer corona. 

This indicates that the similarity between 
Fig. 5 and a minimum corona may be signifi- 
cant. The “brushes” or rays near the poles are 
usually interpreted as given by the magnetic 
field, and this seems to be correct. The most 
prominent phenomena at lower latitudes are 
the “equatorial streamers” or “fans” which 
have a mainly radial fine structure. As long as 
it was firmly believed that the solar field was 
a dipole field it was impossible to identify 
this fine structure as given by the magnetic 
field lines. With the present model it seems 
plausible that practically all the fine structure of 
the corona is given by the magnetic field.“ This may 
be true also during a sunspot maximum but 
then the magnetic field is influenced very 
much by local fields from sunspots and prom- 
inences. Our model of the magnetic field is 
perhaps too simple to account for the condi- 
tions near the sun at the time of a maximum. 
At large distances where the local disturbances 
are less pronounced it may be hoped that the 
model is a better approximation even during 
a maximum. 

If we accept the identification of the fine 
structure of the corona as magnetic field lines, 
we will have a possibility to make a better 
model of the general magnetic field than the 
present model. The field is not necessarily 
symmetrical with respect to the equatorial 
plane; the beams may be shot out at any angles 
and are often much thinner. 


1 Note added in proof: At the conference in 
Guanajuato T. Gorp presented independent arguments 
for the same interpretation. 
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The 24-Hrs Averages of Cosmic Ray Intensity in Stockholm 
1947—48 and 1954 


By ARNE ELD SANDSTROM, Fysiska Institutionen, University of Uppsala 


(Manuscript received May 25, 1955) 


Abstract 


Curves are given showing the variation of the 24-hrs averages of the C.R. intensity as measured 
by means of directional counter telescopes during two registration periods separated by six 
years. During both periods the maximum deviations were 2—3 per cent. In accordance with 
the view expressed by Sımpson (1954) the periods of maximum are to be regarded as recurrent 
rather than being a harmonic function of time. There are some traces of a 27-day component 
having an amplitude of + 0.4 per cent. The C.R. intensity as a function of time is compared 
to the solar cycles, the relative sun-spot number, the sun-spot area and the Ky-index. 


The cosmic ray registrations in Stockholm 
Oct. 1947 to May 1950 (MALMFORS 1949, 
SANDSTROM 1955) offer possibilities for a 
study of the 24-hrs averages as a function of 
time. That as yet no such study has been 
made is partly due to the fact that these 
averages are affected by some systematic in- 
fluence, easily eliminated as regards the diurnal 
variations, but becoming serious as soon as 
absolute intensities are concerned. For part 
of the registrations, however, these difficulties 
can be overcome. 

The 24-hrs averages were derived from the 
sum of twelve bihourly values corrected for 
the barometer effect (LINDHOLM 1950). As a 
function of time these averages display a 
trend due to a systematic change of counting 
rate. After correcting for this trend the varia- 
tions of the 24-hrs values were computed 
in per cent of the average intensity. 

The intensity was registered in the N-, Z-, 
and S-directions by means of four two-fold 
coincidence channels in each direction (MALM- 
FORS 1949). Rather frequently it happened, 
however, that one or several channels were 
Tellus VIII (1956), 1 


out of order. Owing to the individual charac- 
teristics of the counters the counting rate 
varies with the combination of channels. As a 
result the present study had to be confined to 
the period Oct. 1, ’47 to July 10, ’48. Only 
4-channel registrations were considered except 
as regards the N-direction from March 27 
until May 20 1948 for which interval com- 
paratively homogeneous 3-channel registra- 
tions are available. The 24-hrs average is 
approximately 10° impulses derived from a 
total of 12 x10. The corresponding statisti- 
cal fluctuations are of the order of magnitude 
of 0.1 per cent. 

The curves in Fig. 1 representing the three 
directions agree even as to detail. A slight 
difference ought to be expected as the diurnal 
variation has a phase difference of 2.5 hours 
between the N- and S-directions, and 1.8 
hours between the Z- and S-directions (BRUN- 
BERG 1953, BRUNBERG and DATTNER 1954). 
The influence on the 24-hrs values is obviously 
too small to be recognizable. 

A certain periodicity related to the 27-day 
cycles is easily perceived. Such a variation 
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Fig. 1. The three uppermost curves represent the variation of the daily averages in the N-, Z-, and S-direct- 

ions. The fourth curve represents the daily sums of the Kp-indices and the fifth curve the daily maximum 

values of these indices. Below is the relative sun-spot number and (from Jan. 1, 1948) the sun-spot area as 
functions of time. The vertical lines divide the diagram into 27-day cycles. 


has often been reported (HEISENBERG 1953, 
p- 39, SIMPSON, Frror, FONGER, and TREIMAN 
1953, ELLIOTT 1953, EHMERT 1953) but prob- 
ably the present issue is the first time it was 
observed by simultaneous measurements in 
three directions. Through the experimental 
studies of BRUNBERG (1953) it is now possible 
to determine the directional coordinates of 
the incident primary particles. The angle Py 
between the direction of the primaries and the 
equatorial plane did not differ very much 
for the three sets of telescopes. This is illustrat- 
ed in Fig. 2, where P represents the point of 
observation. The north-south aperture of the 
telescopes was big enough to make all the 
three sets of telescopes measure practically 


the same beam of primaries integrated over 
24 hours. 

A comparison between the three curves 
offers some means of control as to the reliabi- 
lity of the registrations. As an example the 
N-direction differs from the other two direct- 
ions on Oct. 15 to 19. The close correspon- 
dence elsewhere indicates that something was 
wrong with the N-telescopes and, as a matter 
of fact, the daily averages were calculated 
from less than 12 bihourly values. Incidentally, 
these days were never included in the studies 
of the diurnal variation. 

After having been adjourned for more than 
three years registrations were resumed in 
1954. The new telescopes differ from the 
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Equatorial 


Fig. 2. The directions in the geographical meridian of 
the beams registered by the three sets of counter 
telescopes. 


former set only as regards aperture and 
counting rate. They do not display any trend 
of the kind that is so disturbing in the ’47—’48 
measurements. Therefore a comparison be- 
tween the two intervals is doubly valuable. 
The effect of the changed aperture can be 
seen from Fig. Days with one or more 
channels out of order were excluded. 

As expected the three directions give 
closely corresponding records. (Fig. 3). We 
find the same variations related to the 27-day 
cycles as in the former case except during 
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Oct. and Nov. They reappear at the end of 
the period of registration. 

The variations average + 0.8 per cent for 
Oct. 47 to July ’48 as well as for April to 
Dec. "54, although the first interval coincides 
with a period of comparatively high solar 
activity and the second with one of minimum 
activity. The maximum deviations from the 
average intensity is between 2 and 3 per cent. 

Although in general one maximum and 
one minimum of C. R. intensity are related 
to each solar cycle a phase shift apparently 
exists. Sometimes the registrations display a 
distribution of maximum and minimum re- 
gions even suggesting a period of twice the 
Tength of the 27- day cycle. There is really no 
reason for the existence of such a period and a 
scrutiny of Fig. 1 reveals that the position of 
the intensity maximum moves through the 

27-day cycles in a comparatively continuous 
Way. Incidentally SIMPSON (1954) prefers to 
speak of recurring periods of maximum in- 
tensity instead of associating them with the 
27-day cycles. 

Comparing all three directions the averages 
of the 24-hrs values of corresponding days in 
all the available 27-day cycles were calculated. 
Finally the first harmonic was derived and 
found to have an amplitude of + 0.4 per cent. 
This is little as compared to the superimposed 
recurring variations. 

The diffusion of particles outwards from 
the sun (ALFVEN 1954) ought to display varia- 
tions related to the solar cycles. The occasional 
corpuscular radiation, originating from disturb- 
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Fig. 3. The variations of the 24-hrs averages in Stockholm during 1954. 
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ances in the sun, will cause irregular changes 
of the C.R. intensity partly depending on 
the position of the earth relative to the emitted 
beam. However, if the latter lasts for more 
than one solar cycle it will cause a 27-day 
periodicity (BRUNBERG and DATTNER 1954). 
It follows that quite regular variations are 
not to be expected. The nature of those 
observed appears to be in accordance with 
this theory. Apparently the + 0.4 per cent 
variation is a characteristic of the diffusion 
of particles from the sun outwards, while the 
less regular but much larger variations are 
due to the action of direct beams. The inter- 
mixture of the two effects makes any calcula- 
tion of the true periodic variation difficult. 
In addition, very probably both the amplitude 
and the phase of the latter vary with condi- 
tions on the sun. 

It is therefore especially interesting to com- 
pare the C. R. intensities with the variations 
in relative sun-spot number (Fig. 1). Another 
noteworthy fact is that the regions of high 
relative sun-spot numbers in Oct. ’47 and in 
April and May ’48 seem to be related to 
regions of low C.R. intensity. To try to 
advance one step further an attempt was 
made to correlate the C. R. variations with 
the total sun-spot area. Data on the latter 
are available from Jan. ’48 (OHMAN 1950). 
We find that the period of low C. R. inten- 
sity in April and May ’48 coincides with a 
sequence of days with large sun-spot areas. 
However, for the N-direction the correlation 
coefficient is only 0.22 (not significant). 

As regards 1954 there are as yet no data 
available either on relative sun-spot numbers 
or on sun-spot areas. Data on geomagnetic 
character figures are also lackirfg. Once again 
we have to resort to the period Oct. ’47 to 
June ’48. One of the curves in Fig. 1 represents 
the sums of all the eight K,-indices of each 
day (HOwE and Weisman 1949), while the 
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one below shows the daily maximum of the 
same indices. The differences between these 
two curves are very small, the regions of 
maximum values being more marked as 
regards the K,-sums. As was mentioned in a 
previous paper (SANDSTRÖM 1955) it does 
not matter, with regard to the present degree 
of accuracy, which of these two values is 
selected for classifying the days geomagnet- 
ically. Comparing the curves we find that, 
even if a cosmic ray minimum does not 
exactly coincide with a maximum of geomag- 
netic activity, the former succeeds the latter 
in a way indicating a common origin. In this 
case it is impossible, however, to recognize 
the rule given by Simpson (1954) that the 
geomagnetic maximum appears two days 
later than the end of the cosmic ray maximum. 
This discrepancy may be due to the fact that 
the present material is statistically small. 

The change from day to day of the 24-hrs 
average constitutes a trend which certainly 
affects the observed time of maximum of 
the diurnal variation. However, the phase 
difference between the C.R. intensity and 
the geomagnetic variations has a remarkable E 
effect. Given a sufficiently large number of ! 

À 


days of the same geomagnetic activity, posi- ' 
tive and negative trends will be randomly 
distributed (SANDsTROM 1955) even as re- 
gards magnitudes. As a result the average daily } 
trend will become insignificant as soon as ! 
the number of days is sufficiently large. The 
periodicity of the general diffusion of particles- 
explains the fact that negative as well as 
positive trends are to be found even on 7 
magnetically quiet days. 

The present paper is part of a programme 
sponsored by the Swedish Natural Science 
Research Council. 

The author wishes to: express his thanks to 
Miss Ulla Björkegren for all the numerical 
calculations. 
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Experimental Evidence of the Directional Components of 


Cosmic Rays 


By ARNE ELD SANDSTROM, Fysiska Institutionen, University of Uppsala 


(Manuscript received June 13, 1955) 


Abstract 


It is possible to resolve the diurnal variation into two components. The first of these has 
its maximum at 18h and corresponds to the particle wind predicted in the theory by ALFVÉN 
(1954). The amplitude ought to be constant both as a function of geomagnetic activity and as 
a function of time. The second component has its maximum at 12h. It represents the sum of 
the various kinds of radial particle flow, and it ought to display variations indirectly related to 
geomagnetic activity. A search for experimental evidence of these components is made possible 
through the determination of the incidental directions of the primary particles by BRUNBERG 
(1953) and BRUNBERG and DATTNER (1953). An analysis of C.R. measurements from Stockholm, 
Nagoya, and Huancayo seems to confirm the two components and their expected properties. 
A very doubtful amplitude change of the 18h-component during major magnetic storms is 
not at variance with theory. The Stockholm measurements were treated by SANDSTROM (1955) 
as to geomagnetic effects, those from Nagoya and Huancayo by SEKIDO and KODAMA (1952). 
A seasonal phase shift reported by the latter for geomagnetically quiet days is compared to 
measurements in Stockholm. If it exists at all the reported phase shift is evidently due to local 
effects. Other secular or seasonal phase shifts (ELLIOT and THAMBYAPHILLAI 1953) have not been 
treated in the present study, as the data did not provide any example of these effects. Several 
points are open to discussion because of the data’s being collected with instruments of widely 
different properties. Directional counter telescopes offer the most valuable results, if they are 


designed for a comparatively small aperture and measure in at least two directions. 


Introduction 


According to a theory by ALFVÉN (1954) 
the main part of all cosmic radiation originates 
in the sun. The following arguments suffice as 
an introduction to the study of available ex- 
perimental results. 

The wellknown daily variation can be re- 
solved into components—-(1) the wind of 
cosmic particles caused by the rotation of the 
sun overtaking the earth in its orbit and (2) 
a radial flow of particles from the sun. The 
latter in its turn consists of two components, 
one being a general outward diffusion of 
particles and the other a direct flow of particles 
inside the corpuscular beams associated with 


active sunspots and solar flares. Outside the 
beams there ought to be a particle flow of 
very low intensity in the opposite direction. 
These conditions are illustrated by the simple 
scheme in High The corpuscular flow is ver 
variable (even changing its sign) while the 
particle wind and the particle diffusion ought 
to display fairly constant features. It is quite 
possible, however, that the nonvariable part- 
icle wind is disturbed whenever a corpuscular 
beam sweeps across the earth. BRUNBERG and 
DATTNER (1954) have shown that the rotating 
eam can retard as well as accelerate the cosmic 
ray particles in the space surrounding the sun, 
the effect extending far beyond the earth’s 
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Fig. 1. Simple scheme showing the main components 


of the cosmic radiation according to Alfvén (1954). 


orbit. They have also pointed out that the 
direct beam will give rise to a 27-day peri- 
odicity in case it lasts for a sufficiently long 
time. 


Corrections for the Geomagnetic Field 


The directions of the incident particles have 
to be computed from the time of maximum 
and the direction of measurement at the point 
of observation. The correction to be applied to 
the observed time of maximum corresponds 
to a simple coordinate transformation in 
the harmonic dials. It is done by turning 
the p- and q-axes through an angle which 
can be read from the experimental curves of 
BRUNBERG (1953) and of BRUNBERG and DATT- 
NER (1953). The average momentum of the 
primary C.R. particles is assumed to be 
20GeV/c (BRUNBERG and DATTNER 1954). 
Fig. 2 illustrates this correction graphically. 


Fig. 2. Graphical correction for the east-west deviat- 
ion of the C. R. particles in the geomagnetic field. 
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The Directional Components in the Harmonic 
Dial 

When the geomagnetic effect has been 
corrected, the new harmonic coefficients p, 
and q, represent the two main components 
discussed in the introduction. The vector p, 
at 18h evidently represents the wind of C.R. 
particles, while q, at 12h is the sum of the 
diffusion of particles and the corpuscular flow 
(Fig. 3). 
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Particle diffusion 
| 


| Corpuscular flow 
| inside beam 


42h 


The vectors representing the components of 
cosmic radiation. 


Fig 


According to theory we might expect pı 
to be constant under all conditions. During 
magnetically quiet days q, ought to be negli- 
gible or positive (equalling the difference 
between a flow inwards and the diffusion 
outwards). As already mentioned, disturbances 
are to be expected. To decrease the influence 
of the statistical fluctuations it is always 
necessary to average over a number of days. 
In this way irregular effects will be diminished 
or even eliminated if they are of alternating 
sign. 

The secular change of the time of maximum 
(ELLIOT and THAMBYAPHILLAI 1953, BRUN- 
BERG 1954) is too large to be accounted for 
solely by changes in qj. It is certainly inde- 
pendent of the daily geomagnetic effects 
(SANDSTROM 1955) and therefore it cannot 
easily be attributed only to changes in the 
corpuscular flow. 


Measurements from Stockholm 


The C.R. data collected in Stockholm by 
Matmrors (1949) Oct. 1947—June 1950 have 
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Fig. 4. The 18h-(p,-) and 12h-(4,-) components in 

Stockholm Oct. ’49—June ’so in their relation to the 

maximum Kp-index. The broken line in the upper 

diagram corresponds to the average of the first three 
points. 


been treated by SANDsTROM (1955) as to geo- 
magnetic effects. The days were divided into 
five groups according to their maximum K,- 
indices. Measurements are available in three 
directions, zenith, 30° to the north, and 30° to 
the south. After having been corrected for the 
deviation of the primary particles in the earth’s 
magnetic field the three directions ought to 
display closely corresponding variations. Ac- 
cordingly in the present study we employ 
the averages of all three, thus diminishing 
the statistical fluctuations by a factor 0.6. 
The limits of errors will still be very big for 
the magnetically quiet days comprising class I 
as well as the most disturbed days comprising 
class V. 

From Fig. 4 we find that p, can be regarded 
as constant for classes I--II, that is the wind 
of primaries remains constant when the K,- 
index increases to 5+. Contrary to theory p, 
appears to change its sign above this K,- 
value. However, when the statistical fluctua- 
tions are taken into account, this may be found 
to be accidental. 


It is very interesting to note the changes of q,. 
Beyond any doubt the amplitude of the 12h- 
component increases with increasing geomag- 
netic activity. The general trend is consistent 
with the assumption that this component can 
be regarded as negligible on magnetically quiet 
days. 


Measurements from Nagoya 


SEKIDO and Kopama (1952) have published 
a study of the relations between the cosmic 
ray diurnal variation in Nagoya 1949—SI 
and the magnetic character figures. After 
correcting for the deviation of the primaries in 
the geomagnetic field it is possible to make a 
comparison between their figure 4 and the 
data from Stockholm. Instead of the K,- 
indices Sekido and Kodama refer their results 
to the Japanese character figure (J.C.F.) dividing 
the days into five groups. To diminish the 
statistical fluctuations they finally work their 
groups A and B into one. Also the fifth group 
is wider than the remaining two. Accordingly 
the groups are not equally spaced. The differ- 
ences in referring to geomagnetic activity 
form no serious obstacle in a qualitative com- 
parison with the Stockholm figures. 


The Nagoya results with a counter telescope 
aperture of + 12° are illustrated in Fig. 5. 
Evidently the tendency is the same as in the 
diagrams of Fig. 4. A very prominent differ- 
ence is the fact that the absolute value of the 
amplitude is approximately three times that 
ever recorded in Stockholm. It is also very 
much bigger than that recorded in Nagoya 
with apertures of + 40° and + 85°. The 18h- 
and 12h-components as measured with these 
apertures are illustrated in Fig. 6. The harmonic 
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The 18h-(p,-) and 12h-(q,-) components in 


Fig. 5: 
Nagoya March ’49—Sept. ”5o in their relation to the 


Japanese magnetic character figure. Small aperture. 
The broken line in the upper diagram corresponds to 
the average of the first three points. 
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Fig. 6. The 18h-(p,-) and 12h-(4,-) components in 

Nagoya March ’49—Sept ’50 and Oct. ’47—Apr. ’49. 

Large apertures. The broken lines in the two upper 

diagrams correspond to the averages of the first four 
points. 


coeflicients were taken from Fig. 2 in the 
paper by Sekido and Kodama. In both cases 
the ı8b-component appears to be constant. 
In accordance with theory the amplitude of 
the 124-component increases with increasing 
geomagnetic activity. The discrepancies be- 
tween the results from telescopes with dif- 
fering apertures will be treated below. 


Measurements from Huancayo 


In their paper SEKIDO and KopaMa (1952) 
have also treated data from Huancayo for 
the years 1937—45. The position of Huancayo 
close to the magnetic equator makes it difficult 
to correct for the deviation of the primaries 
in the geomagnetic field (compare BRUNBERG 
and DATTNER 1953, Diagrams 35—37). We 
must also remenber that the registrations were 
made by means of an ionization chamber and 


% 


4/1/48 17 148 171/49 


20.4 Kp — 3? 


1/1749 


represent the intensity integrated over the 
half sphere. If the average momentum of the 
primaries were assumed to be 20 GeV/c, 
the angular correction for the deviation in 
the earth’s magnetic field would be 105°. 
However, the ionization chamber was pro- 
tected by shields equalling 12 cm of lead, 
which excluded most of the particles of low 
energy. It is impossible to predict the actual 
average value of momentum for the residual 
radiation. For 25 GeV/c the angular correction 
is roughly 80°. This has been selected as a 
first approximation. The harmonic coefficients 
were taken from Fig. 5 in the paper by S£- 
KIDO and KODAMA (1952). In this case they 
employed the international character figure C 
to classify the days geomagnetically. Unfor- 
tunately the relation between the C-figures 
and the K,-indices is complicated (Howe and 
WEISMAN 1949). 

The indirect relation of the two main com- 
ponents to geomagnetic activity is shown by 
Fig. 7. The 18h-component can be regarded as 


Fig. 7. The two C. R. components as functions of 

geomagnetic activity. From measurements as Huancayo 

1937—45. The averages have been calculated from 

the values cited by Sekıpo and KODAMA (1952). The 

broken line in the upper diagram corresponds to the 
average of all four points. 
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Fig. 8. The 18h-(p,-) and the 12h-(q,-) components from the Stockholm measurements versus time for days with 
a geomagnetic activity characterized by 2— Kp 3+and 4— Kp s+. The broken line in the upper two diagrams 
indicates the same average value as the corresponding line in Fig. 4. 
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Fig. 9. The 18h-(p,-) and 12h-(q,-) components from 
the Huancayo measurements versus time. 


constant while the amplitude of the r12h- 
component increases with the geomagnetic 
activity. 

Evidently, despite the uncertainty about the 
correction for the carth’s magnetic field, the 
two components seem to behave in the ex- 
pected way. 


The Components as Functions of Time 


The time dependence of the Stockholm 
measurements can be studied as regards classes 
II and III only. Days belonging to class I as 
well as classes IV and V are too few for such a 
study. Fig. 8 displays the variations of p, 
and q, in these cases. The 18-components is 
comparatively constant as regards days of 
class II but not those of class III. The amplitude 
of the 18*-component need not necessarily be 


constant, as fluctuations, though slight for 
days belonging to one and the same class 
(BRUNBERG and DATTNER 1954), are to be 
expected. 

As far as results from Nagoya are concerned, 
Sekido and Kodama have not published any 
data suitable for studies of the time dependence. 
They report only a seasonal change for mag- 
netically quiet days. This will be discussed 
below. 

As regards the registrations from Huancayo, 
Sekido and Kodama have divided the data 
into four periods. As they were also sorted 
in four groups according to geomagnetic 
activity, we get a good picture of the variations 
(Fig. 9). The diagram for p, does not show 
any systematic variations, q, on the other 
hand displays a marked systematic variation 
with the same trend regardless of geomagnetic 
activity. 


Seasonal Variations 


There exist observations at variance with 
theory and also some contradictory experi- 
mental results. 

A very early time of maximum during 
spring and summer has repeatedly been re- 
ported during the last few years. These reports 
are confirmed by registrations in three direc- 
tions in Stockholm (geomagn. lat. 58°) and 
in Kiruna (geomagn. lat. 68°) during 1954. 
As no K,-indices are as yet available for the 
year 1954 and the number of days is in any 
case too small for the purpose, no study 
based on this material of the phase as a function 
of geomagnetic activity has as yet been made. 


| 1 0.4 0.4 
wag Sp 49 Sp! 


Apr- June 49 


July-Sept 48 
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Fig. 10. The seasonal variation of the daily maximum 

for magnetically quiet days as reported by Sekido and 

Kodama (to the left) compared to the partly con- 

temporary registrations on magnetically quiet days in 

Stockholm (to the right). Corrections have been ap- 

plied for the deviation of the primaries by the earth’s 
magnetic field (compare text). 
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In their paper SEkIDO and Kopama (1952) 
published a harmonic dial demonstrating the 
variations of the daily maximum from winter 
(W) through spring (Sp), summer (S), and 
autumn (A) during 1949 and 1950. When the 
circles indicating standard deviations have 
been omitted and the deviation of the primary 
particles in the earth’s magnetic field has been 
corrected for, the harmonic dial turns into 
the left diagram of Fig. 10. Only magnetically 
quiet days are considered (K, < 2+). Incident- 
ally, the points winter ’49 and winter ’so are 
not uniquely defined. Therefore they have not 
been joined to the other points by lines in- 
dicating the sign of the phase shift. 

For comparison corresponding values from 
Stockholm Oct. ’47 to Aug. ’49 are illustrated 
by the harmonic dial to the right in the same 
figure. After the deviation of the primaries 
by the earth’s magnetic field had been corrected 
for, the averages of the three directions were 
calculated. The statistical fluctuations are great 
owing to the small number of days. Only days 
with K,< 2+ were regarded as sufficiently 
quiet. Nevertheless the Stockholm measure- 
ments cover a far smaller area than do the 
Nagoya results. It even happens that the two 
phase shifts have opposite signs. The magnitude 
of the phase shift is also far bigger as regards 
Nagoya than Stockholm. Considering all 
the discrepancies between these two diagrams 
it seems reasonable to suggest that the seasonal 
variation reported by Sekido and Kodama is 
due to some local effect. 


The Aperture of the Instruments 


The measurements from Stockholm were 
recorded with counter telescopes having an 
aperture of 60° x 120°, the angle in the mag- 
netic meridian plane being + 30°. 

Sekido and Kodama worked with three 
apertures: + 12°, + 40°, and + 85°. Their 
telescopes had quadratic sections. It is re- 
markable that the relationship between the 
12h-component and geomagnetic activity is 
less apparent as regards the two large apertures 
(Fig. 6). A possible explanation would be 
the superposition of daily variations with 


À 3 TT 
phase sbifts covering a part of the period > = 


The curves of Brunberg and Dattner show 
that there is a concentration of initial direct- 
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ions in the NS-plane for varying directions 
of the telescopes. As regards the EW-plane 
curves are lacking, but a considerable spread 
of initial directions is indicated. Obviously 


7 
the angular spread is > -. 
3 


The Components of the 125-Maximum 


It is very difficult to resolve q, into compo- 
nents. However, during magnetically quiet 
days it will certainly be due to particle diffu- 
sion, the backstream of particles outside the 
corpuscular beams being negligible in com- 
parison. To judge from the Stockholm mea- 
surements as well as those from Nagoya 
the latter component is small in itself. The big 
value of q, for magnetically heavily disturbed 
days will then almost entirely correspond to 
the flow of particles inside the beam from 
the sun. 

By analysing the registrations from neutron 
monitors SIMPSON (1954) has shown that the 
prominent periods of maximum intensity are 
recurrent and indirectly associated with in- 
creased geomagnetic activity. This is also 
borne out by registrations in Stockholm 
(SANDSTROM 1955). Thus it appears to be 
well established that a corpuscular flow from 
the sun constitutes a part of the 12"-component. 
The fluctuations registered by means of 
neutron monitors are bigger than those from 
measurements with ionization chambers and 
counter telescopes (SIMPSON 1954). Accord- 
ingly the soft component appears to abound 
in the corpuscular flow. 

From counter telescope measurements SAND- 
STRÖM (1955) has found that there may be a 
small 27-day variation underlying the large 
recurrent increases of intensity. Apparently 
this variation (& + 0.4 per cent) must be 
ascribed to the diffusion outwards of particles. 
As a result we cannot even in this case reckon 
with a component that is constant from day 
to day. We must also remember, that neither 
the regular 27-day fluctuations nor the re- 
current increases of intensity reveal themselves 
through a daily variation. Instead they con- 
stitute a trend which tends to diminish and 
displace the daily maximum. 

When q, has a positive value we are right 
to suppose that the inverted particle flow 
outside the corpuscular beams has been suffı- 
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ciently big to manifest itself. However, it is 
equally probable that these positive values 
are caused by statistical fluctuations combined 
with the errors in the correction for the 
deviation by the earth’s magnetic field. 


Discussion 


By comparing Figs. 4, 5, 6, and 7 we find 
that the 18h-component can be regarded as 
constant except, perhaps, during prominent 
magnetic storms. There are actually indica- 
tions that it is constant in the latter case too. 
As a function of time the amplitude is constant 
(Figs. 8 and 9). In one case only is there some 
contradicting evidence, namely the Nagoya 
measurements with a small aperture. However, 
when the limits of errors are taken into 
consideration, the points do not fall very far 
from the average. The difference in ampli- 
tude is not essential as regards the discussion 
of the constancy or non-constancy of the 
184-component. 

A certain variation of the 18h-component 
correlated to magnetic storms is not at variance 
with theory. According to BRUNBERG and 
DATTNER (1954) a change of the average 
energy of the cosmic ray particles is to be 
expected when the earth is in the vicinity of a 
corpuscular beam. In such a case the angular 
correction ought to be slightly different from 
that derived for an average energy of 20 GeV/c. 
Though the amplitude will then appear to 
change, it really remains constant. 

Another discrepancy is to be found in Fig. 
6, where the displacement of the 18"-compo- 
nent as registered with a large aperture on 
magnetically disturbed days is opposite to 
that found in Figs. 4 and 5. 

However, the disturbance of the 18®-com- 
ponent in connection with major magnetic 
storms is doubtful also for another reason. It 
has been shown (SIMPSON 1954, SANDSTROM 
1955) that the cosmic ray changes do not 
exactly coincide with the increase of geo- 
magnetic activity. Actually there is a phase 
shift of one or two days. As a result there 
should not be a very good correlation in this 
case, when the 18h-component is plotted as a 
function of geomagnetic activity. 

Summing up pros and cons it seems doubtful 
if the observed disturbances of the 18h-com- 
ponent during increased geomagnetic activity 
are real. 


Turning to the 12h-component, Figs. 4, 5, 6, 
and 7 display consistent features. This compo- 
nent undoubtedly varies in a systematic way 
with increasing geomagnetic activity. The 
variations are in full accordance with the theory 
of Alfvén. As a function of time the 12h- 
component does not vary very much, prob- 


ably owing to the fact that we must of 


necessity calculate with averages over com- 
paratively long periods. It appears from Fig. 9 
that the amplitude of the 12%-component 
was small during 1941—42. It cannot be 
mere chance that this period coincides with a 
minimum of solar activity. For quiet days it 
even appears to have changed into a midnight 
component. 

Again, as regards the 12h-components the 
magnitude of the variations shifts from one 
set of measurements to another. Part of 
the differences can be explained by the fact 
that the geomagnetic activity has been charac- 
terized in three different ways, through the 
K, and C indices and the Japanese Character 
Figure. However, this does not explain the 
difference between the measurements of 
Sekido and Kodama with small and big 
apertures. When we consider the comparison 
made in Fig. 10 and the accompanying re- 
marks it seems quite probable that some 
unknown factors cause rather large errors in 
the measurements with the small aperture 
telescope. Such an error might be due to the 
correction for atmospheric effects. The correc- 
tion for the barometric effect has been derived 
from measurements with telescopes having 
comparatively large apertures. It is by no 
means certain that the same correction should 
be applied to instruments with a small aperture. 
In a direction with the zenith angle of 30° 
the cosmic rays have a 1$ per cent longer 
path to traverse than in the zenith direction. 
The corrections for atmospheric effects should 
if possible be determined for each instrument 
individually. It ought also to be remembered 
that in the case of measurements with the 
small aperture telescope SEkIDO and KopAMA 
(1952) applied a correction for the temperature 
effect, which was omitted in the large aperture 
measurements. The temperature correction is 
no doubt less well defined than is the baro- 
metric correction. 

It is evident that the results presented here 
are purely qualitative. Most of the measure- 


Tellus VIII (1956), 1 


DIRECTIONAL COMPONENTS OF COSMIC RAYS 25 


ments indicate that the amplitude of the 18h- 
component is 0.13 per cent (average from 
the dotted lines in Figs. 4, 6, and 7) with the 
single value 0.38 per cent from the small 
aperture measurements of SEkIDO and Ko- 
DAMA (1952). 

More precise measurements are essential 
for quantitative studies of the directional 
components. The atmospheric effect is one 
of the factors limiting precision. Another 
such factor is the uncertainty about the 
average momentum of the primaries. To some 
extent one of these two factors can be over- 
come by measuring in two directions. If 
corrections are applied for the atmospheric 
effects in the conventional way, the two direc- 
tions give the average momentum (BRUN- 
BERG and DATTNER 1954) thus making it 
possible to find the angular correction for the 
time of maximum. 


Conclusions 


1:0 There is nothing in the experimental 
material presented in this study which is at 
variance with the assumption of one 18h- 
component due to a particle wind bound to 
the rotation of the sun and one 12-component 
consisting of a radial flow of particles (ALFVÉN 
1954). 

2:0 The 18%-component is invariable in 
correlation to geomagnetic activity of the 
order of RK, est, 

3:0 From the present material it is imposs- 
ible to decide whether the 18®-component 


is being disturbed or not in the case of major 
magnetic storms. Disturbances are not ne- 
cessarily contrary to theory (BRUNBERG and 
DATTNER 1954). 

4:0 The doubtful case of major magnetic 
storms excepted the amplitude of the 18h- 
component as a function of time appears to 
be constant. Observed variations are irregular 
and of the same order of magnitude as the 
limits of errors. 

5:0 The 124-component is correlated to 
geomagnetic activity. It is nonexistent or of 
reversed sign when geomagnetic activity is 
low. In accordance with theory the amplitude 
increases with increasing geomagnetic act- 
ivity. 

6:0 The 12-component as a function of 
time displays variations outside the limits of 
errors. Measurements from Huancayo indi- 
cate a correlation with the sunspot cycle. 

7:0 A seasonal variation reported by SE- 
KIDO and Kopama (1952) has been compared 
with partly contemporary measurements in 
Stockholm. From this comparison and with 
the accuracy given by the authors it appears 
that the reported seasonal variation is due to 


some local effect. 
* 
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On the Sub-tropical Jet Stream and its Role in the Development 


of Large-scale Convection 


By C. RAMASWAMY, International Meteorological Institute in Stockholm? 


(Manuscript received on October 18, 1955) 


Abstract 


This paper contains mainly the results of a synoptic and climatological study of the large- 
scale convection in northern India and Pakistan during the three months preceding the onset 
of the southwest monsoon. It has been shown that the sea-level and lower tropospheric charts 
give little clue to the development of the large-scale convection and that the latter is overwhelm- 
ingly determined by the divergence in the waves in the sub-tropical jet-stream. It has further 
been shown from detailed synoptic evidence that nor’westers, andhis and the majority of the 
thunderstorms without squalls in northern India and Pakistan in the pre-monsoon period are 
fundamentally the same phenomenon. The role of cold-air advection in the middle and upper 
troposphere in the development of large-scale convection has also been discussed. This study 
has further revealed that the regions of upper-divergence and convergence can be qualitatively 
located by identifying certain typical patterns on the high-level maps more than 12 hours 
before the usual time of commencement of convection and that, consequently, these maps can 
be used as effective tools in the issue of area-warnings against thunder in general and 
nor westers and andhis in particular. 

A general study has also been made of the large-scale convection in southeast Australia, 
Union of South Africa, Bechuanaland, Southern Rhodesia, northeast Argentina, Uruguay, 
southeast Brazil and southeast United States and the similarities between the large-scale con- 
vection in these countries and in Indo-Pakistan have been brought out. On the basis of these 
studies, it has been suggested that the jet stream plays the very important role of producing 
large-scale convection in the subtropics all over the world wherever it over-runs on its equator- 


ward side, moist air possessing a high degree of latent instability. 


1. The nor’wester and andhi problems of 
Northern India and Pakistan. 


From the middle of March to the time 
of establishment of the southwest monsoon, 
northern India (defined broadly, for purpose 
of this paper, as India north of latitude 20° N) 
and Pakistan experience thundersqualls or 
convective duststorms which are occasionally 
of great severity and cause heavy destruction 
of life and property. The thundersqualls which 
occur in Northeast India and East Pakistan, i.e. 
roughly between 83° E and 95° E and between 

1 The author belongs to the Indian Meteorological 


Service. He took part in the meteorological research 
activities at Stockholm when he was on ex-India leave. 


20° N and the foot of the Eastern Himalayas 
are popularly known as nor’westers. They are 
often accompanied by very heavy rain and 
sometimes by hailstorms and they occur gene- 
rally late in the afternoon or in the evening 
in the plains and late in the night or early in 
the morning in the valleys of Assam and along 
the foot of the Eastern Himalayas. Nor’- 
westers are, however, not uncommon over 
these regions in other parts of the day or night. 
The convective duststorms over the rest of 
northern India and over West Pakistan are 
popularly known as andhis. They are charac- 
terized by dustwalls and are invariably accom- 
panied by cumulonimbus clouds. Thunder is 
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also frequently heard during such storms, 
but the rainfall associated with them is usually 
small or negligible. They occur generally late 
in the afternoon or in the evening in the 
plains, but are not uncommon in those regions 
at other times of the day or night. Tornado 
clouds with the typical funnel shape have 
also sometimes been observed with the nor’- 
wester as well as the andhi. Squalls of roo mph 
have been recorded in association with these 
phenomena. 

Calcutta and Delhi, two of the great cities 
in India, experience the typical nor’ wester and 
the andhi respectively. The following sta- 
tistics of squalls of 40 mph or more at Calcutta 
and Delhi (K. P. Ramarrısunan and B. 
GOPINATH Rao 1954) is of interest. 


Number of squalls per 


month Usual date 
Station a Er GEORGE 
eee 52 southwest 
s | = =| = | &| monsoon 
À < 
eet} 
Calcutta 
(Alipore 2.6] 3.6] 5.2] 2.4] 0.2] —| First week 
22° 30° N of June 
88° 30’ E) 
Delhi 
(New Delhi | 0.6] 2.8] 7.0] 4.2} 2.0] 0.4] End of 
20235. N June 
72. 12 EB) 


N.B. When several distinct squalls occurred 
on the same day, a maximum of two have been 
counted with a view to avoid a large contribu- 
tion to the average by exceptional days. 


For purposes of comparison, the periods 
during June when the southwest monsoon 
usually sets in at Calcutta and New Delhi 
have been given by the present writer in the 
above table. These are based on the curves 
showing the normal dates of onset of the 
southwest monsoon published by the Govern- 
ment of India Meteorological Department in 
the Climatological Atlas for Airmen (1943 edi- 
tion, p. 3). As the date of establishment of the 
monsoon in any individual year can differ 
considerably from the normal date as given 
in the Atlas, it will be obvious from the above 
table that large-scale convection in northern 
India and Pakistan, as manifested in the nor’- 
wester and the andhi, markedly decreases with 
the establishment of the southwest monsoon. 
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2. Normal climatological features over North- 
ern India and Pakistan in the pre-monsoon 
period! 


The normal sea-level pattern over this 
region in May which is a typical month of 
large-scale convection is shown in Figure 1. 
It has been reproduced from the Climatological 
Atlas for Airmen (1943 edition, page 7) re- 
ferred to above” The isobars in this diagram 
refer to the mean pressure during the day. 
On this normal map, the present writer has 
superposed isopleths of normal dewpoints 
at 03 G.M.T. (i.e. 0830 hrs Indian Standard 
Time?) based on the data of about 200 stations. 
The trough lines based on the published 
isobars and surface winds have been drawn 
by the present writer. An examination of 
this sea-level map shows the following: 


(a) In association with the low-pressure area 
over West Pakistan and northern India, 
there is a normal incursion of moist air 
from the Bay of Bengal into Northeast 
India and East Pakistan (1e. east of 83° E 
north of 20° N) the moisture, as revealed 
by the dewpoints being a maximum in 
the coastal districts and gradually decreasing 
northwards in the interior. 

(b) In association with the same _ sea-level 
pressure-system, relatively very much 
drier air sweeps over the rest of northern 
India and over West Pakistan. 

(c) There are, broadly speaking, three air- 
streams over northern India and Pakistan 
at sea-level namely (1) westerlies or 
northwesterlies (2) southwesterlies or 
southerlies (3) easterlies. The southwester- 
lies or southerlies and the easterlies are 
both moist while the westerlies and north- 
westerlies are relatively very much drier. 
These three air-streams show a certain 
amount of velocity—convergence along 
the trough lines marked UC, CA and CM. 


1 For the Indo-Pakistan sub-continent as a whole, the 
following may be taken as the four seasons of the year: 

The cold weather period—December to February. 

The pre-monsoon or hot weather period—March to 
May. 

The monsoon period—June to September. 

The post-monsoon period—October and November. 

Unless otherwise specifically stated, the word “mon- 
soon’’ as used in this paper refers to the southwest mon- 
soon. 

2 The Indian Standard Time (I.S.T.) is 51/, hours 
ahead of Greenwich Mean Time. 
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Fig. 1. This map has been reproduced from the Climatological Atlas for Airmen (Government 
of India Publication, 1943 edition) and has been plotted according to the conventions in 
vogue at the time of its publication. On this map, the present author has drawn the trough 
lines and superposed isopleths (dotted lines) of the normal surface dewpoints at 03 G.M.T. 


(d) The normal west to east pressure gradient 
across West Bengal and East Pakistan is 
only about 4 millibars. 


The normal o9 G.M.T. stream lines and 
trough lines at 3,000 and 5,000 feet in May, 
which have not been reproduced here, show 
the following: 


(a) The easterlies and the trough-line CA 
seen on the normal sea-level chart, are 
not seen at 3,000 and 5,000 feet. 


The line UCM seen at sea-level and which 
more or less separates the very dry from 
the moist stream is persisting at the 3,000 
and 5,000 feet levels, but it shifts east- 
wards with increase in height. 


(b) 


(c) The normal winds are only $ to 13 knots 
even at 5,000 feet. 


Above 5,000 feet, the dry westerlies domi- 
nate the scene and steadily increase in speed 
with height up to the tropopause, the increase 
becoming marked above 10,000 feet. The 
normal contour-pattern at the 300 mb level 
in May is shown in Fig. 2. Almost all the 


winds on this map are 09 G.M.T. (1430 I.S.T.) 
wind normals! of resultant wind-direction and 
resultant velocity (i.e. taking into account 
directions) computed by the India Meteorolo- 
gical Department and based on all observations 
available upto 1950 in the case of Indian 
stations and upto 1947 in the case of Pakistan 
stations. The normals of the height of the 300 
mb isobaric surface and of the partial thick- 
ness of the s00~300 mb layer were computed 
by the present author from the daily radio- 
sonde observations at 15 G.M.T. for the years 
1951—ss published in the Indian Daily 
Weather Reports. The normals for Shillong 
are based on the data for the years 1953—55 
only. 


The 300 mb contour map clearly shows 
that there is a maximum of speed of the 
westerlies near about the latitude of New 
Delhi. This observational fact is consistent 
with the findings of Korzswaram and Par- 
THASARATHY (1954) about the mean position 


1 Normals based on less than 10 individual observations 
have not been included in the map. 
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Fig. 2. The figures against some of the stations are the normal heights of the 300 mb 
surface (upper figures) and the normal partial thicknesses of the 500 ~ 300 layer 
(lower figures) at 15 G.M.T. correct to the nearest tens of geopotential feet. The figure in 
the ten thousandth place has been omitted in the upper as well as the lower figures. For 
example, 179 and 260 against Calcutta indicate that the height of the 300 mb sur- 
face is 31,790 geopotential feet and that the partial thickness of the s00~300 mb 
layer is 12,600 geopotential feet. The contours have been drawn for every hundred 


geopotential feet. 


of the axis of the jet stream over India during 
the period April—May. 

In the pre-monsoon period, the highly 
moist air from the Bay of Bengal which is 
continuously invading Northeast India and 
East Pakistan possesses a high degree of latent 
instability (NORMAND, 1938. S. PETTERSSEN, 
1940). This instability cannot however be rea- 
lized on account of an inversion which 
separates the moist layer from the superin- 
cumbent dry westerlies. On the other hand, 
over Northwest India and West Pakistan, 
the lower tropospheric air is normally very 
dry and the lapse-rate in this air is very nearly 
dry adiabatic. The studies made by earlier 
workers on convective duststorms show that, 
prior to the development of these storms, 
there is an injection of moisture in the lower 
troposphere mainly as a result of movement 
of extra-tropical type of cyclonic circulations 
known in Indian meteorological literature as 
western disturbances (RIEHL 1954, p. 269). 
This “semi-moist” air has a high degree of 
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latent instability, the energy from which can 
be released only by a trigger of sufficient in- 
tensity. 


3. Present ideas about the mechanism of the 
nor’wester and the andhi 


The following is a very brief summary of 
the present ideas about the mechanism of the 
nor wester and the andhi: 

The nor’ wester—Almost the entire work on 
this subject has been confined to the study of 
the nor’westers over West Bengal! and East 
Pakistan.t In 1938, K. R. Ramanathan and 
K. P. Ramakrishnan put forward the view that 
the development of the nor’wester is primarily 
determined by cold-air advection above 6 km. 
The authors based their conclusion on the 


1 The region referred to as West Bengal and East 
Pakistan in this paper was politically known as ““Bengal”’ 
before 15 August 1947. Hence the expression “‘nor’ westers 
of Bengal” has gained currency in scientific literature. 
Destructive nor’westers, however, occur in the State of 
Assam and also occasionally in other States in and near 
Northeast India. 
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monthly mean upper-air temperatures com- 
puted from the change of the monthly mean 
winds with height. Surprisingly enough, their 
view about the mechanism receded into the 
back-ground and in spite of an attempt by S. L. 
MALURKHAR (1948, 1949) to bring it again tothe 
forefront, the nor’ wester is considered today as 
an essentially lower tropospheric phenomenon. 
The workers have been seeking for the cause 
of this phenomenon in the accentuation of the 
west to east pressure-gradient across“ Bengal”, 
in insolation, the “convergence-front” between 
moist air and dry air in the lower troposphere, 
cold fronts of the middle-latitude type associa- 
ted with western disturbances, cold air from a 
neighbouring thunderstorm, katabatic winds 
from the hills, and so on (A. K. Roy 1949, 
1950, B. N. Desai 1950). 

The andhi—This phenomenon has been 
studied at only two or three individual sta- 
tions (P. R. KRISHNA Rao 1938, SREENIVA- 
SAIAH and SUR 1939, DEsAI and Mat 1940, 
K. L. SINHA 1952). According to these wor- 
kers, the andhis are associated with the re- 
lease of the latent instability energy in the 
moist air brought by the western disturbances, 
the trigger for the release being insolation or 
cold-fronts of western disturbances or lower 
tropospheric convergence near a discon- 
tinuity. 

At this stage, it would not perhaps be irre- 
levant to mention that the present writer and 
B. L. Bose, on the basis of a study of the mean 
isotherms for the 700 ~ 500 mb and 500 ~ 300 
mb layers for daily situations, had tentatively 
suggested (C. Ramaswamy and B. L. Bose 
1954) that the development of the nor’ westers 
was determined by (1) the dynamics of the 
middle and upper tropospheric west-wind 
belt (2) the thermal advective processes as 
observed from tne change of wind with height, 
in the 700~ soo mb and soo” 300 mb 
layers. The authors had also made a general 
suggestion that similar processes were prob- 
ably responsible for the convective dust- 
storms of Northwest India and West Pakistan. 


4. Technique of analysis in the present 
investigation 
The investigation reported in this paper is 
an examination of the entire probiem de novo 
without any reference to the previous work 
of the author and is based on charts prepared 


and analysed by the ‘author in the International 
Meteorological Institute, Stockholm, Sweden. 
The technique adopted consisted, broadly 
speaking, of analyses of (1) sea-level and upper- 
air maps for 80 daily situations (2) vertical 
time-sections of representative stations (3) 
mean upper-air maps for one unusually long 
spell of large-scale convection and (4) normal 
sea-level and upper-air maps for March, April 
and May in relation to the normal number of 
days of thunder and normal number of squalls 
per month. 

The daily situations for the following 80 
days were studied. The dates? given below 
refer to the soo and 300 mb charts. The corre- 
sponding sea-level and lower tropospheric 
upper-air maps refer to the next morning. 


Daily situations 

12 to 15 April 1952 and 17 April 1952. 

28 and 29 April 1952. 

4 and 5 June 1952. 

6 to 12 March 1953. 

17, 18 and 19 April 1953. 

26 April 1953 (vertical time-sections were 
prepared for the period 24 April to 2 May 
1953 besides those for other periods). 

7. I to 31 May 1953. 

8. 27 March to 14 April 1954. 

9. 27 April to 6 May 1954. 


AWRY NH 


In addition to the above 80 days, the sea- 
level and soo and 300 mb maps which corre- 
sponded to the record squall of ror mph at 
Allahabad (25° 27’, N 81° 44’ E) on 21 March 
1950 were also examined. 

The above days were selected mainly with 
reference to the weather developments east 
of 83° E: they included many severe and 
destructive squalls, a number of which consti- 
tuted records for the stations concerned. They 
also included a large number of cases of severe 
convection between 70° E and 83° E. 

For the majority of these situations, the 
sea-level maps for 03 G.M.T. (i.e. for 0830 
hours Indian Standard Time) published in 
the Indian Daily Weather Reports, which are 
based on a very large amount of data were 
utilized. The upper-air maps were prepared 
mainly from the data published in the Indian 


! Unless otherwise explicitly stated, all dates in this 
paper refer to the 500 and 300 mb charts and not to the 
sea-level or lower tropospheric charts or the weather 
developments. 
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Daily Weather Reports. These reports were 
very kindly placed at the disposal of the 
author by the Director, Swedish Meteorolog- 
ical and Hydrological Institute, Stockholm, 
Sweden. 


On the sea-level maps, trough lines including 
those corresponding to CA, CM, CU in Fig. 
I were marked. Trough lines and stream 
lines were also drawn on the lower tropospheric 
maps in the simple manner suggested by 
RIEHL (1954, p. 189). On all the sea-level 
(03 G.M.T.) and lower tropospheric maps 
(02 G.M.T.) were superposed the thunder- 
storms (including “thunder heard”, ww=17), 
convective duststorms and squalls which devel- 
oped between 03 G.M.T. of that day and 03 
G.M.T. of the following day. These data for 
northern India and Pakistan were specially 
obtained from the India Meteorological De- 
partment so that they may be as complete 
and reliable as possible. Similar data for 
southern India and Burma (the discussion of 
the weather in these regions is outside the 
scope of the present paper) were taken from 
the Indian Daily Weather Reports. The Pa- 
kistan data are not complete. They may not 
also be free from accidental errors, such as 
mutilations, as they were received in India 
through the usual telecommunication channels. 


The soo and 300 mb! maps contained Pilot 
Balloon observations at 09 G.M.T. and Radio- 
wind and Radiosonde observations at ıs G.M. 
T. This difference in time between the two 
sets of observations was unavoidable as there 
was no other way of getting a reasonably 
complete picture of the flow patterns in the 
upper atmosphere. The mean hour of these 
composite maps corresponded TOME AGAIN 
(ic. 1730 Indian Standard Time). Wind-data 
for 18,000 feet were obtained from the India 
Meteorological Department and utilized for 
the soo mb level for stations at which the 
balloons did not reach 20,000 feet. Similarly 
wind-data for 25,000 feet were utilized for 
the 300 mb level when the balloons did not 
reach 30,000 feet. Data in the Indian Daily 
Weather Reports which appeared doubt- 
ful due to printing or other mistakes were 
referred to the India Meteorological Depart- 


1 250 mb charts were also prepared for days on which 
the data for this level were adequate enough to be a 
useful supplement to the 300 mb charts. 
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ment and corrected on the basis of informa- 
tion supplied by that department. 

On these maps, the partial thermal winds 
and partial thickness values for the 700 ~ 
500 mb and $00 ~ 300 mb layers were plotted. 
The thunderstorms (including “thunder heard” 
ww = 17), convective duststorms and squalls 
which developed between 03 G.M.T. of the 
next day and 03 G.M.T. of the following day 
(i.e. the same weather developments as plotted 
on the sea-level and lower tropospheric maps 
for the next morning) were superposed on 
these maps). The main reason for adopting 
this procedure was that the o9 G.M.T. Pilot 
Balloon data on this chart cannot normally be 
available for analysis before 13 G.M.T. (1830 
1.S.T.) by which time the afternoon convection 
in the plains would have started. And what is 
more important, the Rawin and Radiosonde 
data of 15 G.M.T. cannot be available for 
analysis before 18 G.M.T. (2330 I.S.T.) by 
which time, the convection would be over 
or be in the dying phase. Consequently, the 
composite chart of 09/15 G.M.T. which is 
the best one available at present for high level 
analysis can be used for forecasting convection 
in the plains only on the following day. 

As regards the choice of the period 03 G.M. 
T. to 03 G.M.T. for the plotting of convec- 
tion-data, it may be mentioned that the hour 
03 G.M.T. is in many respects, the most im- 
portant synoptic hour in the Indian region 
and, as such, is a convenient hour of reference 
in discussions on forecasting. It is also the 
nearest principal synoptic hour before con- 
vection begins to develop in the plains. 

In the analysis of the soo and 300 mb con- 
tours, considerable attention was given to the 
principle of continuity. Much greater weight 
was also given to the winds than to radiosonde 
data in view of (1) the reasons pointed out by 
Riehl (1954, p. 183) about tropical radiosonde 
data and (2) of the general consensus of opinion 
among experienced analysts of Indian data 
(A. K. Roy and N. C. RAISARCAR, 1955) that, 
in the Indian region, at present, pilot balloon 
data are “very much more dependable” than 
radiosonde data. 

It was often found difficult to make the 
spacing of the contours on the basis of radio- 
sonde data at 15 G.M.T. consistent with the 
pilot balloon winds at 09 G.M.T. or the radio 
winds at 15 G.M.T. This difficulty which was 
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particularly felt at the 300 mb level is per- 
haps partly due to errors in the radiosonde 
data. The exact reasons for these “inconsisten- 
cies” are, however, not clear and should be 
thoroughly investigated as a separate problem 
(RIEHL, 1953, p. 193). In view of these unknown 
factors, it was felt that no special advantage 
would be gained by spacing the lines with a 
geostrophic wind-scale. This slightly increased 
the subjectivity in the analysis which was 
unfortunate but unavoidable. However, as 
the purpose of this paper was merely to 
establish purely qualitative relationships be- 
tween high level patterns and large-scale con- 
vection, this additional subjectivity in anal- 
ysis did not in any way affect the general 
conclusions summarized in the later paragraphs. 


5. Sea-level and lower tropospheric circula- 
tion patterns in relation to large-scale 
convection 


The following is a brief summary of the 
sea-level and lower tropospheric analysis. 
The dates given in brackets against some of 
the paragraphs, refer to the dates of the sea- 
level and lower tropospheric charts. The instances 
cited are not the only ones to illustrate each 
type of situation: many others could also be 
given. 


1. East of 83° E. north of 20° N. 
Sea-level pressure and wind systems. 


(a) The sea-level west to east pressure- 
gradient had no significant association with 
the development of large-scale convection. 
Very weak east to west pressure-gradients 
were sometimes found to be associated with 
widespread thunderstorms and severe squalls 
(14, 16 and 18 April 1952 and 19 and 20 April 
1953) in parts of Northeast India and East 
Pakistan while pressure-gradients in excess of 
the normal for the month were a number of 
times associated with only isolated thunder- 
storms or fair weather west of 90° E south of 
25° N (9 March 1953, 6 April 1954). 

(b) Sometimes, the seasonal trough of low 
pressure appeared as a closed low (isobars 
drawn at 2 mb intervals). No significance 
could be attached to the appearance of this 
closed low or its position or intensity so far as 
the subsequent convective developments were 
concerned. Extensive squalls, a few of which 


caused very heavy destruction, developed over 
Northeast India and the adjoining districts o 
East Pakistan without any closed low at sea 
level (27 April 1953) while fair weather pre- 
vailed over West Bengal and the adjoining 
parts of East Pakistan in spite of closed lows 
over or near these areas (8 March 1953, 9, 
March 1953, 7 May 1953, 5 April 1954, 13 
April 1954). Squalls, widespread thunderstorms 
and moderate convective rain also sometimes 
developed over the State of Assam and the 
contiguous districts of East Pakistan with a 
ridge extending from the Tibetan or upper 
Burman border into East Pakistan through 
Assam (14 April 1954, 5 May 1954). 

(c) The trough lines CU and CM had no 
association with convective developments 
beyond what may be expected by mere 
chance. The trough line CA was, however, 
slightly better associated with the develop- 
ments. The association was, however, far too 
poor to consider it as the determining factor. 

No front of the middle-latitude type or 
tropical cyclone or depression moved across 
the region under discussion. 


Lower tropospheric flow-patterns 


(a) The streamlines at the 3,000 and 5,000 
feet levels showed that the convective weather 
developed mostly in the southwesterlies or 
southerlies or in the easterlies, i.e. in the 
moisture-laden air. 

(b) The trough-lines CU and CM did not 
have any association with the thunderstorm 
field. The number of thunderstorms and squalls 
which developed along and near the CA 
trough line was only a small percentage of 
the total number of thunderstorms, and squalls. 
Severe convective developments took place 
even where the streamlines were anticycloni- 
cally curving at 5,000 feet. A good instance 
of this type was on the 19 April 1953 near 
Gaya and Allahabad. These stations, in spite 
of the anticyclonic curl of the stream-lines 
near their vicinity, experienced heavy squalls 
on the same afternoon. Gaya (24° 45’ N, 
84° $7’ E) experienced a record squall of 100 
mph and Allahabad (25° 27’ N, 81° 44’ E) a 
squall of 60 mph. 

(c) There were one or two occasions when 
feeble western disturbances seemed to have 
moved across Northeast India and East 
Pakistan in the lower troposphere. Also, on 
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some days, cyclonic circulations appeared in 
certain parts of Northeast India in the lower 
troposphere. But they died in situ within 24 
or 48 hours. Neither the western disturbances 
as they appeared at the 5,000 feet level nor 
the locally-developed cyclonic circulations had 
any direct connection with the convective devel- 
opments. 

Figures 3, 4 and 5 show the sea-level isobars 
and trough lines CA, CM, CU (compare 
Fig. 1), the 24 hours isallobars and the 3,000 
and 5,000 feet streamlines at 02/03 G.M.T. on 
14 May 1953. The upper winds at 06 G.M.T. 
have been plotted with broken arrows. The 
hatched areas in Fig. 3 show the region where 
the thunderstorms developed and the cross- 
hatched areas' show the region where thunder- 
squalls or convective duststorms developed 
between 03 G.M.T. of 14 and 03 G.M.T. of 
15 May 1953. According to Press Reports, 
the squalls caused damage to property, inter- 
ruption in communications and loss of several 
human lives in Gangetic West Bengal. Cal- 
cutta itself experienced two squalls, well- 
separated from cach other,—a rather unusual 
development. 


Figures 6, 7 and 8 show the corresponding 
charts on 28 May 1953 and the convective 
developments during the subsequent 24 hours. 
It will be seen that the convective develop- 
ments took place only in Assam and the 
contiguous districts of East Pakistan and that 
fair weather prevailed practically over the 
entire region to the west of 90° E. 


These two sets of charts illustrate some of 
the conclusions summarized above regarding 
regions east of 83° E. The charts look so alike 
and yet were associated with such different 
weather developments. Indeed, in some respects 
e.g., the positions of the trough lines, the 
depth of the moist air, the speed of the winds 
in the moist stream and the 24 hours isallobars, 
the charts for 28 May might, according to 
the conventional ideas, be considered as more 
favourable than those for 14 May 1953 for 
the development of thundersqualls in and 
around Gangetic West Bengal (i.e. between 
21° N and 25° N and 85° E and 89° £). It 
may also be noted that the east to west sea- 

1 The hatched and cross-hatched areas show only 
the general distribution of the storms and squalls. They 


do not imply that each and every observatory in those 
areas reported the phenomena concerned. 
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level pressure-gradients over West Bengal 
and East Pakistan on these two days did not 
differ significantly from the normal gradient 
for the month as seen in Fig-1. 


Il. West of 83° E, north of 20° N 

Sea-level pressure and wind-systems. 

The extensive seasonal low-pressure area to 
the west of 83° E vide Fig. 1, had only one 
closed low at sea-level (isobars drawn at 2 mb 
intervals) on about 43% of the daily situations, 
the closed isobar lying mainly to the west 
of 74° E. On about 9% of the days, the low 
pressure area had two closed lows west of 
83° E and on 6% of the days only, any of 
the lows had more than one closed isobar. 
In the preparation of this statistics, the closed 
lows in Northeast India, which sometimes 
extended to the west of 83° E have not been 
taken into account. The formation, move- 
ment and disappearance of the lows with 
only one closed isobar were often erratic 
(15 to 17 May 1953, 21 to 23 May 1953 and 
2 to 7 May 1954). And neither the lows 
with one closed isobar nor even those with 
more than one nor the trough-lines associated 
with these sea-level systems had any signi- 
ficant connection with large-scale convection 
(18 and 30 April 1952, 18 and 19 April 1953, 
17, 18 and 21 May 1953). 


Lower tropospheric flow-patterns 
(a) The thunderstorms, convective dust- 


storms and squalls developed, broadly speaking, 
in the following three sectors of the wind- 
field at the 5,000 feet level: SSW/WSW, 
E/SE and WNW/NNW. For facility of 
easy reference, these will be referred to in the 
following discussion as the southwesterlies, 
the eastsoutheasterlies, and the northwesterlies. 
The southwesterlies and the eastsoutheasterlies 
are, as is well-known, drawn into the region 
by the western disturbances. Of the total 
number of thunderstorms, convective dust- 
storms and squalls which developed over the 
plain stations west of 83° E north of 20°N, 
39 % developed in the southwesterlies, 16 % 
in the eastsoutheasterlies and 45% in the 
northwesterlies at the sooo feet level. This 
statistics is of course rough and no undue 
importance should be attached to the actual 
percentages. Nevertheless, the relative values 
of the figures do reveal that large-scale convec- 
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tion develops in the northwesterlies much 
more frequently than it is generally believed: 
it probably develops in this stream as fre- 
quently as in the southwesterlies. The im- 
portance of this fact will become obvious in 
the light of the discussion in sub-paragraphs 
(c) and (d) below. 

(b) The large-scale convection developed 
in the three streams mentioned above, most 
often within the body of each of these streams 
and quite far away from the trough-lines 
associated with these systems and very often 
in anticyclonically curving stream-lines (15 and 
30 April 1952 and 6 June 1952, 18 and 19 
April 1953, 15, 16 and 17 May 1953). 

(c) Of the storms and squalls which deve- 
loped in the northwesterlies, the majority deve- 
loped in advance of the ridge associated with 
the incursion of the “semi-moist” south- 
westerlies from the north Arabian Sea. In 
these northwesterlies, there was obviously 
no agency for convergence. The rest of the 
storms and squalls in this stream developed 
in a more or less straight current in which it 
was not possible to find any significant con- 
vergence. Among the cases studied in this 
paper, the most spectacular case of large- 
scale convection in northwesterlies, occurred 
on 30 April 1952 when extensive andhis and/or 
thunderstorms developed over an area of 


Fig. 3. The thick continuous lines are trough lines UC, 
CA, CM and the thin dotted lines are isopleths of 24 hrs 
pressure-changes at 03 G.M.T. drawn at 2 mb intervals. 
The two figures to the left rear o1 each station circle are 
dew-points at 03 G.M.T. Compare with Figure 6 especi- 
ally with the areas of thunderstorms (hatched) and of 
thundersqualls or convective duststorms (cross-hatched). 


Fig. 4. The streamlines have been drawn for the o2 
G.M.T. winds. Compare with Fig. 7 and with the areas 
of convective developments in Fig. 3. 


Fig. 5. The streamlines have been drawn for the o2 
G.M.T. winds. Compare with Fig. 8 and with the areas 
of convective developments in Fig. 3. 


Fig. 6. Same convention in plotting as in Fig. 3. Compare 
with Fig. 3 especially with the areas of convective deve- 
lopments. 


Fig. 7. The streamlines have been drawn for the 02 
G. M.T. winds. Compare with Fig. 4 and with the areas 
of convective developments in Fig. 6. 


Fig. 8. The stream-lines have been drawn for the 02 
G.M.T. winds. Compare with Fig. 5 and with the areas 
of convective developments in Fig. 6. 


Tellus VIII (1956), 1 


about 88,000 sq. miles with Delhi near about 
the centre of this area. This case deserves 
special study not only because it is meteorolo- 
gically of great interest but also because, 
according to Press reports, a Dakota aircraft 
of a scheduled passenger service crashed a few 
miles from Delhi just before landing when 
that station was in the grip of the andhi and 
the aircraft and all the passengers and the 
aircrew were lost. Other good instances of 
thunderstorms, convective duststorms and 
squalls in northwesterlies were on 10, 16, 17, 
18 and 20 May 1953 and 6 May 1954. 

(d) In any discussion on large-scale convec- 
tion in the northwesterlies, it is necessary to 
remember that the very dry stream over this 
region on a large number of days in this 
season is also a northwesterly current and that 
it is characterized by clear or partly clouded 
skies (A. K. Roy, 1946). In view of this and 
also of the statistics and other conclusions 
presented in sub-paragraphs (a) and (c) above, 
it would appear that the forecaster who 
draws his conclusions only from the lower 
tropospheric charts and proceeds on the 
assumption that no deviation from the average 
(RIEHL 1954, p. 178) in the lower tropospheric 
flow-patterns necessarily indicates fair weather, 
runs the risk of making incorrect prognoses 
on a large number of occasions of large-scale 
convection west of 83° E. 

(ec) No cold-front of the middle-latitude 
type (RIEHL 1954, page 237) caused large- 
scale convection in the region under discussion. 

Considerations of space and the need for 
restricting the number of diagrams prevent 
us from reproducing in this paper, the sea- 
level or lower tropospheric charts for any of 
the days mentioned above. However, Figs. 3 
to 8 which had been selected for illustrating 
weather situations east of 83° E also inci- 
dentally illustrate some of the conclusions 
regarding regions west of 83° E. It will be 
seen from Figs. 3, 4 and 5 that the cross- 
hatched areas west of 83° E where there were 
convective duststorms and some thunder- 
storms (there were also showers between and 
around these areas) lay in a region of rising 
sea-level pressures and in a region in which 
there were no cyclonic wind-shifts. In parti- 
cular, it may be noted that the large-scale 
convection between 75° E and 80°E and between 
20° N and 25° N developed in a northwesterly 
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stream in advance of a sea-level and lower tropo- 
spheric ridge. In contrast to this, fair weather 
prevailed over the same region in association 
with northwesterlies and with falling sea-level 
pressures vide Figures 6 to 8. 


6. Analysis of 500 and 300 mb patterns on 
individual days 


Unfortunately, adequate wind-data were not 
available on quite a large number of days to 
fix the 300 mb patterns with reasonable 
certainty. Nevertheless sufficient number of 
situations was available in which the data 
on the soo and 300 mb levels taken together 
could enable the author to draw the following 
purely qualitative conclusions: 


(a) The convective developments (i.e. con- 
vective duststorms and thunderstorms with or 
without squalls) over every part of northern 
India and Pakistan had far better association 
with the soo and 300 mb vorticity patterns 
than with the sea-level maps and lower tro- 
pospheric flow patterns. 


(b) The convective developments took 
place mostly in advance of a trough or in 
the rear of a ridge at the 300 mb level. This 
conclusion was also true in a large number of 
cases, for the 500 mb level. However, as the 
wave-patterns sloped westwards with height, 
more thunderstorms were found in the rear 
of the trough line at the soo mb level than at 
the 300 mb level. In a number of cases, the 
association with the soo mb patterns was poor 
but was quite satisfactory at the 300 mb 
level. 

(c) Other conditions being the same, the 
greater the general intensity of the jet stream, 
the more extensive were the thunderstorms 
and the more severe in general were the squalls. 
It is interesting to note that the record squalls 
were associated with high-speed jets. 

(d) Quite stable weather prevailed on a 
large number of occasions in the rear of a 
trough and in advance of a ridge although 
there was a steady incursion of latently un- 
stable moist air in the lower troposphere over 
those regions. 

(e) Unusual fair weather also prevailed 
over the major part of Northeast India and 
East Pakistan when the winds over that region 
were exceptionally light in the upper tro- 
posphere. This situation arose when the axis 


~ 


of the jet stream was oriented from southwest 
to northeast and had shifted far to the west 
and northwest of the region of light winds. 

(f) As should be expected from the kine- 
matic and hydrostatic properties of waves in 
the westerlies (RossBy 1946), the partial thick- 
ness lines for the 700 ~ 500 mb and 500 ~ 300 
mb layers showed, in general, the wave 
patterns in the jet-stream more conspicuously 
than the soo mb and 300 mb contour-patterns 
respectively. In a large number of cases, the 
partial thermal winds for these layers also 
indicated that the actual winds were backing 
with height in the regions in which there was 
subsequent development of large-scale con- 
vection. 


The high-level contour patterns associated 
with the convective developments have been 
schematically classified into the following 5 
types and shown in Figures 9 (a) to 9 (e). 
For facility of easy reference, the various 
types have been designated by suitable abbre- 
viations shown against each type. 


engl Schematic Pattern is 
No. tion 
I Uniform shear and sinusoidal| Su-Sin 
variation in curvature. 
2 | Uniformshearandrapiddecrease | Su-T 
of curvature ahead of trough- 
line. 
3 | Uniform shear and rapid in-| Su-R 
crease of curvaturein the rear 
of ridge-line. 
4 | Varying shear and rapidly di-| Sv-T 
verging contours mainly ahead 
of trough-line. 
5 | Varying shear and rapidly con-| Sy-R 


verging contours mainly in 
rear of ridge-line 


In this connection the following points 
should be emphasized 


(i) The schematic patterns pre-suppose the 
existence of latently unstable moist air 
in the lower troposphere over the region 
covered by the respective patterns. 

(ii) The schematic patterns are purely. ten- 
tative. In particular, it may be mentioned 
that clear-cut cases of Su-T could not be 
found during this study as it was difficult 
with the usual positions of the troughs 
and the ridges and the available number 
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(iii) 


(v) 
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of observations at the soo and 300 mb 
levels to state categorically that they 
were different from those represented 
by Su-Sin. 

The schematic patterns by no means 
exhaust the possible types of patterns 
which can lead to large-scale convection. 


The patterns are strictly true only for 
the 300 mb level. However, in the absence 
of observations at the 300 mb level, they 
may be taken as very roughly true of the 
soo mb level also bearing in mind that 
the convective developments will gener- 
ally be nearer the trough-line and further 
away from the ridge-line at this level. 


Figures 10 and 11 show the soo and 300 


mb 


contour-patterns at 09/15 G.M.T. on 


18 April 1953 and the convective phenomena 
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LARGE-SCALE CONVECTION IN NORTHERN INDIA & PAKISTAN 
IN RELATON TO UPPER TROPOSPHERIC CONTOUR PATTERNS 
(SCHEMATIC DIAGRAMS) 
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Fig. 9. The relative positions of the thunderstorms, 
duststorms or squalls (symbol q) or their actual number 
as shown in the schematic patterns have no special signif- 
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Fig. 10. Arrows indicated by continuous lines are actual winds while arrows indica- 
ted by broken lines are the partial thermal winds. The letter R indicates that data 
are based on Rawin observations. The figure 18 by the side of some of the wind- 
arrows indicates that the winds are for 18.000 feet level. The radiosonde heights and 
the partial thicknesses have been plotted and the contours have been drawn accord- 
ing to the same convention as in Fig. 2. The symbol q indicates squalls. The 
convective phenomena at some of the stations have been plotted in slightly altered 
positions, for lack of space. The phenomena south of Lat. 20° N have also been 
shown for the sake of completeness but are outside the scope of discussion in the 
present paper. 
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300 MB CONTOURS 
09/15 GMT 18/4/1953 
rR,%,9,V a 
03 GMT 19/4 -03 GMT 20/4/1953 


Fig. 11. The figure 25 by the side of some of the wind-arrows indicates that the 
winds are for the 25,000 feet level. Other conventions in plotting are the same as 
in Fig. Io. 


19 April and 03 G.M.T. of 20 April 1953. 
As convection mostly develops along the foot 
of the Himalayas and in the valleys of Assam 
in the night or very early in the morning 
when observers maintain watch only at a 
very small number of stations, the thunder- 
storms etc. plotted on the maps in these regions 
do not give a complete picture of the actual 
convective developments. In order to make 
up this deficiency, the symbol for showers has 
been plotted at stations which did not report a 
thunderstorm or squall or convective duststorm 
but which reported rainfall during the past 
24 hours in their next morning’s routine 
weather message (i.e. at 0830 I.S.T. or 0300 
G.M.T. on 20 April 1953). Similar procedure 
was also followed in respect of the other 
cases studied in this paper. The patterns seen 
on the soo and 300 mb maps of 18 April 
1953 are of the Sy-Sin type. Among the nu- 
merous andhis and nor’westers seen in the 
diagrams, the reader’s attention is specially 
invited to the record nor’wester squall of 
100 mph (English miles) at Gaya (24° 45’ N, 
84° 57’ E) and a 60 mph andhi at Allahabad 


(25° 27’ N, 81° 44’ E). These stations lie in a 
region, which is normally not susceptible to 
large-scale convection, vide section 8. The 
sea-level and the 5,000 feet patterns for this 
case have already been discussed in section 5. 

Figures 12 and 13 show the soo and 300 mb 
patterns at 09/15 G.M.T. on 7 March 1953 
and the convective developments between 
03 G.M.T. of 8 March 1953 and 03 G.M.T. 
of 9 March 1953. The patterns over Northeast 
India and East Pakistan are of Sy-T type. The 
thunderstorms (and showers) mainly lie ahead 
of the trough-line at the 300 mb level but 
some lie behind the trough-line also in re- 
gions where the contours are diverging and 
the winds are comparatively stronger. The 
day illustrated is one during the spell 6 to 12 
March 1953 when there was extensive con- 
vection in and around the State of Assam. 
It may be of interest to observe that the jet 
stream strengthened after 7 March; on 12 
March, it was intense and was responsible for 
very severe squalls at some places in the area 
of thunderstorms in the above maps, east of 
85° E. Unfortunately, on account of the per- 
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Fig. 12. Same conventions in plotting and drawing of contours as in Fig. 10. 
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Fig. 13. Same conventions in plotting and drawing of contours as in Fig. 11. 
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300 MB CONTOURS 
09/15 GMT 5/5/1954 
R,®,9,V 
03 GMT 6/5 - 03 GMT 7/5/1954 


Fig. 14. Same conventions in plotting and drawing of contours as in Fig. 11. 


sistent low clouds, no pilot balloon obser- 
vations were available east of Calcutta on 
this day (radiosonde observations of Shillong 
were not available for this spell). Hence the 
vorticity-patterns in the region of convection 
on this remarkable day can only be inferred 
from the patterns west of 89° E and from 
their similarity with those for other days 
when observations were available east of 89° E 
during this spell and the similar Sy-T spell 
between 27 March 1954 and 14 April 1954. 

It may be incidentally mentioned that the 
region under discussion consists largely of 
valleys and hills and as such, there can be no 
doubt that orography plays some part in 
producing convection in this region. However, 
it has to be noted that similar weather did 
not develop on days on which the lower 
tropospheric flow patterns were far more 
favourable for orographic developments. Fur- 
ther, an examination of a relief map of India 
shows that the hills in and around Assam are 
so situated that any convective clouds which 
develop over them due to orography would 


not drift into the two important valleys in 
that State (with which we are mainly con- 
cerned) but would move along the hills on 
account of the middle and upper tropospheric 
flow patterns. The discussions and the diagrams 
in Section 7 of this paper are also quite relevant 
in this connection. In the circumstances, there 
cannot be any doubt that, in the large number 
of cases of convection in and around Assam 
studied in this paper, the contribution by 
orography was much less than that by the 
jet-stream waves. Incidentally, these remarks 
are equally true of large-scale convection along 
the foot of the western Himalayas. 

The convective phenomena between 75° 
and 80° E, north of 30° N in Figs 12 and 13 
have partly occurred in the hill-stations in the 
western Himalayas and partly in the submon- 
tane districts. While orography has obviously 
played a part in the development, the Sy-R 
patterns seen in that region (see particularly 
the winds on the soo mb chart in Fig. 12) 
have also contributed to the development, 
especially in the submontane districts. 
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Fig. 14 shows the 300 mb pattern at 09/15 
G. M.T. on 5 May 1954 and the convective 
developments between 03 G.M.T. of 6 May 
1954 and 03 G.M.T. of 7 May 1954. It is a 
typical instance to show how convective deve- 
lopments in very distant regions which may appear 
on the sea-level and lower tropospheric charts as 
quite disconnected and even random developments 
(see Section 5) may appear as well-connected and 
orderly developments on the vorticity pattern at 
the 300 mb level. It also illustrates the paradox 
that under suitable conditions severe weather 
can develop “at a ridge”. The convective 
developments in extreme northwest India are 
of type Su-R while those in and around 
Gangetic West Bengal are of type Sy-R. 
The weak SSW or SW winds at Visakha- 
patnam and Masulipatam which are in lower 
latitudes and the strong W winds at Calcutta 
which is in a higher latitude, indicate con- 
verging contours in the rear of a ridge. With 
regard to the northwesterly wind at Jamshed- 
pur (west of Calcutta) which does not fit in 
with the contours, it may be mentioned that 
the wind refers to 25,000 feet and that the 
winds over Calcutta also, on this day at 09 
G.M.T., were 320 degrees at 25,000 feet but 
backed with height and became 260 degrees 
at 30,000 feet (300 mb-level). 

Among the cases studied in this paper, those 
relating to 13 and 29 April 1952, 26 April 
1065 8, 134020, 24,) a7 and 31 May +1953, 
28 March 1954, 5 and 28 April 1954 and 
2 May 1954 are of great interest!, as the upper 
tropospheric patterns on these days led to 
“sudden” and/or severe convective develop- 
ments (andhis and/or nor’westers) or to the 
development of thundersqualls after a long 
break or to exceptional fair weather. Conside- 
rations of space do not permit us to go into 
these cases here. We would, however, like 
to describe below at least in a few sentences 
what was perhaps the most unique case in 
the above list. 

26 April 1953.—This led to the most exten- 
sive development of squalls in Northeast India 
studied during the present investigation and 
caused very heavy destruction at Barrackpore, 
a suburb of Calcutta. The 300 and 250 mb 
charts show a Su-Sin pattern. The velocity- 


1 The sea-level and the lower tropospheric situations 
for a number of these cases have been very briefly referred 
to in section 5. 
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profile on this day at the 250 mb level along 
the meridian 80° E shows a remarkably sharp 
“knife-edged” jet. The lower tropospheric situa- 
tion on the morning of 27 April 1953 has 
already been referred to in Section 5. The 
vertical time-section over Calcutta on this 
day may be seen in Fig. 15. 

Figures 15, 16 and 17 show the vertical 
time-sections over Calcutta during three 
typical spells. Fig. 15 represents a spell in 
which a trough lay with its axis to the west 
of Calcutta at the 250 mb level with varying 
intensities and slightly different positions. Fig. 
16 represents the situation during a part of 
the spell of vigorous convection in and around 
the state of Assam between 27 March 1954 
and 14 April 1954 (vide section 7) when the 
trough at the soo mb level and higher levels 
lay with its axis to the east of Calcutta. Fig. 17 
shows a spell characterized by weak vertical 
wind shear. In studying these three diagrams, 
the convective developments on each day have 
to be compared with the upper winds, mixing 
ratios and 24 hours pressure changes of the 
previous evening and also with the 24 hours 
pressure-changes at 03 G.M.T. on the day of 
the development. The winds, mixing ratios, 
etc. in each spell should also be compared with 
the corresponding elements in the two other 
spells with a view to assess their importance 
in the development of convection. It will be 
observed that. 

(a) The upper winds in the middle and 
upper troposphere had by far better association 
with the convective developments than the 
winds in the lower troposphere or the mixing 
ratios or the 24 hours pressure-changes. 

(b) In Fig. 15, there was a general tendency 
for the convective activity to increase with 
the strengthening of the southwesterlies or 
westerlies at the 300 mb level and aloft. 
Likewise, there was a decrease in the convec- 
tive activity with a marked decrease in the 
strength of the southwesterlies and westerlies 
in the upper troposphere. 

Fig. 18 shows the vertical time-section over 
Delhi during the passage of a jet-stream wave. 
It will be observed that the general conclusions 
stated in the above cases for Calcutta with 
regard to the wind, mixing ratios and pressure 
changes, are true in this case also. The moisture- 
content of the air was, however, less than 
that over Calcutta even during the period 


CALCUTTA IN ADVANCE OF A JET STREAM TROUGH 
24 APRIL — 2 MAY 1953 
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Fig. 15. The figures near the wind-arrows and near the base-line below the 2,000 feet winds are mixing ratios. 

The vertical section below 10,000 feet has been drawn on an enlarged scale. The 24 hrs pressure changes refer 

to 03 and 12 G.M.T. respectively. Compare particularly the upper tropospheric winds on each day with the 
convective developments on the following day. 


CALCUTTA IN THE REAR OF A JET STREAM TROUGH 
4 APRIL —II APRIL 1954 
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Fig. 16. Same conventions in plotting as in Fig. 15. Compare particularly the lower and upper tropospheric 
wind systems and lower tropospheric mixing ratios with the corresponding elements in Fig. 15. 
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IN A REGION OF WEAK VERTICAL WIND-SHEAR 
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Fig. 17. Same conventions in plotting as in Fig. 15. Compare particularly the lower and upper tropospheric 
wind systems and lower tropospheric mixing ratios with the corresponding elements in Fig. 15. 


A JET STREAM WAVE MOVES ACROSS DELHI 
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Fig. 18. Same conventions in plotting as in Fig. 15. Compare particularly the general 


ratios with those in Figs 15, 16 and 17. 
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Fig. 19. Mean pattern at 3,000 feet during the spell of 

large-scale convection in and near the State of Assam 

(i.e. roughly north of 23° N between 90° and 95° E). 
Compare with the normal pattern in Fig. 23. 
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Fig. 23. Normal pattern at 3,000 feet in April. Compare 
with Fig. 19. 
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Fig. 20. Mean pattern at 5.000 feet during the spell of 
large-scale convection in and near the State of Assam. 
Compare with the normal pattern in Fig. 24. 


Fig. 24. Normal pattern at 5,000 feet in April. Compare 
with Fig. 20. 
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Fig. 21. Mean pattern at 500 mb level during the spell 

of large-scale convection in and near the State of Assam. 

Same convention in plotting as in Fig. 10. Compare 
with Fig. 25. 


Fig. 22. Mean pattern at 300 mb level during the spell 
of large-scale convection in and near the State of Assam. 
Same convention in plotting as in Fig. 2. Compare with 
Fig. 26 especially with the spceds of the actual winds and 
the horizontal and vertical wind-shear west of 90° E. 


Fig. 25. Normal pattern at 500 mb level in April. Compare 
with Fig. 21. 


Fig. 26. Normal pattern at 300 mb level in April. 
Compare with Fig. 22. 


Fig. 27. The rainfall figures refer to the plain-stations. 

The numerator in each figure refers to the actual rainfall 

(in inches and correct to one place of decimals) while 

the denominator refers to normal rainfall over the 

same area. AA’ is the axis of the jet stream trough at 

300 mb level during the spell. BB’ is the normal position 
of the axis of the trough in April vide Fig. 26. 
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of intense convection. The discussion in section 
2 is relevant in this connection. 


7. Mean patterns during an unusually long 
spell of large-scale convection 


As already mentioned in section 6, the State 
of Assam and the contiguous districts of East 
Pakistan (i.e. roughly cast of 90° E and north 
of 23°N) experienced in association with 
Sy-T situations, an abnormal spell of wide- 
spread thunderstorms, frequent squalls and 
heavy convective rainfall during the period 
27 March 1954 to 14 April 1954! (dates refer 
to the soo and 300 mb charts). Figures 19, 
20, 21, and 22 show the mean streamline 
pattern at 3,000 and 5,000 feet at 09 G.M.T. 
and the mean contour-patterns at the soo 
and 300 mb levels and the partial thermal 
winds for the 700 ~ soo mb and soo ~ 300 
mb levels at 09/15 G.M.T. during the period 
30 March to ı2 April 1954. 

In order to make the charts as representative 
of the mean conditions as possible, no station 
with less than 6 days’ daily data (out of 
the total number of 14 days) was used in the 
computations. This unfortunately resulted in 
no mean wind-data for theregion east of 90° E 
at the soo and 300 mb levels. However, the 
radiosonde-data of Shillong were available 
for all the 14 days of the spell. These and also 
the radiosonde data for all the other Indian 
stations for the same period taken from the 
Indian Daily Weather Reports were specially 
checked by reference to the India Meteorolo- 
gical Department. No data for any station 
were corrected by the present writer from 
synoptic considerations. 

A comparison of these mean patterns with 
the corresponding normal patterns for April,? 
vide Fig. 23, 24, 25, and 26 shows the following: 

(a) With regard to convergence, it is very 
difficult to find any significant difference 
between the mean streamline patterns for the 
spell and the corresponding normal? streamline 


1 The spell continued beyond this date. The present 
study, however, has been confined up to 14 April 1954 
only. 

? These are based on 09 G.M.T. normal pilot balloon 
winds and 15 G.M.T. normal radiosonde heights as in 
the case of May, vide Section 2. 

® No normals are available for Mohanbari (27° 29’ N 
95° or’ E) in the extreme northeast of Assam. It is there- 
fore not known whether there is also not a trough-line 
near that station at 3,000 feet and 5,000 feet, in the 


= 


patterns for the 3,000 and 5,000 feet levels in 
the region of large-scale convection. If at ale 
contrary to the conventional ideas—the stream- 
lines in the mean patterns for the spell show 
an anticyclonic instead of cyclonic curving 
in the region between 90° E and 94° E and 
between 23° N and 26° N where convection 
was most vigorous. It may be added that 

the mean pattern for 5,000 feet level, for 02 

G.M.T. is also very similar to the normal 

pattern for that hour. These have, however, 

not been reproduced here. 

(b) In contrast to the above, there is pro- 
nounced difference between the mean patterns 
for the spell and normal“ patterns for the soo 
and 300 mb levels. In particular it may be 
noted that in the 300 mb level 

(i) the mean jet stream during the spell is 
more concentrated than the normal jet: 
the increase in the horizontal as well as 
vertical shear and the strengthening of 
the wind at every station along and to 
the south of the jet axis may be partic- 
ularly noted. A rough computation shows 
that there has been an increase of more 
than 35 % in the general wind-speed, 
irrespective of direction, north of lat. 
20° N. As will be shown in section 9, 
this increase in wind-speed is equivalent 
to an increase in upper convergence (in 
the rear of the trough line) and increase 
in upper divergence (ahead of the trough 
line) of about 80%, provided other condi- 
tions are the same. 

(ii) The jet stream trough whose axis is near 
86° E on the normal chart (see Fig. 26: 
the partial thermal winds show the 
trough a little more clearly) is much more 
pronounced in the mean pattern for the 
spell and has shifted eastwards. The mean 
upper winds over Calcutta and the mean 
heights of the 300 mb surface over Shillong, 
Calcutta and Allahabad (see Fig. 22) show 
that the region of large-scale convection lay 
in a trough and that the trough-line of this 


normal pattern, as in the mean pattern for the spell. 
Incidentally, the very low speeds of the winds at Mohan- 
bari at 3,000 feet (or knot) and 5,000 feet (04 knots) 
in the mean pattern Figs. 19 and 20 may be noted. 

* The pilot balloon ascents at Delhi did not reach the 
300 mb level on a number of days on which high winds 
were recorded by the radar. The latter observations were 
available on all the 14 days and, as such, may be con- 
sidered as more representative of the mean conditions. 
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large-scale convection in May. The figures in thin type 
over Burma and the Bay and Arabian Sea Islands are the 
number of days of thunder at individual stations. Over 
India and Pakistan, lines of equal frequencies of thunder 
have been drawn. The figures in thick type are the 
average number of squalls per month at individual stations. 


system lay to the east of Calcutta. There 
cannot be any doubt on these two points no 
matter how the contours are drawn. And it 
would perhaps be justifiable to place 
the trough-line somewhere near 89° E or 
90° E as has been done in Fig. 22. It there- 
fore follows that the unusually long spell 
of large-scale convection in and near Assam 
was associated with an unusual strengthen- 
ing of the jet stream and the shifting of 
the trough to an unusual position and its 
remaining there in a quasi-stationary state 
during a whole fortnight. 


Fig. 27 shows the actual and normal rainfall 
at plain stations in various parts of Northeast 
India and East Pakistan between 03 G.M.T. of 
31 March 1954 and 03 G.M.T. of 14 April1954. 
This period of rainfall corresponds exactly 
to the period of the mean 500 mb and 300 mb 
patterns on the basis of the convention adopted 
earlier in this paper for comparing the high 
level patterns with the subsequent weather 
developments. For facility of discussion the 
trough-line AA’ in the 300 mb mean pattern 
for the spell (Fig. 22) and the trough-line BB’ 
in the normal pattern (Fig. 26) have been re- 
produced on the rainfall chart. It is obvious 
from the chart that 
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(i) the region of abnormal excess of rainfall 
lay ahead of the trough-line that is, in a 
region of pronounced upper divergence 
(vide Section 9). 

(ii) the region of abnormally deficient or no 
rainfall lay in the rear of the trough-line, 
which is normally a region of pronounced 
upper convergence (vide Section 9). 


Thus the above study brings out the im- 
portance of upper divergence and conver- 
gence in the production of convective rainfall 
which in a State like Assam can produce floods 
even before the monsoon sets in. 


8. Normal 300 mb flow-patterns in relation 
to normal convective activity 


Fig. 28 shows the normal number of days 
of thunder in May reproduced from the India 
Meteorological Department publication Cli- 
matological Atlas for Airmen (1943 edition, page 
25). On this map have been superpose 

(a) the normal contour-patterns at 300 mb 
level in May shown in Fig. 2. 

(b) the normal number of squalls of 40 
mph or more in May at different stations in 
India and Pakistan as published by various 
workers (P. R. Krısuna Rao 1938, C. RAMAS- 
wamy and K. C. MAJUMDAR 1950, K. P. 
RAMAKRISHNAN and B. GoPINATH Rao 1954). 


An examination of this composite diagram 
for regions north of latitude 15° N, i.e. for 
regions under the influence of the westerlies, 
leads us to the following conclusions:! 


(a) The area east of 83° E which experiences 
more than 8 days of thunder, more or 
less coincides with the area in which there 
is decreasing cyclonic and increasing 
anticyclonic vorticity. The maxima of 
thunder of more than 8 days in the State of 
Assam may to some extent be due to oro- 
graphic effects in that area, but neverthe- 
less the pattern does suggest that the 
vorticity variations contribute very signif- 
icantly to the development of convection. 

(b) There is a remarkable decrease in convec- 
tive activity west of the trough-line AA’ 
which is a region of increasing cyclonic 
vorticity. It is very interesting to see 


1 In normal patterns, the shear variations are less than 
in the daily patterns. The effect of variation in the shear 
on vorticity has, therefore, been neglected in this discus- 
sion, for the sake of simplicity. 
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that Allahabad (25° 08’ N, 79° 07’ E) and 
Nagpur (21° 09’ N, 79°07’ E) which lie 
close to the trough-line AA’ but in its 
rear experience 1.6 and 0.4 squalls respec- 
tively, while Calcutta which lies ahead 
of this trough-line experiences 5.2 squalls. 

(c) There is a marked decrease in the number 
of squalls as we proceed from Calcutta to 
Madras along the coast: Calcutta has 5.2 
squalls, Visakhapatnam 1.4 and Madras 
0.6. This decrease occurs in spite of the 
fact that the entire coastal region is swept 
by moist air in the lower troposphere and 
the depth of the moist air is greater over 
Madras than over Calcutta. The reason 
for this is that as we proceed further and 
further to the south, the wind decreases 
considerably in speed. As will be shown 
in the next section, a pronounced increase 
of the westerly wind with height is one 
of the essential conditions for the develop- 
ment of upper divergence upon which 
depends large-scale convection. And se- 
condly the upper divergence is propor- 
tional to the square of the wind-speed in 
jet-stream zones. Both these conditions 
are, in the mean, less and less fulfilled as 
we proceed from Calcutta to Madras 
along the coast. Consequently large-scale 
convection as manifested in squalls be- 
comes progressively less as we proceed 
towards the south from Calcutta. 

(d) There is a secondary maximum of thun- 
derstorms between 77° E and 80° E north 
of 27°N. This region falls east of the 
trough-line BB’, ic. in a region of de- 
creasing cyclonic vorticity. 

Although the association between the num- 
ber of days of thunder and the normal contour 
pattern is not satisfactory to the west of 75° E, 
the phenomenon of thunder in this region 
shows in general a decrease as we proceed 
from 34°N to 20° N. This seems to be for 
the same reason as the decrease in the number 
of squalls as we proceed from Calcutta to 
Madras namely that we are proceeding from a 
region of strong winds to a region of lighter 
winds in the upper troposphere. It may be in- 
cidentally noted that New Delhi which is 
very near the axis of the mean jet experiences 
the maximum number of squalls (7.0) in 
Northern India during this month. 

With regard to the unsatisfactory association 


between the number of days of thunder and 
the contour-pattern west of 75° E, the following 
tentative explanation is offered: 


(a) In this region, the predominant convective 
phenomenon is the duststorm in which 
thunder is often not heard due to the 
noise of the violent winds. Hence the 
frequency of thunder as reported by the 
observers- cannot give a true picture of 
the convective acitivity in this region: 
frequency of cotivective duststorms alone 
will give a correct picture. Unfortunately 
no statistics of convective duststorms are 
available as the phenomena reported as 
duststorms by the observers include not 
only the convective duststorms but also 
those! caused by steep pressure gradients 
which are not determined by convective 
processes. 


(b) From the study of the individual cases, 
one gets the impression that, in this re- 
gion, there is a greater variability of wave- 
patterns than to the east of 75° E or 80° E. 
Consequently, perhaps, no pattern typical 
of convective situations emerges as the 
normal for the month of May. 


The 300 mb contours for March and April 
have also been similarly studied with refer- 
ence to the normal number of days of thunder 
and of squalls. They also show similar charac- 
teristics but the wave-pattern east of 80° E is 
less marked in these two months. 

A discussion about the maximum of thunder 
in the extreme south of peninsular India is 
outside the scope of the present paper as it 
occurs outside the normal sphere of influence 
of the sub-tropical jet stream. However, it 
may be noted that the relative proportion 
between the number of squalls and the number 
of days of thunder in this region is very differ- 
ent from that in Northeast India where also 
the number of days of thunder is large. This 
would suggest that the mechanism of thunder 
in the southern tip of the peninsula is very 
different from that in the region of influence 
of the sub-tropical jet and that the phenom- 


! These are common in this region and are referred to 
in Indian synoptic literature as dust-raising winds. They 
are popularly known as the loo. It is quite easy to di- 
stinguish them from convective duststorms on the 
synoptic charts. It need hardly be added that only con- 
vective duststorms have been plotted on the synoptic 
charts and diagrams in this paper. 
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enon is probably of a much milder type, re- 
sembling the thunder in Northeast India during 
the southwest monsoon period. 


9. Mechanism of large-scale convection over 
northern India and Pakistan 


It is well-known that the horizontal diver- 
gence in a zonal wave moving without 
change of shape in a broad current is given 
by the equation 


divx v= = 


where u is the zonal wind-speed 
c is the speed of the zonal wave 
£ is the relative vorticity 
f is the coriolis parameter 


and x and y are the rectangular coordinates 
measured in the west to east and south to 
north directions respectively. In the derivation 
of this equation, the w component of the 
wind-speed and the north to south variation 
of the absolute vorticity have been neglected. 

Synoptic experience indicates that the wave- 
patterns in the westerlies over India and 
Pakistan move very slowly in comparison 
even with the wind at the 700 mb level. 
And, as we have already seen, the westerlies 
markedly increase with height above the 700 
mb level in the pre-monsoon period. Hence 
u — c in the above equation is generally positive 
above the 700 mb level and increases with 
height at higher levels. Now let us assume 
that the level of non-divergence over Indo- 
Pakistan is near about the 700 mb! level. 
Let us also assume, as is implied in the above 
equation, that the contribution to the variation 
in the vorticity is mainly due to the curvature 
term. It then follows from the above equation 
that in the eastern half of the trough and the 


western half of the ridge where = is negative, 


v 


there is a steady increase of divergence with 
height above the 700 mb level, the diver- 
gence reaching a maximum at the tropopause 


1 In middle latitudes, this level has been assumed to 
be near about the 600 mb level but we are justified in 
assuming a slightly lower level over Indo-Pakistan in 
view of the variation of wind with height observed in 
tropical cyclones. However, it is not important from the 
point of view of the later discussion in this paper, whether 
the level of non-divergence is nearer 600 or 700 mb. 
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which lies near the 200 mb level over Indo- 
Pakistan in the pre-monsoon period. Like- 
wise, in the western half of the trough and the 


eastern half of the ridge where . is positive, 
Ix 


there must be convergence above the 700 mb 
level which would steadily increase with 
height and reach a maximum at the tropopause. 
This upper divergence and convergence in 
the different parts of the zonal wave must 
result in a corresponding compensation con- 
vergence and divergence below the 700 mb 
level according to the Dines model and such a 
compensation over Indo-Pakistan must be 
even more complete than in middle latitudes 
in view of the well-known fact that the 
pressure-changes at sea-level in the sub- 
tropics in the absence of tropical cyclones is 
much less than the pressure changes in middle- 
latitudes. 
As a result of this upper divergence and 
lower convergence, there would be, in the 
eastern half of the trough or the western half 
of the ridge 
(i) strong upward motion at all heights in 
the troposphere 

(ii) lateral contraction and vertical stretching 
of the air below the level of non-diver- 
gence 

(iii) lateral spreading and vertical shrinking of 
the air above the level of non-divergence. 

Now let us confine ourselves for the time- 
being to the synoptic situations in Northeast 
India and East Pakistan where warm, latently 
unstable air is separated from the upper dry 
westerlies by an inversion. As a result of (1) 
and (ii) above, the latently unstable moist air 
below the inversion will be vigorously forced 
upwards. In this process, it will break the 
inversion lid and shoot up to the levels where 
the energy from latent instability can be 
released. Such a process would obviously 
manifest itself in a rapid development of a 
vigorous squall-producing thunderstorm cell. 
It may be added that the process as visualized 
above is fully supported by observations of 
the actual development of the nor’wester. 

Per contra, in the western half of a trough 
and the eastern half of a ridge where there 
would be upper convergence and lower diver- 
gence, there would be 

(i) strong downward motion at all heights 

in the troposphere 
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Fig. 29. Schematic representation of dynamical processes in Northeast India and East Pakistan in the pre- 
monsoon period. Thick arrows indicate air movement in moist air and thin arrows air movement in 
dry air. Below and very near the inversion in Fig. 29 (a), short continuous arrows indicate forced lifting 
due to orography and long continuous arrows, velocity-convergence in lower troposphere in associa- 
tion with trough-lines etc. Broken ellipses indicate local convection due to insolation and long broken arrows 
indicate nocturnal Bleeker-Andre convergence in the valleys of Assam and along the foot of the 
Himalayas. Note particularly the magnitude of the dynamical effects associated with insolation, lower- 
tropospheric trough-lines etc. in (Fig. 29 (a) in comparison with the dynamical effects associated with 
upper divergence. Note the absence of the inversion in Fig. 29 (b). This schematic diagram also holds 
good with minor modifications, for the andhis in Northwest India and West Pakistan (see text). 


(i 
(iii) 


lateral spreading and vertical shrinking of 
the air below the level of non-divergence 


lateral contraction and vertical stretching 
of the air above the level of non-diver- 
gence. 


In such a development in Northeast India 
and East Pakistan, the moist air below the 
inversion will shrink and spread under the 
inversion. The inversion itself would also 
probably intensify and lower as a result of 
the strong downward motion above it. Con- 
sequently quite stable, fair weather will pre- 
vail. These fair weather conditions have been 
repeatedly observed in Northeast India in the 
rear of troughs, vide discussion in Section 6. 

The dynamical processes discussed in the 
above paragraphs have been schematically 
depicted in Figures 29 (a), 29 (b), and 29 (c). 
The effect of surface-heating and of the 
other lower tropospheric agencies such as 
convergence associated with the trough-lines 
or with the nocturnal sinking of isobaric sur- 
faces in valleys (BreeKER and ANDRE 1951) 


are also shown schematically below the inver- 
sion. These agencies which are of minor 
importance in comparison with upper diver- 
gence and convergence are discussed in greater 
detail in the later part of this section. 

The above discussions are equally valid for 
the rest of northern India and West Pakistan 
also (ie. west of 83° E, north of 20° N) if 
only we take into account one additional 
factor, namely a prior injection of moisture 
below the region of upper divergence. This 
injection usually takes place, as pointed out by 
the earlier workers, by the western disturb- 
ances. The inversion which we have mention- 
ed in the case of Northeast India and East 
Pakistan does not exist in Northwest India or 
West Pakistan, presumably because the in- 
cursion of moist air below the dry upper 
westerlies is a temporary phenomenon and as 
such there is not sufficient time for the devel- 
opment of the inversion. Nevertheless, latent 
instability develops in this semi-moist air which 
means that the moist air is subject to resistance 
to lifting. And this resistance is overcome by 
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the vertical stretching and the strong upward 
motion induced by upper divergence. In the 
circumstances the schematic diagrams in Fig. 
29 (a), (b) and (c) can, in a sense, be also taken 
as valid for the development of large-scale 
convection over any part of northern India 
and Pakistan. 


It should be specially stressed that the 
schematic diagrams are only crude representa- 
tions of the actual dynamical processes and 
that they are not intended to give any quanti- 
tative estimate of the vertical stretching, lateral 
spreading etc. The author’s purpose in present- 
ing them is mainly to show graphically that 
the dynamical effects produced by upper di- 
vergence and convergence are much more 
pronounced than those produced by the lower 
tropospheric agencies shown in Fig. 29 (a). 


If the zonal wave described above is of the 
sinusoidal type, the upper divergence would 
be a maximum not immediately to the east 
of the trough-line (or to the west of the ridge- 
line) but somewhat downstream in the case 
of the trough and upstream in the case of the 


ridge where ae is a maximum. Hence large- 
Ox 


scale convection would tend to be a maximum 
in such regions. This has been actually observed 
on quite a number of occasions in our studies 
and has been schematically shown as Sy-Sin 
in Fig. 9 (a). 

Expressing the above in more general 
terms, if the wave-pattern is of any other 
type, large-scale convection would tend to 
be a maximum in the region where the rate 
of change in curvature of the contours is a 
maximum. Two simple cases of this type are 
shown in schematic patterns Su-T and Su-R 
in Figures 9 (b) and 9 (c). In Su-T, there is 
an abrupt decrease in cyclonic vorticity just 
to the east of the trough-line with little or no 
change in vorticity further downstream. 
Consequently, large-scale convection would 
be practically confined to the area shown in 
that pattern. Likewise in the case of Su-R type, 
there is little or no change in vorticity followed 
by an abrupt increase in anticyclonic vorticity 
downstream in the rear of the ridge-line. 
Consequently, large-scale convection will be 
practically confined to the area shown in 
that pattern. It is interesting to note that a 
large number of cases of Su-R type have been 
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found in our studies confirming the conclu- 
sions arrived at here purely from theoretical 
considerations. 


In the above discussions, it had been assumed 
that the contribution to the change in the 
vorticity by the shear-term was negligible. 
However, this assumption would not always 
be justified. S. PETTERSSEN, Dunn and MEANS 
(1955) have, recently, in another connection 
solved this difficulty, so far as jet stream zones 
are concerned. They have pointed out that 
in the jet-stream zones, the vorticity-advection 
A can be written as equal to 


_y (= + KK) (2) 


x 


where K, is the curvature of the contour K, 
is the orthogonal curvature of the contour 


and V is the wind-speed. 


Now let us see the implications of the above 
equation with reference to the schematic 
patterns Sy-T and Sy-R in Figures 9 (d) and 


9 (e). 
OK, 
In the pattern Sy-T, the term Zu 8 well as 


K,K, will be negative and, therefore, the 
vorticity-advection A will be positive and 
large. As positive vorticity-advection is an 
indication of upper divergence, it follows 
that there will be pronounced upper diver- 
gence and consequent vigorous convection 
ahead of the trough line. In the rear of the 
= will be positive but K, K, 
will be negative since the contours are diverg- 
ing. Hence the vorticity advection A will 
be positive or negative depending upon the 


K, : 
algebraic sum of = and K,K,,. It is, therefore, 
x 


quite reasonable to expect positive values of 
vorticity-advection in the rear of the trough- 
line if the contours are markedly diffluent. 
And since the speed of the wind enters as a 
second power term in the above equation, 
it will be readily seen that, with high-speed 
jets, there can be appreciable upper diver- 
gence in the rear of the trough-line also. 
Further, along and near the trough-line, 
the temperatures would be lowest and this 
would further contribute to the development 
of instability. Consequently, there could be 
development of large-scale convection in 


trough-line, 
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the rear of the trough-line, the severity of 
the convection decreasing as we proceed 
further and further away from the trough- 
line. This is exactly what we find in weather 
situations of the Sy-T type, the westward ex- 
tension of the large-scale convection in the rear 
of the trough-line depending upon the extent of 
fanning of the contours and the wind-speed in 
the jet stream in that region. 

Arguing on the same lines, it can be seen 
that in the case of the Sy-R pattern, the upper 
divergence would be pronounced in the rear 
of the ridge-line but could also be appreciable 
ahead of the ridge-line, the eastward extension 
of the divergence area ahead of the ridge-line 
depending upon the confluence of the contours 
and the wind-speed. In this case, however, 
the higher temperatures along and near the 
ridge-line may slightly counteract the devel- 
opment of instability. These conclusions are 
fully supported by the schematic pattern shown 
as Sy-R in Fig. 9 (e). 

During the course of our investigation, we 
have found that thunderstorms occasionally 
occur ahead of the trough-line in regions 
where the contours show a certain amount of 


. DR : 
confluence. In this case un is negative and 
x 


K,K, is positive but the vorticity advection 
may still be positive if the algebraic sum of 
these two quantities is negative. And if the 
wind-speeds in such a region are high, the 
net positive vorticity-advection, i.e. net upper 
divergence, may be quite large and lead to 
violent convection in spite of the confluent 
contours. 

The effect of the curvature and orthogonal 
curvature of the contours has been discussed 
above on a purely qualitative basis. A quanti- 
tative estimation of vorticity advection alone 
can make it absolutely certain in any individual 
case whether the net vorticity advection is 
positive or not, i.e. whether there is net upper 
divergence or not. Such quantitative estima- 
tion of upper divergence should preferably 
be made at the 250 mb level in view of the 
height of the tropical tropopause but this can 
perhaps be attempted with confidence only 
when more radio-wind observations become 
available over Indo-Pakistan. 

It may be specially emphasized that equation 
- (2) above is of general applicability and as 
such holds equally good for all the schematie 
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patterns in Fig. 9. This implies that, other 
conditions being the same, the severity of the. 
large-scale convection would depend greatly on 
the wind-speed in the upper troposphere. As this 
conclusion has a direct bearing on the problem 
of area-warnings for squalls it will be referred 
to again in section 10. 

It can also be readily seen from equation (2) 
that there can be numerous combinations in 


the values of V?, — and K,K, which can 


correspond to a specified positive value of 
vorticity-advection i.e. very roughly to a speci- 
fied degree of large-scale convection. It would 
therefore appear that the schematic patterns 
in Fig. 9 which represent the types actually 
observed by us on the upper air maps are only 
some of the various possible patterns which 
can lead to large-scale convection. These other 
patterns may come to light as a result of 
further investigations based on more observa- 
tions. 

In the above discussions, we have made a 
fundamental assumption that there is a steady 
increase of the westerly wind with height. 
When, however, this condition is not satisfied, 
there may not be upper divergence and there 
may even be upper convergence if # becomes 
less than c in equation (1). In such situations, 
there would obviously be no convective 
developments. We have met with such situa- 
tions (which, however, seem to be rather 
rare). Vide Fig. 17 and discussion in section 6. 

With regard to the general magnitude of 
upper divergence or convergence, it may be 
mentioned that Petterson and his collaborators 
(1955) have recently pointed out that vorticity 
maximain the range of2 x 10-*sec-! to 5x 1074 
sec? are regularly found at the trough-lines 
in the jet stream over the U.S.A. at the 300 
mb level. Since the troughs in the sub-tropical 
jet are usually less pronounced than those 
in the polar-front jet, we cannot expect such 
high values of vorticity maxima as 2 x 10-4 
sec"! to 5 x 107* sec! at the trough-lines over 
Indo-Pakistan. However, RIEHL (1954, p. 204) 
has shown that a convergence of only 
1x 10°® sec? in moist air can produce 
much bad weather when lasting for several 
hours. These figures, apart from the evidence 
we have already adduced, lead us to think that 
we have in upper divergence and convergence 
over Indo-Pakistan a factor much more po- 
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tent and of a higher order of magnitude 
than lower tropospheric convergence in winds 
which are often light even at 5,000 feet in 
the pre-monsoon period. 

PALMEN (unpublished!) has shown that the 
ageostrophic wind-component across the con- 


tours is approximately given by the following 
expression 


Vna= 
g do [ah 0 (oh\ w IT 
= = . US - + — — 
elf ov,\ Lon \ dt As \on Kr an 
SAS NET Brunt- Confluence- Thermal 
Douglas Diffluence wind 
term term term 
where g is the acceleration of gravity 


f is the coriolis parameter 

v, is the speed of the wind along 
the contour 

Vra is the ageostrophic component 
perpendicular to the contour 

h is the height of the isobaric sur- 
face 

w is the vertical velocity 


is the temperature gradient on the 
isobaric surface 

s is the direction along the contour. 
n is the direction normal to the 
contour. 


The above equation has been derived from 
1p 
P in 
0 ds 
natural coordinates, using the thermal wind 


v, 
the equation of motion er ena = 


3 Ov; : 
equation to evaluate os and the geostrophic 


dz 
: : Ov; dv, 
approximation to evaluate the terms — and —- 


ot os 
In the jet-stream zone to the south of the 


Ur: Ws à 
=n À large and therefore f 5, à 


small. And consequently v,. is large. Further, f 
is progressively less and less in lower latitudes. 
Hence the ageostrophic motion represented 
by v, is greater in the sub-tropical jet than 
in the polar-front jet and greater to the south 
of the subtropical jet than to the north of it. 


jet axis, 


1 The author is specially indebted to Professor E. 
Palmén for drawing his attention to this equation which 
has not yet been published. 
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As ageostrophic motion is an indication of 
divergence and convergence and as the axis 
of the subtropical jet is situated over Indo- 
Pakistan near lat. 29° N in the typical large- 
scale convection months April and May and 
tends to shift to near about 35°N in June 
before it eventually shifts to the north of the 
sub-continent (with the establishment of the 
monsoon), it will be seen that northern India 
and Pakistan lie on the equatorward side of 
the jet stream and close to the jet axis during 
the large-scale convection period. Conditions 
are, therefore, most favourable in these regions 
for pronounced upper divergence and con- 
vergence and for the consequent pronounced 
dynamical reactions on the lower troposphere 
as already discussed in the preceding para- 
graphs. 

We have so far considered only compen- 
sation-convergence in the moist air. Besides 
this, there may be velocity-convergence in 
this air, as for instance along the trough-line 
CA in Fig. 1. The moist air may also be subject 
to upward impulses due to insolation or 
forced lifting due to orography. In the valleys 
of Assam and along the foot of the Himalayas, 
the nocturnal sinking of the isobaric surfaces 
(BLEEKER and ANDRE, 1951) may lead to 
low-level convergence in the moist air late 
in the night and very early in the morning. 
This type of convergence is of a higher order of 
magnitude than that due to katabatic winds 
and as such is a more important factor in 
the development of nocturnal thunderstorms 
in the valleys and along the foot of the Hima- 
layas. All these lower tropospheric processes 
undoubtedly contribute to the development 
of large-scale convection. But they are either 
invariably present, as for instance, orography 
or are very often present in fair weather as 
well as convective weather as, for instance, 
insolation! and Bleeker-Andre valley conver- 
gence or show only slight association with 
convective developments as, for instance, 
velocity-convergence. On the other hand, 


1 The normal time of development of the large-scale 
convection in the plains (P. R. KrisHNA Rao 1938, 
INDIA Mer. Derr. 1944, C. Ramaswamy and K. C. 
MAJUMDAR 1950, Y. P. R. BHALOTRA 1954) cannot be 
explained on the basis of insolation alone as vigorous 
convection develops rather too long after the time of 
maximum temperature. Some other local destabilizing 
factor, which is at present unknown, seems to come 
into play and determine the time of development. 
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we have seen that in regions where latently 
unstable, moist air is present in the lower 
troposphere, large-scale convection seems to 
be invariably connected with upper divergence 
while persistent fair weather is connecte 
with upper convergence. We have also seen 
that this upper divergence or convergence is 
of a higher order of magnitude than lower- 
tropospheric velocity-convergence. We have, 
therefore, to conclude that among the various 
triggers so far discussed for the release of the 
latent instability energy, upper divergence is 
the most important one. 

With regard to the role of cold fronts of 
the middle-latitude type, we would mention 
that among the 80 cases distributed over 
three years which we have studied and in 
which there were destructive squalls of rec- 
ord intensity, we have not seen any case 
which could be attributed to a cold-front. of 
the middle-latitude type. It may also be noted 
that meteorologists who have investigated 
into tropical conditions in other parts of the 
world (RIEHL 1954, p. 237) have, on the basis 
of other evidence, seriously questioned the 
existence of fronts in the tropics—fronts as 
this term is commonly understood. And we 
would add that, in the case of Indo-Pakiscan, 
the great Himalayan range to the north and 
the Karakoram and Hindukush ranges to the 
northwest form a mighty barrier against 
surges of polar air reaching the plains of this 
sub-continent. It therefore seems a priori very 
unlikely that in the pre-monsoon period, espe- 
cially in April, May and early in June, polar 
air even of the modified type could reach the 
plains of northern India and Pakistan in 
sufficient depth as to retain its characteristics 
in spite of the strong surface-heating and pro- 
duce a middle-latitude type of frontal effects. 
And even if such a cold-front were found, it 
would be worth while to see whether the upper 
divergence associated with the 300 mb trough 
of the western disturbance does not have amore 
direct connection with the convective devel- 
opments than with the cold-front at the 
ground. For the present, therefore, we shall 
proceed on the assumption that such cold- 
fronts are very rare and, as such, need not 
enter into a general discussion on large-scale 
convection in northern India and Pakistan. 

The thermal advective processes as seen 
in the backing of winds with height between 


the 700 and 300 mb layers referred to in 
Section 6 should now be discussed in relation 
to the mechanism of large-scale convection. 
In middle latitudes, such a backing of winds 
could automatically be taken as evidence of 
cold-air advection ‘whether or not! it caused 
a drop in temperature at a particular station. 
However, in view of the following reasons, it 
was difficult to proceed strictly on the same 
basis in our present studies: 


(i) The winds were’ probably ageostrophic 
to some extent over Northern India and 
Pakistan. 


(ii) The cold troughs and warm ridges over 
Indo-Pakistan were less deep than similar 
troughs and ridges in middle latitudes. 
Consequently, the local temperature- 
changes in association with these systems 
would be less pronounced than in middle 
latitudes and at times may even be less 
than the errors in tropical radiosonde- 
data at these high levels. 


In view of the above and also of the fact 
that the region of maximum convective ac- 
tivity in trough-areas lay ahead of the trough- 
lines at the 300 mb level and not in the neigh- 
bourhood of these lines themselves where 
the temperatures would be lowest, it has to 
be concluded that although cold-air advection 
(as inferred from wind-observations) does 
occur in a large number of cases and as such, is 
important,” the dynamical effect discussed in the 
earlier paragraphs is very much more important. 

It will, thus, be evident that large-scale con- 
vection over northern India and Pakistan in 
the pre-monsoon period is overwhelmingly 
determined by upper-divergence, i.e. by a dy- 
namical process which begins to operate above 
10,000 feet and attains a maximum at the level o 
the tropopause, namely, near about the 200 mb 
level. 

It would also appear from this study that 


* That cold-air advection does not necessarily imply 
a drop in temperature at a station can be seen from the 
equation for the “individual change’’ of the potential 
temperature of a moving air-parcel. The vertical com- 
ponent of the wind can be such that the “local change” 


2 Ä B 7) 
of potential temperature 1s zero x = 0}. 
dt 


* The importance of cold-air advection in the upper 
levels in the development of instability in the sub-tropics 
has also been stressed by Riehl (1947) in his studies of 
the thunderstroms in sub-tropical U.S.A. in summer. 
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andhis, nor’westers and the majority of the thun- 
derstorms without squalls in northern India and 
Pakistan during the pre-monsoon period are funda- 
mentally the same! phenomenon, the difference 
between the first two being merely due to a 
difference in the moisture-content of the air 
in the lower troposphere and the difference 
between the last two being due to different 
degrees of convection produced by one and 
the same mechanism. 


It is of interest to know the genesis of the 
et-stream waves which so profoundly control 
the large-scale convection over northern India 
and Pakistan. This has, however, to be taken 
up as a separate problem. Nevertheless the 
author would like to record here his general 
impression that there is probably a stationary 
(or quasi-stationary) wave-pattern in the jet 
over Indo-Pakistan with a wave-length of 
3,000 to 4,000 kilometres and with the axis of 
the trough in this pattern somewhere between 
82° E and 87°E (see Fig. 2). In addition, a 
series of waves seem to move across the sub- 
continent producing, as a result, varying 
vorticity patterns over the northern half of 
the sub-continent. That there can be two 
such wave-patterns superposed over the same 
region has already been recognized by the 
earlier workers on the polar-front jet (AROwA- 
PROJECT 1953). It may be further added that 
the stationary (or quasi-stationary) wave 
referred to above does not appear to be 
merely an effect of the contour of the Hima- 
layas. The normal contours at the 300 mb 
level in May and April (Figs. 2 and 26) may 
be seen in this connection. 

As mentioned in Section 1, large-scale 
convection in northern India and Pakistan as 
manifested in the nor’wester and the andhi 
markedly decreases with the establishment of 
the monsoon. The reason for this is that the 
sub-tropical jet-stream shifts to the north of 
the Himalayas with the establishment of the 
monsoon and consequently the upper diver- 
gence associated with the waves in the jet- 


1 In an article on unusually dusty weather at 
Delhi in 1952 and 1953, (S. C. Roy, 1954) has also 
mentioned that the andhi and the norwester are essentially 
of the same origin. His article does not however contain 
any specific synoptic evidence in support of this statement. 
And, further, it is clear from his later statement in the 
same article that, in his opinion, the determining 
factor in the development of the andhi is in the lower 
_ troposphere. 
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stream which are so essential for large-scale 
convection is absent in northern India and 
Pakistan! during the monsoon-period. With 
the withdrawal of the monsoon, the sub- 
tropical jet-stream returns to the south of 
the Himalayas (KoTEswaraM 1954) and there 
is a revival of the convective activity over the 
sub-continent. (A. K. Roy 1949, K. P. Rama- 
KRISHNAN and B. GopinaTH RAO 1954). 
However, on account of the highly mon- 
soonal character of the weather systems over 
Indo-Pakistan, the sea-level isobars and the 
lower tropospheric flow patterns in the post- 
monsoon and cold-weather periods are such 
that most of the sub-continent is invaded by a 
very dry air stream of land origin and there 
are only very occasional incursions of moist 
air from the Bay of Bengal and the Arabian 
sea into the country in association with moving 
pressure-systems. And the moist air which 
thus penetrates into the country has increasingly 
stable characteristics with the advance of the 
year. Consequently large-scale convection over 
northern India and Pakistan is much less 
frequent in the post-monsoon period than in 
the pre-monsoon period and virtually ceases 
in the cold-weather period. 


10. The 500 and 300 mb patterns as tools in 
forecasting large-scale convection in the 
pre-monsoon period 


The soo and 300 mb vorticity patterns over 
northern India and Pakistan (as judged from 
the contour-patterns) are persistent: they move 
very slowly and undergo only moderate 
variations in their configuration during a 
24-hour period. The 300 mb charts for the 
period 28 April 1954 to 6 May 1954 when 
there were fairly adequate wind-data for this 
level for a continuous period of 9 days may be 
cited as an example. The upper tropospheric 
charts also show a high degree of consistency 
with the later convective developments, 
which is supported by the theoretical reasoning 
outlined in the preceding section. In this 


1 The Northwest Frontier Province of West Pakistan 
and the adjoining areas (regions near about 34° N and 
to the north of this latitude) are an exception: they 
continue to be predominantly under the influence of the 
westerlies in the middle and upper troposphere through- 
out the monsoon period (S. P. VENKITESWARAN 1950) 
and consequently continue to experience squalls through- 
out this period (C. Ramaswamy and K. C. MAJUMDAR 


1950). 


56 C. RAMASWAMY 


connection it may be recalled that the 500 
and 300 mb charts refer to 09/15 G.M.T. 
and that the convective developments super- 
posed on these charts refer to the period 
03 G.M.T. of the next day to 03 G.M.T. of 
the following day. The analysis of the high- 
level charts in the qualitative manner done in 
this paper can be easily completed by 20 
G.M.T., i.e. long before convection begins to 
develop in the plains. And it is to be remem- 
bered that convection usually reaches a vigo- 
rous stage only after 09 G.M.T. and at a large 
number of places only after 12 G.M.T. Thus 
the persistency of the middle and upper 
tropospheric vorticity patterns and the com- 
paratively late development of convection, 
have placed at our disposal an effective tool for 
forecasting more than 12 hours ahead, the broad 
regions in which thunder in general and thun- 
dersqualls and convective duststorms in particular 
would develop. Obviously, this new technique 
is still only in its first stage of development 
and can be used in its present form only for 
area-warnings. It should be improved further 
when more radio-wind observations become 
available and further researches are carried 
out on small-scale motions in the middle and 
upper troposphere. It need hardly be added 
that the success of the technique with the 
present net-work of stations will depend upon 
the number of high-level wind observations 
available on any individual day and upon 
the judicious use of high-level radiosonde 
data in regions in which wind-observations 
are inadequate. 

The new technique would, broadly speaking, 
be to identify first the regions where the existing 
supply of moist air is likely to continue or a 
fresh supply is likely to be made available on 
the next day as a result of changes in the 
lower tropospheric flow patterns. Then the 
regions of upper divergence and convergence 
have to be identified. This can be done very 
qualitatively on the basis of the schematic 
patterns presented in this paper, keeping in 
mind that, other conditions being the same, the 
speed of the wind is a factor of major importance 
and that the maximum upper convergence or 
divergence occurs in general in regions where 
the contours undergo the most rapid changes 
in curvature. The 300 mb chart is the one 
where these factors are to be specifically looked 
for but in the absence of adequate data at 


this level, the soo mb chart can be used as 
the next best, keeping in mind the westward 
tilt of the wave-patterns with height. The ver- 
tical wind-shear between the 700 and soo mb 
and between the s00 and 300 mb layers 
(strength and direction of the relative topog- 
raphy winds) is another important factor 
to be considered by the forecaster. The con- 
tribution by the lower tropospheric velocity- 
convergence as for instance along the CA 
trough-line (see Fig. 1) may also be taken 
into account but very much less weight should be 
given to it than to upper divergence and con- 
vergence. The time of development of the 
convection has to be forecast on the basis of 
the additional agencies known at present, such 
as insolation in the plains and nocturnal low- 
level convergence in the valleys, which would 
help the upper divergence field to overcome 
the resistance to lifting of the moist air. And 
having thus decided upon his prognosis, the 
forecaster should continuously be on guard in 
regions of upper divergence even if the lower 
tropospheric analysis does not clearly reveal 
moist air supply because it is in these regions of 
upper divergence that there can be “‘sudden’’ and 
violent convective developments, as there actually 
were, for instance, with the upper divergence 
situations on 13 April 1952 (nor’westers) and 
29 April 1952 (andhis). 


11. Large-scale convection in the sub-tropics 
in other parts of the world 


An examination of the mean streamlines 
and isobars (RIEHL 1954, p. 9 and 10) and 
daily and mean upper air soundings (SERRA 
1943, LOEW 1943, 1945, Roux 1954, WALTER 
GEORGH 1952, KENDREWS 1953, pp. 123, 394, 
398) of the sub-tropics in other parts of the 
world shows that the equatorward side of 
the jetstream lies above latently unstable, 
moist air (NORMAND, 1938, PETTERSSEN, 1940) 
in the following regions in the summer half 
of the year. 

(1) Southeast United States 

(2) New South Wales and Southern Queens- 
land in Australia 

(3) Union of South Africa outside West 
Cape Province, southeast Bechuanaland and 
South Rhodesia south of 20° S in Africa 
and 

(4) Southeast Brazil (south of 21° S), Uruguay 
and northeast Argentina in South Africa. 
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It is also seen that, in many of the regions 
listed above, the moist air in the lower levels 
is separated from the upper dry westerlies 
by an inversion (or a minimum of lapse-rates) 
as in Northeast India and East Pakistan. 

It is clear from the available literature that 
all the above regions experience large-scale 
convection during the summer half of the 
year (EDWARD BROOKS 1951, GRIFFITH TAYLOR 
1917, LOEW 1943, KENDREW 1953, pp. 118, 127, 
Sil, S12, CARVAIHO 1913, LE Roux’ 1953, 
UNION OF SOUTH AFRICA WEATHER BUREAU 
1949, COMMONWEALTH METEOROLOGICAL SER- 
VICES 19$I, SERRA 1938). 

From the mean positions of the equatorial 
trough as given by RIEHL (1954, page 13) we 
can presume that the convergence in this 
trough does not contribute materially even at 
the sea-level to the development of large-scale 
convection in the above regions. We may, 
therefore, leave out the equatorial trough in 
our further discussions. 

Southeast United States—We shall consider 
here only the tornadoes. An examination of 
the climatic Atlas of U.S.A. by STEPHENS 
VISHER (1954) shows that March, April and 
May are the months in which the States east 
of 110° W and south of 35° N get the maximum 
number of tornadoes. It is interesting to note 
that these are the very months in which 
northern India and Pakistan also experience 
severe thundersqualls or convective duststorms. 
It is also seen from Visher’s Atlas that the 
States in more northerly latitudes in U.S.A. 
experience the maximum number of torna- 
does progressively later in the year compared 
to the States in more southerly latitudes, 
suggesting their association with the pro- 
gressive northward shift of the mean jet- 
stream (Arowa Project, 1953). 

It is of course well-known that the upper or 
surface cold-front is of great importance in 
the development of tornadoes and that 80 % 
of the tornadoes occur in association with a 
well-developed parent low at the surface. 
Nevertheless the fact that even 20 % of them 
develop without a surface low or a front or 
even in a region of surface high-pressure 
(EpwarD Brooks 1951) would suggest that 
upper divergence associated with vorticity 
variations in the upper troposphere might be 
contributing to the development of the tor- 
nadoes. 
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FE New South Wales and Central Queensland — 
We shall discuss the convective developments 
in this region only south of 20°S and east 
of 140° E. 

The axis of the mean jet-stream over Austra- 
lia at the 200 mb-level is slightly to the 
north of 30° S in winter and slightly to the 
south of 30° S in summer (GiBBs 1953). It is, 
therefore, justifiable to assume that the mean 
position of the jet for the year as a whole is 
along 30° S. An examination of the isobronts 
of annual number of days of thunder published 
by H. BARKHLEY (1934) shows that there is a 
pronounced maximum of thunderstorms in 
Southeast Queensland (i.e. roughly between 
25°S and 30°S). According to Barkhley, 
this region of maximum number of thunderstorms 
is traversed by squall-lines. 

Among the line-squalls which affect Southeast 
Australia, the spectacular southerly bursters 
deserve special mention. The city of Sydney 
(33° 52° S, 151° 12’-E) experiences the typical 
southerly burster. Fig. 30 shows the normal 
monthly distribution of this line-squall at 
Sydney based on 32 years’ data. (GRIFFITHS 
1917). On the same diagram has been super- 
posed the monthly distribution of the squalls 
(andhis) at Peshawar (34° 01’ N, 71° 35’ E). 
These are based on Dines P.T. Anemograph- 
data for 17 years (C. Ramaswamy and K. 
C. MAJUMDAR 1950). As the northern hemi- 
spheric summer corresponds to the southern 
hemispheric winter, the abscissa for the curve 
starts from July in the case of Sydney and from 
January in the case of Peshawar. As the speci- 
fications of wind-speed defining an andhi 
and a southerly burster are not the same, 


SOUTHERLY BURSTERS 
AND 
ANDHIS 


SYDNEY -33°52'S ISI°IZE 


PESHAWAR —34°02'N 71°37E 


NO.OF SOUTHERLY BURSTERS AND ANDHIS 


JU À SON D J F M A M J SYDNEY 
S 0 N D PESHAWAR 
Fig. 30. The ordinates represent the normal number of 
southerly bursters at Sydney and andhis at Peshawar in 
each of the months of the year. 
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the absolute number of squalls at the two 
stations are not comparable. However, the 
similarity in the monthly trends at the two 
stations is of great interest. These two stations 
have nothing in common except that they 
are at about the same latitude (although in 
different hemispheres) and that they are in 
the vicinity of the jet stream in the upper 
troposphere. In the circumstances, the ob- 
served similarity in the monthly variation of 
the squall-frequencies does, in the light of the 
detailed evidence already adduced in the case 
of Indo-Pakistan, lend support to the view 
that the jet stream may be the common and the most 
important factor in the development of large-scale 
convection in the sub-tropics in both the hemi- 
spheres. 

Union of South Africa and the neighbouring 
regions—According to Gibbs (1953) the mean 
position of the jet stream over Australia may 
also be taken as true of the rest of the southern 
hemisphere. Gibbs has also shown that the 
mean position of the jet stream in the spring is 
slightly more towards the equator than in 
the autumn. From Gibbs’ paper as well as 
from the later analysis of the radar winds over 
Pretoria (25° 46’ S, 28° 14’ E) and over Maun 
(19° 59'S, 23°25’E) by Hormeyer (1953, 
1954) it is quite clear that even southeast 
Bechuanaland and southern Rhodesia are 
dominated by the jet stream throughout the 
summer. 

The seasonal and annual distribution of the 
thunderstorms over these countries as pre- 
sented by J. J. LE Roux (1953) shows a close 
association with the mean positions of the jet 
stream in the different seasons as given by Gibbs. 
Most of the thunderstorms occur on the 
equatorward side of the jet in all the cases. 
And when we remember that these thunder- 
storms are large in number during the period 
October to April and are accompanied by 
severe squalls especially in the summer months, 
when latently unstable moist air is most 
available in the lower troposphere, the associat- 
ion between the jet stream and large-scale 
convection over subtropical South Africa 
becomes even more evident. 

ı Southeast Brazil, Uruguay and Northeast Ar- 
gentina—The well-known Pamperos in these 
regions occur on the equatorward side of the 
jet stream mostly in spring and summer 
(KENDREW 1953, pp. 511 and 512). 


The appearance öf the approaching Pampero 
and the meteorological changes it causes as 
it passes over a station (U.S. Navy DEPART- 
MENT WEATHER SUMMARY OF SOUTH AMERICA, 
1945) bear a strong resemblance to those 
associated with the Indo-Pakistan nor’ wester. 
The Pampero like its Australian counterpart, 
the southerly burster (GRIFFITHS 1917) is 
supposed to be caused by cold-fronts. WALTER 
Grorcn (1952) has, however, recently shown 
that the cold-fronts which cause the Pampero 
are different from the cold-fronts of middle- 
latitudes and that they pass over a station I—3 
hours before the development of the Pampero. 
As the 300 mb trough in such cases would 
obviously be tilted westwards with respect 
to the surface-front, it would be interesting 
to see whether the Pamperos (and also the 
southerly bursters) do not show a more direct 
connection with the vorticity variations at the 
300 mb level (ic. with upper divergence) 
than with the cold-fronts at the surface. 


GENERAL. It will be seen from the above 
that the jet stream is in the right position to 
cause the large-scale convection in the above 
mentioned countries as in Indo-Pakistan. The 
only important difference between these 
countries and Indo-Pakistan is that, in the 
case of the former, the large-scale convection 
continues throughout the summer while in 
the case of the latter (i.e. Indo-Pakistan) 
the large-scale convection abruptly decreases 
with the establishment of the monsoon. And 
for the reasons already fully discussed in 
Section 9 this difference is due to the fact 
that the jet stream itself moves away from 
Indo-Pakistan during the monsoon. And we 
have further shown that in the case of Peshawar 
in the extreme northwest of Pakistan, which is 
more or less continuously under the influence of 
the jet stream even during the monsoon period, 
the monthly variation of andhis bears a striking 
resemblance to the monthly variation of the 
southerly bursters at Sydney, which is at the 
same latitude in the southern hemisphere. It 
would thus appear that the sub-tropical jet 
stream does not produce large-scale convec- 
tion only in Indo-Pakistan: it seems to produce 
similar convection in the sub-tropics all over 
the world, wherever it overruns on its equator- 
ward side, moist air with pronounced latent 
instability. 
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Abstract 


A few suggestions for improving the barotropic model presented in a previous article (BOLIN, 
1955) have been tested. The results show that replacing the geostrophic relation by a more 
general balance equation improves the forecasts over 24 hours in certain respects. . Further- 
more the incorporation of the effects of a stratosphere in a very approximate way gives a 
better description of the behaviour of the largest scales of motion in the atmosphere. 

The application of the balance equation in numerical forecasting is extended to three-dimen- 
sional baroclinic flow. It is shown that the criterion of integrability of the balance equation to 
obtain the stream function from a knowledge of the pressure field has to be satisfied in the 
course of the computations to maintain elliptic character of the equations. 


I. Introduction 


In a recent article by the author (Bou, 
1955, hereafter denoted by III) a number of 
forecasts for 24, 48 and 72 hours with the 
barotropic model were presented. The results 
showed that this model has a definite prognostic 
value up to and possibly beyond three days. 
Some difficulties appeared, however, when 
trying to extend the forecast period in this 
way. For example large anticyclones (and to 
some extent also cyclones) were forecast, which 
had little correspondence in reality and it was 
noticed that in general phenomena of middle 
scale were better forecast than those of a 
planetary scale. An activation of very small- 
scale waves in the course of the computations 
was often observed partly being a result of the 
non-linear character of the forecast equation, 
partly depending upon truncation errors. These 


1 The research reported in this document has been 
sponsored in part by the Geophysics Research Directorate 
of the Air Force Cambridge Research Center, Air 
Research and Development Command, United States Air 
Force, under contract No. AF 61 (514)-648-C, through 
the European Office ARDC. 
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as well as other errors were apparently not due 
to the neglect of baroclinic effects and led 
us to the conclusion that unsolved fundamental 
problems still remained within the frame- 
work of the barotropic model and that some 
further improvements of the forecasts still 
might be obtained without incorporating the 
baroclinicity of the atmosphere. Furthermore, 
such difficulties as the effect of the boundary 
assumptions are preferably studied with the 
aid of the simplest model. A number of 
modifications of the barotropic model were 
discussed in this article. We shall here report 
on the results of incorporating some of these 
modifications one by one into the simple 
barotropic model. 

The attempts to eliminate the geostrophic 
approximation resulted in a more genera 
balance equation between the pressure and 
wind fields (cf. CHARNEY 1955; BOLIN 1955). 
Recently THOMPSON (1955) has generalized this 
balance equation to apply to the more general 
three-dimensional flow of the atmosphere. The 
author will give a somewhat different analysis 
of this problem in the last section of this paper. 
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2. The use of the balance equation between 
wind and pressure in the non-divergent 


barotropic model 


It was suggested in III that replacing the 
geostrophic relation by a more general balance 
equation might improve the barotropic fore- 
casts, in particular in areas with large cyclonic 
or anticyclonic vorticity or in regions with a 
strong deformation field. Representing the 
(non-divergent) wind field by the stream 
function y and the pressure field by ¢ (= gz, z 
being the height of a pressure surface) we 
arrived at 


V2 y+ 2 (PxxPyy > Pry + VP vf = vd (1) 


where fis the Coriolis parameter and v denotes 
the two-dimensional (horizontal) nabla oper- 
ator. It is now our problem to determineyand& 
so that equation (1) is satisfied and the observed 
wind and pressure fields are described in the 
best possible way. It was considered that the 
pressure field is quite well described by the 
contour field of a constant pressure map in 
spite of the fact that the observed winds are 
used as a guide in the analysis. We can then 
determine the corresponding stream function 
field with the aid of (1). It was shown that this 
equation is of the elliptic type if 


Prescribing some boundary values for y we 
can in such cases obtain the interior distribution 
of y by an iterative method similar to Lieb- 
mann’s scheme for solving a Poisson equation 
(cf. ID. 

One important modification in the scheme 
outlined in II turned out to be necessary. 
The non-linear term in (r) should be evalu- 
ated by the finite difference expression (m is 
a scale factor, cf Bouin, 1955) 

Pax Pyy — xy? À 

m! 

I wk [ (vi a TIRE ef Yi, j) x 


x (Yir1,g-1+ Yi-1,j41 — 29,5) — 


= (Pi, jet + Di js — Yin ; — pis)? | 


~ 


Otherwise a systematic error is introduced 
due to the fact that the first and second term 
are evaluated as finite differences with a 
different gridsize. 

It was first decided to investigate the im- 
portance of the non-linear terms in (1) and 
neglect the variation of the Coriolis parameter. 
This is, however, not permissible. It is true 
that the simplified relation 


vy = 77% 


is used for evaluating the vorticity, when 
applying the geostrophic approximation, and 
that both the non-linear terms and the effect 
of the variation of the Coriolis term are 
neglected here. The corresponding y-field is, 
however, not used for obtaining the wind. 
Instead the geostrophic wind v, = — f-1v¢ xk 
is applied. Denoting the stream function of 
the non-divergent part of the geostrophic wind 
by y* we notice that 


frapr + vf cb = vi 


Comparing this equation with (1) we find 
that they are very similar except for the non- 
linear terms. The geostrophic approximation 
has merely been used in evaluating the term 
in (1) depending upon the variation of the 
Coriolis parameter. In view of the fact that 
vf-vpxftvf-+vd is proportional to the 
zonal velocity and thus predominantly positive 
it is clear that a systematic and large error is 
introduced if the variation of the Coriolis 
parameter is neglected in (1) when determining 
y and the wind field. 

Having obtained the initial y-field, the 
forecasting proceeds in the usual manner with 
the aid of the vorticity equation in the form 


“vtpas(opef y) (3) 


Here J denotes the Jacobi operator. At the end 
of the forecast ¢ can be evaluated with the aid 
of equation (1), which merely means solving a 
Poisson equation with proper boundary con- 
ditions. The procedure thus briefly described 
has been tested and the results are reported here. 
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The soo-mb surface is used in barotropic 
forecasts approximately to depict the mean 
motion of the atmosphere. For the same reason 
the soo-mb surface will be used in these 
considerations?. It was found that the criterion 
(2) usually is not fulfilled everywhere on an 
ordinary soo-mb map. Regions where (2) is 
not satisfied on an average constitute 3—5 % 
of the total area in middle or high latitudes 
(based on maps from September and October 
1954). The iterative method given in III does 
not converge if (2) is not satisfied, which was 
readily seen by actually attempting a solution. 
Some modification has to be introduced. Since 
very little is known about methods for solving 
a differential equation which is elliptic in some 
regions and hyperbolic in others, it was decided 
to change the ¢-field slightly to eliminate the 
hyperbolic areas. By this we also prescribe 
that &> — f everywhere. Our main purpose 
here is to get a preliminary idea of the im- 
portance of using the balance equation (1) 
instead of the geostrophic approximation. This 
slight modification of the data field is probably 
not important for getting a first answer to this 
question; certainly not in areas of strong cyclonic 
vorticity. In more than so% of the points 
it merely meant smoothing out small-scale 
irregularities due to imperfect interpolation 
in reading off the data from the map. In no 
point a height value had to be changed by 
more than 15 m, usually only 5 m or less. 


In applying the iterative scheme outlined in 
II it is not sufficient that (2) is fulfilled, but the 
approximate solutions y”) obtained in the 
course of the relaxation must always be such 
that Cl) = v2ym) > — f in order to secure 
convergence. This was continuously checked 
during the computations which is particularly 
important if over-relaxation is used. As a first 
guess we used yl") given by 


v7 (4) 


y) was obtained by solving (4) applying the 
same boundary conditions as for (1) (cf. below). 


1 In justifying the use of the soo-mb surface in the 
barotropic model the idea of an equivalent barotropic 
atmosphere is introduced in which the wind (or pressure- 
field) varies linearly with height. A similar discussion can 
be carried through for the balance equation. 
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It is obvious that p®) satisfies the criterion 
vey) > — fif (2) is satisfied.1 

y must be prescribed on the boundary of 
the region in order to solve (r) for y in the 
interior, i.e. the flow across the boundary 
should be given. We shall here apply the 
geostrophic relation in the following way. Let 
us assume that 


dy 1 dd 
Pa a (5) 


where s denotes a coordinate along the bounda- 
ry and y is a constant, which we determine 
from the additional requirement that v - (k x 
x vy) = 0, Thus 


y= -3 fF - (6) 


S denotes the total length of the boundary. 
Having determined y from (6) we can integrate 
(s) and obtain y except for an irrelevant 
constant. Usually y is very small and the proce- 
dure in effect means an assumption of geo- 
strophic flow across the boundary. The 
boundary assumptions during the following 
forecast will therefore be practically the same 
as in the forecasts made with the geostrophic 
approximation. It should be remarked here 
that this procedure is somewhat inconsistant 
due to the variation of the Coriolis parameter 
as can be seen by integrating the basic equa- 
tions of motion over the area. This will be 
left for a more detailed discussion by Dr 
J. Adem at the International Institute of 
Meteorology in Stockholm. 


Computations have been made for an area 
21x24 gridpoints with a gridsize of 300 km 
at 50° lat. (a stereographic map projection was 
used). Because of the fact that the criterion of 
ellipticity (2) is barely satisfied in some regions 
the rate of convergence is quite slow when 
solving (1). The maximum residual decreases 
by a factor of two in 25—40 iterations (no 


1 It turned out that a better first approximation 
would have been obtained by using 


yy) =- VE nt À vo 


since then the major part of the first harmonic of the 
final solution is obtained in w(t), whereby the rate of 
convergence of the following relaxation becomes better. 
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500 MB 
2 OCTOBER 1954 
03 GMT 


Fig. 1 a. soo mb contours on October 2, 0300 GMT, 1954. The heights are given in decameter. 


500 MB 
2 OCT 1954, 03 GMT 


Fig. 1 b. The streamfunction in unit 1.28 X 10° m? sec-1 

as derived with eq. (1) from the contour field given 

in fig. 1a. The observed winds are given in con- 

ventional notations: a full barb denotes 10 knots and 
a triangle so knots. 


over-relaxation) and it therefore was necessary 
to perform 100—I50 iterations to obtain 
satisfactory accuracy. This takes 40 minutes 
with the Swedish Computor BESK. It is 
quite important not to discontinue the itera- 
tions too early. Already small errors will 
have an appreciable effect in the following 
forecast, since they will effect the first harmonic 
and give rise to an erroneous circulation around 
the area. Since the rate of convergence 
decreases quite rapidly when extending the 
area, the present scheme is not suitable for 
routine forecasting for an area of the size 
used in III and also in sec. 3 of this article. For 
the limited number of forecasts made here the 
procedure worked satisfactorily. 

Three 24-hour forecasts have been made 
using the procedure outlined above. They were 
based on maps from 28 and 30 September and 
2 October, 03 GMT, 1954, and covered an 
area somewhat larger than the one shown in 
fig. 1 b. We shall denote these forecasts by 
BB in the following discussion. Ordinary 
barotropic forecasts using the geostrophic 
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Table 1 
Date 2 VB 3 €B y TB Fe % VB # 
m m m m m m m m 
24 =forecastss OR a | 
2SNSCDI SA NOR are. 0.91 | 0.96 47 39 87 93 104 3 — 8 Io 
BOSSEPE- SA, O3 de = Nevarıne 0.96 | 0.92 30 43 61 52 80 5 — 22 03 
BROCK AHA OS See 0.89 | 0.96 49 38 108 93 98 19 o 26 


24" FORECAST 
500 MB 
2 OCT. 03.00-3 OCT. 03.00 
1954 


+ 


Fig. 2 a. 24-hour barotropic forecast of soo mb contour field from October 2, 0300 GMT to October 3, 0300 GMT, 
1954 using the geostrophic approximation. 


approximation (denoted by B) were also 
prepared for this area. These forecasts were 
exceptionally good (cf. III). This is of some 
importance, since we thereby can study the 
effect of the present modification without 
serious interference with other sources of 
errors. The verification was made inside an 
area covered by 10 x 12 gridpoints and some- 
what shifted towards the eastern portion of 
the area. 


The results are shown in table 1. Here r 
denotes the correlation coefticient between 
observed and computed changes using the 
geostrophic approximation, while rg in- 
dicates the corresponding coefficient for the 
changes computed with the aid of the balance 
equation. e and eg denote the errors in the 
B and BB forecasts respectively and o, and 
og the root mean squares of the computed 
changes, while o is the root mean square of 
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Fig. 2 b. 24-hour barotropic forecast of 500 mb contour 
field from October 2, 0300 GMT to October 3, 0300 
GMT, using the streamfunction given in fig. 1 b. 
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Fig. 2c. soo mb contours on October 3, 0300 GMT, 1954. 


the observed changes. y and yg denote the 
average changes in the two forecasts and x, 


finally the observed change. 


We observe: 

a) The small scale features of the map (900— 
1,500 km) are better forecast using the balance 
equation. This may be quite important in some 
cases. An interesting example is shown in the 
forecast from October 2. Figure 1 a shows 
the initial d-field and figure 1b the corre- 
sponding stream-function y obtained from 
eq. (1). Notice here in particular the over- 
estimation of the winds in the vicinity of the 
low in the Atlantic if using the geostrophic 
approximation. In the B-forecast for 24 hours 
(fig. 2a) this caused the development of 
comparatively high &-values just to the east 
of Iceland. The ridge becomes unsymmetrical. 
The BB-forecast (fig. 2 b) shows considerably 
better agreement with reality (fig. 2 c). It is of 
some interest to point out that this error in the 
B-forecast sf ee still more during the 
following 48 hours (cf. III) and gave quite a 
bad forecast in most parts of the Norwegian 
Sea. The details of the forecasts from Sep- 


tember 28 and 30 were also better forecast 
using y, than was obtained with the geo- 
strophic approximation, which is shown by 
considering the changes of the wind. 


b) The large-scale flow pattern is quite 
sensitive to good accuracy in the solution. 
Systematic errors due to for example an 
inaccurate evaluation of the non-linear term 
in (1) may sometimes cause significant errors 
(as was found by actual computations). This 
would be still more the case if forecasting 
over a larger area and also if extending the 
forecasting period because of the activation 
of small scale features in the flow. 


We may conclude that the use of the balance 
equation instead of the geostrophic relation 
between wind and pressure, improves the 
forecast of certain details of the map that may 
be important. It seems therefore desirable 
to develop a more rapidly converging method 
of solving the balance equation. Furthermore 
the boundary conditions need further con- 
sideration as well as the question of the im- 
portance of the divergence (cf. sec. 3). 


Tellus VIII (1956), 1 


BAROTROPIC MODEL AND THE BALANCE EQUATION 


3. An approximate way of incorporating the 
effect of the stratosphere on the mean 
motion of the troposphere 


One of the most obvious defects of all models 
for numerical forecasting advanced so far is 
the very crude description of the vertical 
structure of the atmosphere. For example, the 
existence of a tropopause is usually completely 
neglected. It is well-known that the high 
vertical stability aloft has a marked influence 
upon the processes in the troposphere. We 
indicated in III how the effect of such a stable 
stratosphere could be incorporated into the 
barotropic or one-parameter model in a very 
simple way. On top of the homogeneous fluid, 
which is supposed to represent mean conditions 
in the troposphere, we place an infinitely deep 
fluid of somewhat smaller density. The ratio 
(x) of the density difference between these 
two layers and the density of the lower fluid 
is chosen in such a way that the variations in 
the height of the interface depicts the variations 
in the height of the tropopause in reality. For 
simplicity the upper fluid is assumed to be in 
rest. By this last assumption the problem 
remains one-parametric. It is obvious that this 
model is extremely crude, but it at least 
permits an average convergence or divergence 
in the troposphere, which is of great importance 
for the study of the longest waves. 

A homogeneous fluid with a free surface 
has previously been studied by RossBy (1945) 
and YEH (1949), and the procedure outlined 
by these authors can be followed in detail in 
order to arrive at the present forecast equations 
(cf. ID). We obtain the following formula: 


(e824) = Jl, D) 6 


where 


2 
A= f (8) 
x gDy 
Here D is the depth of the bottom fluid and 
D, is the average value of D. The model 
studied by Yeh (model B, free surface) is 
obtained as a special case by putting x’ = 1 
and the ordinary barotropic model corresponds 
to x’ = oo. To depict actual height variations 
of the tropopause we shall put x’ = 1/8. 
Before presenting the results of numerical 
forecasts with this model, we shall discuss some 
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principle features of it with the aid of the 
corresponding linearized equations. We assume 
a basic current U, neglect variations in the 
north-south (y) direction and denote the 
perturbations of the interface by D’. (This is 
obviously somewhat inconsistant. U = — 
—g-f-t+dD/dy and is finite and D must 
therefore become zero if going sufficiently far 
north. Permitting a finite extension of the 
wave in the north-south direction would not 


change the results in principle; cf. CHARNEY 


and ELIASSEN, 1949.) We get 
SD’ FD dD’ ID’ 
Ly pee 
ae AN Tri a) 


The phase velocity c and the group velocity 
Cg are given by 


U - pik? 
Tre te 
U+ß/k?+242c/k? 
= 1+42/k2 ey 


B = Of/dy and k = 2x/L, where L is the wave 
length. c and c, are given as functions L in 
fig. 3. The following values of the parameters 
have been used: Ui 15 mjsec, x’ = 1/8, Dy = 
= 8,000 m, g = 9.81 m/sec”?, f and B are 
evaluated at lat. 45°. For comparison c and c, 
for waves in a zonal current governed by the 


r 


307 


Fig. 3. Phase velocity c and group velocity cg as a 
function ot wave length for divergent and non-diver- 
gent model. 
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non-divergent barotropic model are also 
shown. It is clear that only waves larger than 
some 3,000 or 4,000 km are essentially in- 
fluenced by the divergence term. The values 
obtained by Yeh and those given here differ 
considerably. We obtain for example | c, |< 20 
m/sec which is very reasonable in comparison 
with the values given by Yeh. 

Yeh also studied the behaviour of a solitary 
wave in course of time to simulate a blocking 
tidge. We have here recomputed some of his 
pattern with x’ = 1/8. The initial solitary wave 
represented by D’(x, 0) = A exp (— x?/I?), 
where / was chosen in such a way that D’(10° - 
V2, 0) = 1/2, D'(0, 0). Fig. 4. shows: the 
character of the wave after 18 and 36 hours 
and also the changes of such a wave using the 
non-divergent barotropic model (x > ©). It 
is clear that the dispersion has been decreased 
considerably in our present model. Qualita- 
tively Yeh obtained the same result but because 
of the large value of x the effect was only no- 
ticeable for very long waves, indeed. — It is 
interesting to compare these results with those 
obtained by numerical forecasting with the 
aid of the non-linear equation (7). 

In using (7) for numerical computations we 
notice that 9D/dt and D can be replaced by 
dz/dt and z, where z is the height of any 
‘pressure surface in the lower layer of our 
model. We shall again assume conditions at 
soo mb to be representative for mean condi- 
tions in the troposphere and use them as 
initial conditions for our computations. We 
then proceed as follows: 

Let i and j denote the gridpoints in a rec- 
tangular grid covering the area over which a 
forecast is going to be made. 2, and mr, are 


given at time t as well as one time step (At) 
earlier, (r — 1). Here 
: RUES. À 
Ni = & Vi? +fij (11) 
Ji 
my is a scale factor of the map projection used 
andj; denotes the finite difference Laplace 
operator. Then Jj; (n’, 2”) is evaluated at all 
interior points, J being the finite difference 
form of the Jacobi operator. Transforming 
equation (7) into finite difference form gives 
02 A oz; 


5 ij 


I 
v2 ; 1 (yt. 2t 
"Ot mj ot En ulm", 2") 


(12) 


Fig. 4. The development of a solitary wave in an in- 

finite basic current incorporating the effect of the 

divergence above and according to a non-divergent 
model below. 


which can be solved by Liebmann’s method 
of iteration using Ze dt as a first approxima- 
tion. On the boundaries we assume 92"/dt = 0 
in accordance with the assumptions for the 
ordinary barotropic model. Extrapolation in 
time yields 

ayt 


ott+1 — +t-1 Se) 
= a Al; Er (13) 


yj" 1 4s finally evaluated from z7,*1 with the aid 
of (11) in all interior points, while n=; on 
the boundaries. We are then ready for the 
next time step. At the initial time t = 0, we 
have to use uncentered differences in the time 
extrapolation. To get better accuracy, the 
first time step was only half of the value used 
in the following steps. 

We notice that the new vorticity field at 
each instant is evaluated from z (eq. 12); 
which in turn has been obtained by solving 
(12). In ordinary barotropic forecasts, on the 
other hand, the vorticity change is directly 
given by the vorticity advection. To avoid 
cumulative errors in the present procedure it 
is therefore necessary to have somewhat more 
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Table 2 


Date à 
m 


2 | 


En 9, D 
m m 


A IH I a a ST ae  B re EEE 


24"-forecasts 


DOC 45 03 Mae rides 0.84 | 0.73 144 62 212 83 102 150 44 34 
SEDEC 5A OS mie de 0.86 | 0.81 75 IOI 139 105 124 —4 6 | —29 
DSMDECESE 03e cates soe aes 0.75 | 0.84 92 57 130 95 100 At 8 20 
THe DNA Ss 03. PRE. RL: 0.54 | 0.60 106 92 109 78 116 34 6 17 
Mean | 0.75 | 0.75 104 78 148 90 110 
Mean error -— == 66 18 — 
48"-forecasts 
Ses 4S OS. si. = dre 0.67 | 0.60 207 137 292 152 172 153 83 55 
SD EC 234,203... ee rt 0.80 | 0.63 134 144 211 153 179 54 8 | —33 
TIME 51,703... messe 0.58 | 0.73 212 112 250 123 143 153 73 19 
EAMHED IRF, 703: ot ye ee let 0.35 | 0.59 169 120 170 118 143 | —60 17 15 
Mean | 0.60 — == 


0.64 | 


Mean error 


| 231 ‚| 159 | = 


176 | 128 


severe demands on the accuracy of the solution 
of (12) than on the solution of the Poisson 
equation in the ordinary barotropic model. In 
solving (12) the maximum residual was not 
permitted to be more than half of the value 
previously used. On the other hand the rate 
of convergence is more rapid for a Helmholz 
equation than for a Poisson equation. The 
gridsize used in these computations was 300 
km in a polar stereographic projection at 
so° lat. and a time step (At) was one hour. 

Four forecasts have been made for 24 and 
48 hours covering an area of the same size as 
in III. Two of those were chosen to illustrate 
the large scale behaviour of the atmosphere 
(February 17, 1948 and December 13, 1951). 
The results are summarized in table 2, where 
also the results of ordinary barotropic forecasts 
for the same area are given. r denotes the 
correlation coefficient between observed and 
computed changes with the ordinary baro- 
tropic model (again denoted by B), rp is 
the corresponding coefficient for the forecasts 
with the divergence included (in the following 
discussion denoted by BD). ¢ and en are the 
mean errors in these two forecasts and o, and 
op are the root mean squares of the changes 
respectively. o denotes the root mean square 
of the observed changes. y and yp are finally 
the average changes over the verification area 
with the two models and % is the correspond- 
ing observed value. The verification has been 
made inside the dashed line as shown on fig. 
1 in II. In all these cases x’ = 1/10. Later the 
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forecasts for December 8, 1954 were repeated 
with x’ = 1/8. 

It is quite clear that a significant improve- 
ment of the forecasts has been obtained with 
this modification of the barotropic model. A 
significant difference of the correlation coef- 
ficient is hardly obtained on an average, but 
the average errors are considerably smaller if 
using the BD model. Also the root mean 
square of the changes are in much better 
accord with reality as obtained with the BD- 
model. We observe, for example, that the 
ratio op/o, is almost constant, while o,/o,, 
varies considerably. op/ox is, however, less 
than one indicating.that the value for x’ used 
in these computations is somewhat too small. 
Repeating the forecast with x’ = 1/8 for one 
case brought op in close agreement with o, 
but did not change rp and ep appreciably. Too 
few tests have been made to determine the 
best value of x’, but x’ = 1/8 is probably close 
to the optimum value. The mean change over 
the verification area, finally, is considerably 
better forecast with the BD model. All this 
clearly indicates that essentially the largest 
scales of motion are influenced by this modi- 
fication of the barotropic model as is also 
obvious from inspecting the change pattern 
themselves. 

It is of interest to study some of these 
forecasts in detail in particular in view of the 
discussion of dispersion for the linearized model 
that has been given above. 

The synoptic situation of February 17, 1948, 
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Fig. 5 b. 48-hour B-forecast of s00 mb contour field 
from February 17, 03 GMT to February 19. 03 GMT, 
1948. 


has previously been described in detail by the 
author (BERGGREN, BOLIN, Rosspy, 1949). A 
very intense blocking ridge had been formed 
over and to the west of the British Isles and 
was still intensifying on the initial day of the 
forecast. The initial map used in the compu- 
tations is reproduced in fig. 5 a. Of course, 
the few radiosondes over the Atlantic makes 
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17 FEBR 1948, 03 GMT + 48 HOURS —— 


BD — MODEL 


Fig. 5 c. 48-hour BD-forecast of soo mb contour field 
from February 17, 03 GMT to February 19, 03 GMT, 
1948. 


the analysis somewhat uncertain in particular 
in the area of the trough at long. 40° W. The 
conclusions to be drawn are based upon the 
difference of the forecasts made with the 
B-model and the BD-model and are pre- 
sumably valid even if the initial map is changed 
somewhat. Fig. 5b shows the 48-hour B- 
forecast and fig. 5c the corresponding BD- 
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19 FEBR 1948, 03 GMT 


Fig. sd. soo mb contour field on February 19, 


03 GMT, 1948. 


forecast. The actual flow on February 19, 03 
is given in fig. 5 d. It is obvious that the neglect 
of baroclinic processes causes considerable 
errors in both forecasts. Thus the intensification 
of the trough, that is located over Labrador on 
February 19, is completely missing. This in 
turns may be the reason for the fact that the 
anticyclogenesis between Iceland and Southern 
Greenland is not forecast to be as intense as in 
reality. Also over middle Europe some errors 
very likely depend on the neglect of baroclinic 
processes. A comparison of how the blocking 
ridge is forecast by the two models is, however, 
interesting. A marked improvement is here 
obtained with the BD-model. The contour 
value in the centre is forecast to decrease by 
about 30 m, while it remains unchanged in 
reality. The B-forecast, on the other hand, has 
errors of more than 200 m in this area. The 
difference can be interpreted in terms of 
different rates of dispersion as discussed above 
for the linearized model. The scale of the 
solitary wave was chosen to be about the same 
as the size of the ridge on February 17. A 
definite difference in the rate of dispersion 
was obtained with the two models. A closer 
agreement between the one-dimensional linear 
case and the two-dimensional non-linear 
treatment is not to be expected. The similarity 
is, however, suggestive. 
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Also in the forecast from December 13, 1952 
a similar difference in the prognosis of a 
blocking ridge was obtained, but not as 
pronounced as in the case studied above. 


The barotropic forecasts from December 3 
and 8, 1954 were characterized by strong 
anticyclogenesis in the southern part of the 
area which had little counterpart in reality. 
This is the most disastrous error when fore- 
casting with the barotropic model and often 
makes it impossible to extend the forecasts 
beyond 48 or 72 hours. Undoubtedly the fact 
that constant values of y and & are imposed as 
boundary conditions is one main source of 
these errors (cf. II). Some experiments in 
which the values prescribed on the boundaries 
were altered and attempts to forecast for the 
whole northern hemisphere north of the 
subtropic anticyclones (GATES, personal com- 
munication) show that the boundary assump- 
tions are not solely responsible for these errors. 
Forecasts with the BD-model show a marked 
improvement in this respect. Partly this may 
be the effect of a more correct description of 
the dynamics of large scale phenomena in the 
atmosphere, partly it probably is the result of 
the smaller speed of propagation of influences 
of large-scale systems that is valid for this model 
and that prevents the large scale boundary 
effects to spread rapidly into the centre of the 
forecast area. Indeed, the improvement may be 
very large. For example the maximum height 
of the subtropic anticyclone on December 3, 
1954 was forecast to increase from 5,920 m to 
6,180 m two days later with the B-model, 
while the actual height was 5,940 m. A 
maximum height of 5,981 m was forecast with 


the BD-model. 


From the discussion above it seems highly 
desirable to extend the forecasts to the whole 
northern hemisphere and thereby to a large 
extent eliminate the effect of erroneous as- 
sumptions on the boundaries. However, already 
from the few forecasts made so far it seems 
justifiable to draw the following conclusion. 
It seems necessary to include the effect of a tropo- 
pause and its variations into any model for numerical 
forecasting to be able to forecast phenomena of a 
planetary scale satisfactorily. The method used 
here is admittedly very approximate, but still 
seems to give an important correction in the 
right direction. 


72 BERT BOLIN 


4. On the generalization of the balance equa- 
tion between wind and pressure fields 
to three-dimensional flow 


We have previously studied the balance 
equation between the wind and pressure fields 
for two-dimensional non-divergent flow in 
some detail. It is possible to show that a similar 
balance equation is valid also in the three- 
dimensional case (cf. Tompson, 1955). The 
underlying hypothesis is that gravity-inertia 
waves are not of any interest in the study of 
the large scale behaviour of the atmosphere and 
should be eliminated. Thompson then has 
shown that it is necessary and sufficient to 
omit the total derivative of the horizontal 
divergence in the divergence equation to filter 
out these oscillations. This was shown for a 
linearized model of the atmosphere. The same 
is true in the non-linear case without any 
further approximations as will be shown 
below. However, the question still remains, 
if it is desirable to eliminate these oscillations 
under all circumstances. 


The period of gravity-inertia oscillations is 
equal to or less than twelve pendulum hours, 
i.e. 12—24 hours to the north of lat. 30° N 
(in the northern hemisphere). This is about 
one order of magnitude less than the period 
of the large scale planetary waves in middle 
latitudes. On the other hand, one frequently 
observes wave-cyclones with a characteristic 
time-scale of a little more than 24 hours (cf. 
part I, BERGGREN, Bout, RossBy, 1949). 
Furthermore, the speed of propagation of 
the zero mode of the gravity-inertia waves is 
about 300 m/sec., which is at least one order 
of magnitude greater than that of the weather 
systems. The internal waves, however, move 
considerably slower (cf. BOLIN, 1953) and their 
speed is again only a few times larger than for 
example the speed of waves in the westerlies. 
These comparisons indicate that the elimination 
of gravity-inertia waves probably is well 
justified for the zero (or barotropic) mode and 
for the largest scales of motion in the atmos- 
phere. On the other hand, when studying 
phenomena of the scale of ordinary cyclones 
with the aid of multiple-parameter models, in 
particular in cases of rapid development, it is 


more questionable if the internal gravity- 
inertia waves are irrevalent. With these 


reservations in mind we shall for the time 


being accept the hypothesis that the possibility 
for gravity-inertia oscillation to occur in our 
equations may be excluded. 

Let us use a semi-Lagrangian coordinate 
system with ©, the potential temperature, as 
our vertical coordinate. In this way we 
automatically incorporate the assumption of 
adiabatic flow. This transformation is possible 
provided hydrostatic balance exists, 1e. the 
atmosphere should have a stable stratification. 
The equations of motion then are 


ou Ov 
a, — fv = N (14) 
dv OW 
Se Oe op (15) 


where 


(16) 


the so-called Montgomery’s streamfunction. 
The hydrostatic equation is 


Jig lee 
STE 


and the continuity equation becomes 


CA AA ie fone ve 
535) *? (5)- 


u and v are the two components in the x and 
y-axes respectively of the wind vector v in an 
isentropic surface, f is the Coriolis parameter, 
T is temperature, g is acceleration of gravity, 
p is pressure, z is the height of an isentropic 
surface, & = ¢,/cy where ¢ and ¢ are the 
specific heat of air at constant pressure and 
constant volume, respectively. All differentia- 
tions with respect to f, x and y are made 
keeping © constant and v denotes the two- 
dimensional nabla operator in an isentropic 
surface. Let us decompose the wind vector 
into one rotational and one divergent field 


(19) 
(20) 
(21) 


¢ is the relative vorticity. The two components 
of the deformation field A and B are defined as 
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(17) 


(18) 


v= -vyxk+va 
Thus 
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| _ ou ov 

9x dy 

| u av (22) 
2 Ok 


They are functions of both y and «. 
Equations (14) and (15) can be replaced by 
the divergence equation and the vorticity 
equation 
d dv du ou dv 
HU vera (EF 55) 


+fu- fiv - fo = - VW (23) 
ECHN+LEHNV-v=o (24) 


Subscript indicates differentiation. In the linear 
case it was necessary and sufficient to neglect 
the first term in equation (23) in order to 
eliminate gravity-inertia oscillations from our 
system of equations: (17), (18), (23) and (24). 
We shall see that this is also possible in the 
non-linear case. It is clear from (23) that this 
approximation is permissible if 


Vive (25) 
where the inequality sign means at least one 
order of magnitude difference. To the extent 
this inequality is fulfilled gravity-inertia oscil- 
Jations are unimportant. Already CHARNEY 
(1948) pointed out that this is an important 
feature of the large-scale motion of the atmos- 
phere. With the further assumption that ¢<f 
Charney could show that the geostrophic 
approximation is permissible in relating the 
wind and pressure fields to each other. This is, 
however, a more serious restriction and is 
usually not satisfied in moderate or intense 
cyclones and anticyclones. Only eliminating 
gravity-inertia oscillations, i.e. applying (25) 
transforms (23) into the following relation 


fo = (A+? 0?) fut fv I (26) 


The equations (17), (18), (24) and (26) now 
constitute a complete set of equations in x, y, 
@ and t with the dependent variables y, «, Y 
and p. Our problem is to design the proper 
boundary and initial conditions for solving 
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this system of equations and devise a feasible 
(numerical) method for this solution. It is 
furthermore of interest to show that inertia- 
gravity waves are excluded from this set of 
equations. Thereby the maximum speed of 
propagation of waves is presumably reduced 
permitting larger time steps in the finite 
difference time extrapolation.1 

It is interesting to notice that the state of 
the atmosphere is determined with relatively 
good accuracy by the three quantities y, Y and 
p, which can be determined by observations, 
while « only means a comparatively small 
correction to the major part of the wind field 
given by y and is thus of the same order of 
magnitude as the errors of the observations. 
This does not mean, however, that the «-field 
and thus the (v - v)-field is unimportant for 
determining the changes of the state of the 
atmosphere as can be seen from our set of 
equations. Hence first we are faced with 
problem of determining the initial field of « 
from the quantities that are observed. Knowing 
y", WF, p" and a’ at a given instant (Tr) we can 
obtain (9¢/dt)’ and thus (dp/ot)* from (24). 
Extrapolation in time yields p°*1. An approxi- 
mate value of Y*! can be obtained from (26), 
if neglecting « in evaluating the left side of 
the equation, by solving a Poisson equation 
with proper boundary conditions. p°*1is then 
evaluated with (17) and we are back to the 
problem of determining « from a knowledge 
of the three other variables. To do this one 
can proceed as follows: 

Let us introduce the vertical velocity w = 
= dp/dt as a new variable instead of dp/ dt. 
Equation (18) can then be written 


Or ee: Jw 
Poot Vo (3) ES | 47 


Thus a knowledge of a w will permit us to 
evaluate V?x and by proper integration «. 
For simplicity we shall in the present discussion 
neglect the small effect of the variation of the 
Coriolis parameter in (26) and thus use 


fo-2 (A+ BP-2)= 09% (26) 


1 Charney has suggested to solve the simultaneous 
system of equations directly. The question of integrab- 
ility as discussed below is still of interest. 
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This is not necessary but merely done for 
convenience. For the following derivation it 
is useful to replace the vorticity equation (24) 
in our system by another equation also derived 
from the two basic equations (14) and (15). 


By proper differentiation and making use of 


the definitions of A and B according to (22) 
we obtain. 


MAA —f{B+AV-v—f,v —fyut 
aad EN 
xan ala x 


rvvBrfAsBv v-fw+ 
2 


dxdy 


+ ft +2 =0 (30) 


Multiplying (24), (29), and (30) by (¢ +f), 
— A and — B respectively and adding we 
obtain with the aid of (28) 


CU ghey Gas Ge -v+F,+F,=0 
ot 


where 


G==[(C+f)- 4? - BY = pe (32) 


(31) 


( BEN =) or 
F,= - -2 


ox? oy? dxdy ce 


ay (33) 


and 


vf AJ KR) (3a) 


J denotes the Jacobi operator and K = 1}, 
(u? + v2) being the kinetic energy per unit 
mass. With the aid of (27) and (32) we can 
transform (31) into 


Furthermore we differentiate (17) with respect 
to t and obtain 


(36) 


We can now eliminate w or 2W/0t between 
these two equations and obtain 


ZA 22 Nios Phas (2 pe 
3-26 (52) le? 20 à 


‘9 kc, OF 
(3) Se") © 


0 op\ da | 
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We shall restrict the following discussion to 
equation (38) since a knowledge of w will 
enable us to evaluate « with the aid of (27). 
Equation (38) is a second order differential 
equation in w. The right side is a function of 
YW, y and p, but also of «, which is related 
to w by equation (27). One can, however, 
verify that no other terms are of second order, 
than those retained on the left side. Thus the 
character of the equation is directly determined 
from (38). We find that it is elliptic provided 


ep 
(€ (3) = ©) 
We have already assumed that the atmosphere 


should be stable, i.e. 2p/20 < o and in that 
case the criterion of ellipticity becomes 


(39) 


ae 
re ze (40) 
in view of (32). This is the same criterion as 
the one that has to be satisfied in order to solve 
(28) for y as a boundary value problem, when 
Y is known and « is neglected in evaluating 
A and B (cf. Carney, 1955 and Bot, 1955). 
It is easy to show that (40) implies that & > — f 
but not vice versa. 

It was already pointed out that the right side 
of equation (38) is a function of « Now 
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| vx |<| vy| as implied by (25), and thus the 
terms depending upon « are small compared 
to those only being functions of y, Y and p- 
Since the character of the equation is de- 
termined by the left side it should therefore 
be possible to determine w by successive ap- 
proximations. First the right side is evaluated 
by putting « = o and we solve the equation 
(38) for w. The corresponding «-field can be 
used for a better estimate of the right side, etc. 
Probably already the first approximation is 
sufficient for our purpose. 


Equation (38) is a differential equation for 
© that does not contain any derivatives of time. 
It is well-known that the vertical velocities can- 
not be obtained by merely knowing the mass 
and windfields at a given instant if gravity- 
inertia waves can occur in the system, but it is 
necessary to know also the rate of change of 
these fields. The very fact that we can get an 
equation for @ of the type (38) proves that 
gravity-inertia waves are excluded from our 
system of equations (cf. BOLIN, 1953). 


It was pointed out in section 2 that (40) is 
usually not satisfied all through the atmosphere 
in middle and high latitudes, in particularly not 
at higher levels. Thus equation (38) is hyper- 
bolic in such regions. It is well-known that 
boundary conditions cannot be prescribed for 
a hyperbolic equation. In this case, however, 
the equation is hyperbolic only within small 
areas and it is not known to the author in 
which way the boundary conditions can be 
prescribed. Furthermore, even if (40) is 
satisfied everywhere at one instant the question 
arises whether it will remain so in the course 


of timet. It is clearly seen from (31) that G 
may become negative even if being positive 
everywhere at one instant because of the 
existence of the two terms F, and F,. Thus it 
also follows that hyperbolic regions may develop. 
This is a basic difficulty inherent in using the 
balance equation (26) in numerical forecasting. 
The physical meaning of the criterion (40) is 
not clear but it is interesting to notice that in 
the linearized case or in the case of a circular 
vortex it is similar to the criterion of inertia 
(or dynamic) stability. 

It might be worth while to attempt an 
integration and force the criterion (40) to be 
satisfied in the course of the computations. 
It seems to the author, however, that we 
need a better understanding of the physical 
processes involved to use the balance equation 
in treating the three-dimensional flow of 
the atmosphere. 
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Successive Numerical Map Analyses 
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Abstract 


Two series of barotropic forecasts have been computed from numerically analysed maps. 
The first series consists of 31 forecasts and the second one of 4. 

Due to the fact that it is possible to obtain the forecasts directly from observed data, which 
are punched immediately after they are received, it is possible to save both time (3 4 4 hours) 


and personnel. 


The reliability of these forecasts seem to be at least as good as forecasts computed from con- 


ventionally analysed maps. 


I. Introduction 


During the autumn of 1954 the meteoro- 
logical section of the Royal Swedish Air 
Force in cooperation with the Meteorological 
Institute of the University of Stockholm 
undertook a series of routine computations of 
barotropic forecasts. The results of these 
computations (BOLIN 1955) demonstrated the 
desirability of making such forecasts over a 
longer period. For this reason 72 forecasts 
were computed during the periods December 
I—17, 1954, January 17—February 25 and 
April 12—May 24, 1955. 

The results obtained during the first two 
periods have already been presented (BERG- 
THORSSON, D66s, FRYKLUND, HAUG and Linp- 
Quist 1955). In this connection it should be 
mentioned that the final mean values of the 
correlation coefficients obtained were in- 
significantly affected by the introduction of 
the third period, April 12—May 24, 1955. 
These routine computations of barotropic 
forecasts were resumed in the Autumn of 
1955. 

The development of fast electronic com- 
putors does not only make it possible to 


apply numerical methods in forecasting upper 
flow patterns in the atmosphere but will also 
mean a revolution with respect to data han- 
dling and processing in general. In view of 
the very large amount of data that passes 
through a forecasting center it seems that 
automation will be particularly important 
here. 

The development of methods for weather 
map analysis with a machine (BERGTHORSSON, 
Döös 1955, hereafter denoted by NWMA) 
should be looked upon as a step in this 
direction. It is of some interest to see to 
what extent it already now is possible to 
automatise the analysis and forecasting of the 
soo mb surface and the following experi- 
ment was carried out during the period No- 
vember 14—December 22, 1955. 


2. The Analysis 


When computing the analyses, a preliminary 
value is required consisting of a weighted 
mean value of a forecast (12-hour or 24- 
hour barotropic forecast) valid for the same 
time as the analysis and the normal height 
of the 500-mb surface. The first analysis for a 
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Fig. 1. Location of the grid used in the analyses and 

the forecasts. The analyses were made over the whole 

area, while the forecast computations were made over 

the area limited by the inner solid line. The small 

rectangle indicates the area over which the verification 
of the forecasts was made. 


series of forecasts must thus be made through 
conventional means. 

It has proved sufficient to use a 24-hour 
forecast in constructing the preliminary field. 
Hereby the necessity is avoided of analysing 
two maps a day and also of computing an 
extra 12-hour forecast, which is required in 
the case of the 12-hour alternative. No 48- 
hour and 72-hour forecasts were required on 
Saturdays and Sundays. On the other hand, 
the routine numerical analyses and 24-hour 
forecasts were made on these days. In this 
manner it was possible to achieve continuity 
without resorting to conventional analysis. 
The map for 1500 GMT was analysed 
throughout this period. 

The observation data were punched on a 
tape as they arrived. The majority of the data 
had been received by 2400 GMT after which 
the analysis was computed, usually in 30—40 
minutes, depending on the number of data 
received. The forecast computation—which 
required some 90 minutes—was made immedi- 
ately afterwards. All forecast maps were 
completed by approximately 0300 GMT 
(0400 local time), and were then transmitted 
to the various Air Force units. The gain in 
time obtained by computing the analysis 
numerically amounts to some 3 hours. 

It was possible in most cases to eliminate 
erroneous data in the first test, by comparing 
observed heights and winds with the inter- 
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polated ones in the preliminary field. It 
should be remembered, however, that the 
difference between observed and interpolated 
value may occasionally exceed the prescribed 
tolerance without necessarily warranting the 
conclusion that the observed data are wrong. 
The preliminary field may of course have a 
higher or lower degree of exactness. This is 
particularly true in the areas near the bound- 
aries when developments have been strong, in 
which cases the forecast for these areas is less 
exact. It may also be due to the fact that the 
true height strongly deviates from the normal. 

The possibility was exploited of removing 
the influence of erroneous data when the 
observed data were compared with the inter- 
polated values in the completed analysis after 
which the corrected values could be inserted. 


3. The Forecast 


The procedure indicated by CHARNEY and 
Pricries (1953) has largely been followed for 
the integration of the barotropic equation. 

In order to avoid amplification of small 
scale disturbances, describable neither by 
means of the barotropic model nor by the 
grid-system (one grid unit = 300 km at 
50° lat), the initial and thereafter every 
24-hour field have been smoothed. 

The smoothed field is defined by the 


formulas se 

fag hve 

E=C+kmt 
where & denotes the height-value, ¢ the relative 
vorticity and V? the finite difference Laplacian 
operator. 


; I < 
With In satisfactory results have been 
I 


obtained. 

Further, a method was introduced which 
made re-starting possible from the previous 
time-step in the event of a break-down of 
the machine. Thus the field is doubly stored 
for each timestep, and on the computation of 
the immediately following time-step the 
computed relative vorticity is compared with 
the relative vorticity, extrapolated from the 
preceding two time-steps. If in any point 
the difference exceeds a given tolerance the 
machine will stop, after which it will resume 
operations from the immediately preceding 
time step. 
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24-hour forecast 48-hour forecast 72-hour forecast 


Date q | D [u | — + — ——  — — © 
r [or | ov | & [ele r [ox |oy| € |eloc} 7 Or | Oy € 


.06] 0.79] 113| 154] 159 
.0310.77| 84] 185] 169} : 
0.59] 81] 173] 151 
24] 0.53} 118} 207] 167 
.71| 0.60] 123] 173} 135 


14/1I—55 15 | 234] 25 |0.88| 50] 57] 36|0.72|0.88| 84) 99! So 
15/II—55 15 248| 22 |0.89| 47| 46] 27/0.58 0.83} 90} 98} 93 
16/11—55 15 | 244| 33 |0.92| 62] 48] 30| 0.49] 0.83| 80} 118) 87 
17/11—55 15 | 233] 31 | 0.92] 70] 75] 34|0.48|0.71| 75) 127) 93 
18/11—55.15 | 233] 30 |o.87| 60] 58] 31|0.52|0.80| 90} 88) 64 
19/11—55 15 | 209] 24 |0.95| 82] 85] 29] 0.36 
20/II—55 15 | 189] 34 ] 0.88] 63] 97] 56] 0.89 . 
21/11—55 15 | 241] 25 |0.80] 97] 107] 63] 0.65|0.89| 151] 254] 155] 1.03 
22/II—55 15 | 223] 26 |0.97| 116| 146} 52]0.45 - 
23/1155 15 | 218] 25 |0.95| 96| 177| 103|1.07| 0.76] 107| 250] 199| 1.86 
24/11—55 15 | 212| 24 |o.77| 82] oo! 105| 1.28| 0.92| 141] 193] 219| 1.55|0.94| 163| 282| 277 
25/11—55 15 | 224| 27 |0.88| r0o0| 77| 63| 0.63| 0.83] 135] 119| 126] 0.94] 0.89] 146| 168| 141 
26/11—55 15 | 226] 27 |o.82| 51] 58) 41/0.81 
27/11—55 15 | 225] 29 |o.87| O61] 63] 32)0.52 — | — 
28/11—55 15 | 238| 29 |o.64| 57| 53| 63) 1.11] 0.82] 108| 135] 124 
29/11—55 15 | 210] 33 |0.93| 94] 73] 46] 0.49| 0.94| 119| 130] 65 
30/11—55 15 | 205] 27 |o.94| 87| 83] 31|0.36|0.81| ror] 114] 83 
1/12—55 15 161| 27 |o.76| 91] 95| 5610.62] 0.66) 111| 124| 93 
2/12—55 15 | 242] 31 |o.8ıl 76| 72] 57|0.75| 0.69| 121; 72| 88 
3/12—55 15 | 229] 31 |0.97| 100| IOI] 49] 0.45 
4/12—55 15 | 239] 31 |0.95| 64; 67| 55] 0.86 
5/12—55 15 | 266] 28 |o.9gıl 84] 67] 36] 0.43] 0.46] 138] 225] 169 
6/12—55 15 | 230] 31 [0.94] 91] 120] 47] 0.51] 0.81] 148] 287) 169 
7/12—55 15 | 251| 28 |0.94| 98] 155) 80] 0.82|0.91| 125] 262] 148 
8/12—55 15 | 264] 37 |0.55| 60| 113] 71|1.18|0.84| 120] 225) 133 
0/12—55 15 | 231) 29 | 0.92] 73] 123] 77|1.05|0.81| 88) 214] 179 
10/12—55 15 | 252] 29 |o.03| 71} 82] 75] 1.06) — | — | — | — 
11/12—55 15 | 214] 25 |0.80} 64] 65] 50] 0.78 = 
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12/12—55 15 | 235| 21 |0.79| 59| 65| 39] 0.67| 0.65] 99] 119| 1o2| 1.03] — | — | — | — 
13/12—55 15 | 268] 32 [0.79] 73] 61] 48] 0.65]0.70] 122] 108] 1r0| 0.90 

14/12—55 15 | 248] 26 |o.70| 67| 58] 53] 0.79]0.68} 97] 107| 92] 0.94|0.54| ıı5| 203] 167 
15/12—55 15 | 219} 28 |o.67| 70| 76) 53]0.76]0.18} 95] 167] 150] 1.58] 0.12] 131] 359] 356 
19/12—55 15 | 230] 27 |o.83| 56} 88! 59] 1.05]0.50] 61| ııı| 113] 1.85}0.60] go} 129] 138 
20/12—55 15 | 234| 28 |0.74| 51] 50] 56] 1.10|0.76| 103] 117] 84] 0.81] 0.58] 117] 183] 139 
21/12—55 15 | 238] 27 |o.85| 73] 91] 45] 0.62]0.64] 109| 162] 121] 1.11] 0.73} 83] 184] 139 
22/12—55 15 237| 26 |0.63| 80] 95] 90] 1.13}0.56| 70] 127] 101|1.44|0.12| 99] 136| 178 


Mean | 231] 28 |0.85 75| 84! s4l 0.74 0.74] 107] 154| 120 200] 175] 


1.15] 0.62| 123 


4. Results the observed and computed changes, o, and 
The forecasts were made from two series %y are the root mean squares of the observed 
of numerically analysed maps, for the period and computed changes and ¢ is the root mean 
from 1500 GMT November 14 until rsoo square of the difference between the observed 
GMT December 15, 1955 and for the four and computed changes (the mean error). The 
days 1500 GMT December 19—22, 1955. Verification was made over an area shown in 
The first series consisted of 31 forecasts. 48 fig. I. 
and 72-hour forecasts were made practically One may observe that the computed changes 
daily except on Saturdays and Sundays, (0,) usually are larger than the observed ones 
whilst 24-hour forecasts were made every day. (6,3. th “the Gases “Where er Et strongly 
Both these series begin with a conventional a 
analysis of the 0300 GMT map. A 12-hour from unity the average observed and computed 
forecast was then computed as the basis for changes x and y usually are quite different. 
the 1500 GMT analysis. Thus an erroneous rise over a major part of 
The results of the forecasts are given in the forecast area has been obtained from the 
the table. ris the correlation coefficient between computations in those cases. This means that 
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there is no systematic error in the computed 
displacements of the troughs and ridges in the 
westerlies. 


Such errors seem to be considerably reduced 
by incorporation of a total horizontal diver- 
gence in the troposphere picturing the vary- 
ing height of the troposphere (Born 1956). 


The geographical distribution of the reliabil- 
ity of the forecasts has been investigated. The 
correlation coefficients for 24-, 48- and 72-hour 
forecasts have been computed in 27 points. 
The results are shown in fig. 2. 


The unexpectedly high correlation coeffi- 
cient for the 24-hour forecasts in the Atlantic 
area is due to the influence on the analysis 
of the 24-hour forecast from the previous day 
being comparatively large in areas where the 
station network is sparse. The correlation 
coefficient for the 24-hour forecasts is therefore 
not a suitable gauge for the reliability in 
such areas. 


When the actual changes are small, rather 
unreliable forecasts are obtained. This parti- 
ally explains the low average correlation in the 
northern regions of the Scandinavian Peninsula 
and also in regions west of this area. During 
the period mid-November—mid-December 
the observed changes were comparatively 
small here. In this connection it is interesting 
to compare the distribution of the correlation 
coefficients with the mean circulation at the 
500 mb surface for this period. 

In the 24-hour forecasts the best results 
were obtained for an area east of the Baltic 
(> 0.90), whilst the best results in the 48- 
hours forecasts were obtained for an area 
running from Iceland across the British Isles 
and Central Europe up to southern Sweden 
(> 0.75). In the 72-hours forecasts the best 
results were for the region covering the 
British Isles, western France and northwestern 
Spain. 

In the table are also given the number of 
stations (q) used in each analysis and the root 
mean square of the differences between ob- 
served and the analysed heights (D) at the 
stations. These values demonstrate that the 
degree of smoothing of the analysis is of the 
correct magnitude, as per the standard error 
in height observations (NYBERG 1952). In com- 
parison with the numerical analyses where 
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Fig. 2. Geographical distribution of the correlation 
between observed and computed 24-hour changes (a), 
48-hour changes (b) and 72-hour changes (c). 
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12-hour continuity was used (NWMA), the 
smoothing is somewhat stronger. 


As two independent conventional analyses 
of the 1500 GMT maps were not available, 
no comparison was made, such as that pre- 
sented in NWMA. 


An estimate of the reliability of the numeri- 
cal analysis as compared with the conventional 
may be obtained by comparing the reliability 
of the numerical forecasts computed from 
conventionally analysed maps for the same 
season previous year. A certain improvement 
would then be noted. However, this does not 
prove that the numerical method is better than 
the conventional. Firstly, the volume of the 
material is too small to permit firm conclusions 
from being drawn and secondly, an improve- 
ment of the forecasts can be ascribed to the 
smoothing procedure introduced into the 
forecast computations. 


5. Conclusions 


With the present size and orientation of 
the grid the area along the boundaries is 
particularly sensitive. This is partly due to the 
fact that the forecast used in computing the 
preliminary field is inaccurate along the 
boundaries because of the boundary condition 
(96/9t=0). This applies particularly to the 
upper boundary where the height changes 
may be considerable. The error is further 
aggravated by the necessity of extrapolating 
one grid unit along the boundary of the field, 
in compensation for the unit lost during the 
forecast computation. In addition there are 
fewer stations within this area. Variations in 
time of the topography of the soo-mb surface 
will then be damped. This may in turn cause 
abnormally strong gradients. 


Thus a too intense anticyclone was formed 
between the North Pole and the upper bound- 
ary towards the end of the first period, presum- 
ably due to the unrealistic boundary condi- 
tions in the forecast computations. Because 
of the incapacity of the analysis on the follow- 
ing day to correct the erroneous preliminary 
field the anticyclone was further intensified 
in the next forecast. This resulted in a com- 
pletely impossible analysis for this area three 
days later (December 16). 


It is therefore particularly important to 


obtain as many data as possible in the vicinity 
of the upper boundary. The way in which the 
programme is written at present makes it 
impossible to include observations from sta- 
tions outside the boundaries. These data, 
however, are both important and necessary 
and can be exploited. By using nearby stations 
it is possible to sketch the analysis for these 
areas. An imaginary station can then be 
inserted on the boundary. The error thus made 
has been found negligible if the distance 
between the real station and the boundary 
does not exceed 500 km. 


Naturally it is possible to apply this method, 
to insert imaginary stations, inside the bound- 
aries too. 

This is particularly useful when data have 
not been received from areas where the sta- 
tion network is sparse. 


The difficulties in obtaining a good analysis 
for boundary areas will obviously be reduced 
when it becomes possible to extend the grid 
so that the flow across the boundaries will be 
as weak as possible. 

The analysis for the Atlantic and the Pacific 
Ocean have in most cases been satisfactory 
although upper air data only have been used. 
Data obtained from reconnaissance flights 
have been used when received. 

It is possible, however, to improve the 
analysis for these areas by extrapolating the 
height of the s00 mb level from carefully 
chosen surface data. The stratification can 
then usually be so accurately estimated, that 
the error in the extrapolated height does not 
exceed 40 m. As a rule some five such reports 
are obtained from the Atlantic. 

Another means of improving the analysis 
would be the use of the predicted 30-day 
means of the 700 mb surface which are issued 
two times a month by the Extended Forecast 
Section of the U.S. Weather Bureau. How- 
ever, these must first be extrapolated for the 
soo mb level. The simplest method of doin 
this is to add to the normal height of the 
soo mb level 5/7 of the predicted anomalies 
of the heights of the 7oo mb surface. 

Particularly in high latitude regions where 
the normal has great influence (The Pacific 
Ocean) a better result would be achieved 
than when the ordinary normal heights are 
used. 
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It is possible to summarize the results as 
follows: 

The numerical analysis is better than the 
conventional for areas where the station net- 
work is dense due to the fact that the smoothing 
necessitated by observation errors is more 
consistent. 

Where the station network is sparse, as 
for instance over the Atlantic, the numerical 
method gives in most cases comparable results. 
By resorting to extrapolated surface stations 


there will be minor differences only in the 
results obtained from either method. 

It has been found possible to save both 
personnel and time, some 3—4 hours. 

Barotropic forecasts from numerically ana- 
lysed maps are at least as good as those based 
on conventional analysis. In addition the 
numerical analysis is more objective than the 
conventional as the human factor in inter- 
pretation is eliminated. 
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Abstract 


The paper contains a study of the dynamics of two-dimensional straight parallel flow with a 
continuously curved profile. The differential equation for the stream-function of a neutral wave 
disturbance is used as an analytical tool to investigate the properties and the evolution of waves 
having a certain initial configuration of the flow. The behaviour of the waves is explained by 
using an extension of the method applied by HormBoE (1953) in his study of the Kelvin and 
Rayleigh waves. In order to do that, the problem has been approached from a point of view 
which deviates somewhat from the customary treatment. The method is first illustrated by a 
model with a hyperbolic tangent profile (Th-flow), which is shown to have a very simple 
analytical solution for the stationary wave L,. The form of the solution makes it possible to 
interpret it for other wave lengths (L>L,)as a forced stationary wave which is the resultant ofa 
non stationary physical wave and a Rayleigh wave induced by the proper infinite sliding vortic- 
ities at the central level. When this foreign field is removed from the solution, the remaining 
field, which now is a physically possible initial configuration of the flow, will propagate with 
different speed from level to level. As the tilt develops the wave will begin to amplify if L> L,, 
and it will be damped if L < L,. So the stationary wave represents the transition from short 
stability waves to long instability waves as had been anticipated by Fjortoft (1). The method is 
then extended to more general types of profiles. Some general results on the conditions for 
the existence of stationary waves, which were proved or anticipated by RAYLEIGH (1880), 
FJeRTOFT (1950) and HÖILAND (1951), are given a unified treatment. It is finally shown that, 
if a system has stationary waves, they represent the transition from damped to amplified waves. 
If the number of stationary waves is even, both the shortest and the longest waves are damped 
waves, and the intermediate intervals of the spectrum are alternate bands of amplified and damped 
Iwaves. If the number of stationary waves is odd, the longest waves in the spectrum are amplified. 


Introduction 


The dynamics of two-dimensional, straight 
parallel flow in barotropic models with linear 
profiles has recently been studied in detail by 
HoLMB0E (1953). His physical analysis of the 
problem, with the aid of the fundamental 
theorems of circulation and vorticity has 
brought the theory to a degree of clarity 
which had never been attained before. The 
present paper is an extension of the theory to 
the case of flows having a curved velocity 
profile. 


The models considered by Holmboe are 


systems with continuous velocity profiles 
having a number of kinks and straight line 
segments between them (multiple layer Couette 
flow). The simplest of these models is the 
double layer Couette flow: two unbounded 
layers with different values of the vorticity, 
separated by an interface where the velocity is 
continuous, so that the vorticity has a dis- 
continuity equal to the difference between the 
shear of the two layers. When the interface is 
given a sinusoidal deformation, the resulting 
field of motion will generate a neutral wave 
(a wave without change in shape and intensity) 
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moving towards the concave side of the kink 
in the Couette flow profile: the Rayleigh wave. 
In the treatment of this model given in Horm- 
BOE (1953), the field of motion is separated 
into two components, one of which is the 
double-layer Couette flow with the plane 
undeformed interface; the second component 
is irrotational everywhere except in the space 
between the undeformed and deformed posi- 
tion of the interface. All the properties of the 
Rayleigh wave are thus derived in terms of 
the interface deformation and the jump in the 
vorticity at the interface. In the case of three- 
layer Couette flow the same method was 
applied, by considering the total field (for 
arbitrary deformations of the two interfaces) 
as the resultant of three component fields: 
the Couette flow with plane (undeformed) 
interfaces, and two periodic fields which are 
defined by the vorticities in the areas between 
the two positions of the upper and lower 
interfaces respectively. The evolution of the 
field can, therefore, be described by the inter- 
action of the Rayleigh waves which are produc- 
ed at the two interfaces. 


The models we are going to study in this 
paper are of a more complex nature, due to 
the fact that the profiles to be considered have 
a continuous variation of vorticity. The ana- 
lytical treatment of these systems gives rise to 
mathematical difficulties, but the physical inter- 
pretation of the problem and the qualitative 
description of the field can be carried through 
without effort. The whole treatmen is restricted 
to the cases in which the profile of the mean 
flow has inflection points. In all other cases 
the method cannot be applied. The solutions 
obtained through the use of the neutral wave 
differential equation are much too restricted 
to describe any physically possible disturbance 
in these more general cases. Therefore, to 
solve these problems, the stream function of 
the periodic field must be treated as an arbitrary 
function of time, and the vorticity equation 
has to be solved by other methods. 
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1. — Kinematic. We shall consider only 
homogeneous, incompressible (non-divergent), 
inviscid fluids. The types of motion to be 
investigated belong to a rather restricted class. 
The system is unbounded and the trajectory 
of any parcel is entirely contained in a plane 
which will be called its plane of motion. The 
planes of motion of all parcels define a set of 
parallel planes such that there is no variation of 
the motion in the direction normal to them. 
We can describe the motion at any time either 
by giving the velocity vector at each point in 
the plane of motion, or by giving the isopleths 
of the vorticity. In accordance with the field 
theorem, in an unbounded field, one can be 
obtained from the other, so that both represen- 
tions are equivalent.) We shall make use of 
the vorticity representation for reasons that 
will become immediately apparent. In the 
models we are going to deal with all vorticity 
isopleths are assumed to be, at some prescribed 
time, sine curves oscillating around straight 
lines whose direction, being unique for the 
whole plane, will be called the direction of the 
mean flow. 

The motion will be described with reference 
to a right-handed system of mutually per- 
pendicular unit vectors (i, j, k) such that j is 
perpendicular to the plane of motion and i is 
along the direction of the mean flow. The lines 
Z = const. will be referred to as levels at the 
height z. 

For any given level z, there will be a sinus- 
oidal vorticity isopleth, whose value we shall 
designate by Q = Q(z), with the inflection 
points on the line z = const. Let A be any point 
at this level, and B the point of the Q-isopleth 
on the vertical through A (fig. 1). The differenc 
between the vorticities at B and A is given by 
a Taylor’s series as a function of the distance 
€ and the derivatives of the vorticity at either 
A or B. We shall assume the amplitude of the 
Q-isopleth to be so small that only the first 
term in the Taylor’s expansion need to be 
considered in order to have the desired ap- 
proximation. The vorticity at any point on 
the level z can therefore be separated into two 
terms, one of which (Q) is the mean vorticity 


1 There are a few obvious exceptions to this rule: (i) A 
field of translation may always be added without changing 
the vorticity field; (ii) in unbounded one layer Couette 
flow the direction of the flow is not specified by the constant 
vorticity of the flow. 
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at that level, and the other (q) is the departure 
from the mean value. We have, then, for 


each level z 
Q+q=Q-QO% 
q= — QC 


or 


(1.1) 


In accordance with the field theorem, the 
total velocity vector V can be considered as 
the sum of a mean flow U, whose vorticity is 
Q, and a departure from it, v, whose vorticity 
is q. Consequently, if we write 


V=Üli+v 
we have 
Qj= vUxi=-U’j 
and 
gj=’xv 


We can represent the periodic flow either 
in terms of its rectangular components or of 
its stream-function 

v=ui+wk=jxvy, 


(1.4) 


(1.5) 


2. — Dynamics. The physical properties of 
the fluids, listed at the beginning of the previous 
section, have been chosen to give the simplest 
possible dynamics. It is well known that any 


which gives 
V xv =qj=Vpj 


Fig. 2. 


~ 


motion of a fluid with these properties can be 
entirely described by means of a unique dynam- 
ic principle expressing the conservation of the 
vorticity of the fluid parcels along the path. 

In fig. 2 we have assumed the existence of a 
motion of the described type with profile 
of the mean flow U = U(z) represented at the 
left. At the right are shown two sinusoidal 
vorticity isopleths at some arbitrary levels. 
They are considered to be in phase at all levels. 
The circles with the arrows indicate the sense 
of the periodic vorticity q. If we focus our 
attention on the level z,, the velocity field 
associated only with the periodic vorticities 
qo at that level has the configuration of a 
Rayleigh wave (see HOLMBOE 1953, p. I5). 
If there were no other velocity fields except the 
mean flow and the periodic flow of this Ray- 
leigh wave, the parcels would move rhythmic- 
ally up and down with the same periodicity. 
The Q,-isopleth will, therefore, change with 
time in such a way that it will appear as a 
wave moving towards the right, with reference 
to the mean flow at the level z,, with a phase 
speed proportional to the value of the vertical 
velocities at the inflection points (which in 
turn are entirely defined by q,). However, 
the velocity field associated with the periodic 
vorticity g,, at any other level z,, will have 
the same characteristics as the previous one. 
It will therefore produce at the level z a 
transversal velocity component whose value 
will be the intensity at the level z, multiplied 


a 
by the factor e-klz-zıl, where k-- is the 


wave number and L is the wave length of the 
sinusoidal vorticity isopleth. The total intensity 
of the transversal velocity field at z,, which 
determines the motion of the Q,-isopleth, will 
be given, consequently, by the contributions of 
the periodic vorticities from all levels. 

The previous analysis makes evident that each 
layer of the fluid contributes with a different 
“weight” to the total motion of the Q-waves 
at each level. The expression (1.1) shows clearly 
that it is the curvature of the profile of the 
mean flow that, for a given distribution of 
the Q-isopleth deformation, determines the 
intensity of the periodic vorticity field at that 
level. We can therefore speak of any level as 
being more or less “active” than another one 
according to the values of the curvature of 
the mean profile. 
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The integrated effect of all periodic vortic- 
ities within the fluid will change from level 
to level, and we shall have the Q-waves moving 
with different speeds at different levels. The 
crest and trough lines, initially vertical, will 
in general become tilted with time. The result 
of such a process has been analyzed, from a 
general point of view, by Fyorrorr (1950) 
(see also HOLMBOE 1953 p. 9—10), showing 
its relation to the amplification or damping of 
the waves. We shall return later on to this 
question. 


3. — The stationary wave. If the wave motion 
we have considered in the previous sections is 
such that the streamlines coincide with the 
vorticity isopleths, the field of motion will not 
change the configuration of the vorticity field. 
We can then resort again to the field theorem 
and conclude that there will be no local change 
of the velocity at any point. The wave is said 
to be stationary. The Q-isopleth deformation 
also represents, in this case, the vertical displace- 
ments of the fluid parcels, and is related to the 
stream function y by the expression 

D = =» (3.1) 
which leads to the well known differential 
equation of the stationary wave: 


: U” 
y= (e +5) Y 


If we now consider a frame moving with 
the constant speed C with reference to the 
one used above, (3.1) and (3.2) become respec- 
tively 


(3.2) 


(i= C) Cy, (3-3) 
and 


(3.4) 


Equation (3.3) is called the neutral wave 
condition since it expresses the relationship be- 
tween the values of the stream-function and 
the displacements of the parcels in a wave 
propagated without change of shape and amplitude 
with the phase velocity C. 

Before attempting any general discussion of 
equation (3.2) or (3.4) and its conditions of 
solubility, we shall turn now, in the following 
sections, to a detailed consideration of a par- 
ticular case. 
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4. — Flow with hyperbolic tangent profile (Th- 


flow). We shall consider a model having a 


profile of the mean motion which in many 
ways gives the simplest analytical results and 
allows, at the same time, a clear physical 
analysis of the underlying mechanism which 
governs the evolution of the waves. The mean 
speed varies from level to level proportionally 
to the hyperbolic tangent of the height: (the 
fluid is assumed to be unbounded.) 


U= -thz (4.1) 

The mean vorticity has the value —1 at the 
origin and o at infinity. At the level z = 2 it 
has decreased already to 0.1. It follows that 
in any wave motion with this mean flow, the 
most active part, in the sense defined in section 
2, will be concentrated in two relatively shallow 
layers symmetrically located with reference to 
the central level. 


5.— The stationary wave in Th-flow. For 
reasons that will become apparent in section 
10, we shall look only for waves at rest with 
reference to the fluid at the level of the inflec- 
tion point of the mean profile. The periodic 
velocity field will then have a stream-function 
satisfying the differential equation 


w’’ = (k? — 2 sch?)y, 


(5.1) 


with the general solution 


p =[e, (k + thz)e-** + e,(k — thz) e®*] cos kx 
(5.2) 


To avoid infinite values of the field at infinity, 
we take the first term as the solution for z >z,, 
and the second term as the solution for z < 2p, 
where z, is any level with a finite value of the 
height. 

At the level z,, the expressions for both y 
and y’ must satisfy the requirement of con- 
tinuity, which gives 


k? =1, 


regardless of the level z which we select. 
Therefore, the only waves stationary with 
reference to the central level correspond to a 
wave length L = 2x. The simplest case is 
given by 2, = 0, since it will have complete 
symmetry with reference to the central level. 
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Fig. 3. The stationary wave in Th-flow. 


We shall describe in detail this particular case, 
which has the stream-function 


wy =(1+|U]) e-!#! cos kx, (5.3) 
where the amplitude factor ¢ has been omitted.! 

In fig. 3 we have drawn the curve (5.3), at 
point of abscissa nz, showing the variation of 
the periodic field with height. Two typical 
vorticity isopleths above and below the central 
level are shown in the center of the diagram. 
The field of periodic vorticity q is represented 
by the circles with arrows. The Q wave at 
any level, if driven only by the local periodic 
vorticities at that level, would move, with 
respect to the mean flow, towards the concave 
side of the profile. The effect of the periodic 
vorticities at all other levels is to increase the 
speed of the wave relative to the mean flow 
until its relative phase speed reaches a value 
equal and opposite to the corresponding veloc- 
ity of the mean flow. The wave so obtained has, 
therefore, the character of a stationary b-wave 
as defined by HOLMBOE (1953). 


6.— The forced stationary wave in Th-flow. 
In order to find the stationary wave it was 
necessary to adjust the general solution of the 
differential equation for the stream-function 
so as to satisfy the boundary conditions. How- 
ever, by giving up one of the restrictions, 
namely the continuity of w’, it is possible to 
1 The amplitude factor & is still arbitrary, subject only 
to the restriction imposed by the assumption of small 
amplitude of the waves. It can be seen that near the origin 
the field of vertical velocity approaches the value —e, 
whereas the mean flow itself goes to zero. The theory 
breaks down, then, at this level, since in equation (3.1) 
it was assumed, in order to obtain (3.2), that the periodic 
field is negligible as compared to the mean speed. It 
can be shown, however, that equation (3.2) is valid at 
any level. 


get valuable information about the behavior of 
periodic motions with different wave lengths. 

If y’ is discontinuous at some level, it simply 
means that the fluid slides on itself. If the fluid 
is homogeneous, no possible motion within 
the fluid itself is able to set up that type of 
discontinuity. For this reason, all solutions 
containing them have been eliminated from 
the outset. A closer scrutiny of the meaning of 
these solutions shows, however, that they have 
a deeper significance. The solution 


wy =(k+|U])e-*l#! cos kx, (6.1) 
when k +1 represents a periodic motion that 
cannot be produced by displacements alone. 
We shall show, however, that the addition of a 
periodic sliding vorticity of appropriate strength, 
at the level z = 0, to a field produced by displace- 
ments, will suffice to get the desired stationary 
wave. By subtracting from (6.1) the field of 
sliding vorticities we can find the behavior of 
the remaining field for values of k + 1. 

The solution (6.1) which gives a continuous 
field of y, but has at the level z = o the jump 
in y’ whose magnitude is 

Ay’ =2 (1 — k?) cos kx, (6.2) 
is a stationary wave produced by the addition 
of two fields: (i) a continuous field, which we 
shall call y,, whose vorticity arises from the 
displacements of the fluid parcels, in agreement 
with the law q = —U’’¢; (ii) a Rayleigh field 
Wr, irrotational everywhere except at z =0 
where it has the sliding vorticity (6.2), ie. 


vr=(k-kt)e-klzl cos kx. (6.3) 
When (6.3) is subtracted from (6.1) we get 
wy = (k-+ + |U]) e-*l#l cos kx. (6.4) 


The stationary wave (6.1) whose analytical 
representation is now 


YP=Yot YR, 


will be called a forced stationary wave, to signify 
that it is produced by some foreign field which 
is artificially introduced into the system. The 
removal of the Rayleigh field yp will change 
the field of motion in such a way that the Q- 
waves no longer remain stationary, but the 
knowledge of the form of the remaining field 
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Fig. 4. The b-wave in Th-flow. 


will allow us to get information about the 
subsequent evolution of the field. This field is 
analogous to the b-wave in skew-symmetric 
triple layer Couette flow (see fig. 10a in 
HOoLMBOE 1953). We shall therefore call it the 
b-wave in Th-flow. 


7.— The b-wave in Th-flow. In fig. 4 is 
shown an initial state of Th-flow with forced 
stationary waves for two values of the wave 
length. In order to fix ideas the values k = 2 
and k = 1/, have been selected as typifying 
“short” and “long” waves, respectively. 

A direct inspection of fig. 4 shows that the 
field of motion associated with the periodic 


vorticities of short waves needs the support of 


additional sliding vorticities of the same sense, 
i.e., such that they provide a field of motion 
which speeds up the Q-waves at all levels. 
It is obvious that the removal of the sliding 
vorticities will have the effect of slowing 
down the Q-waves at all levels, producing a 
relative displacement of the waves such that 
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the crest and trough lines, originally vertical, 
become tilted with the same sense as the shear 
of the mean flow. In accordance with the 
Reynolds-Fjértoft’s energy transfer principle 
that type of motion will produce a damping 
of the wave. This result can also be directly 
obtained by a direct analysis of the field of 
motion of the tilted wave. In fact, if we focus 
our attention on a Q-wave at any level in the 
upper layer we see that it is damped by the 
periodic fields associated with all Q-waves 
below that level, and amplified by the periodic 
field associated with all Q-waves above. On 
the other hand, if we consider a Q-wave at 
any level in the lower layer, we see that it 
is damped by the field of motion associated 
with all Q-waves above its level, and ampli- 
tied by the periodic motion associated with the 
Q-waves below. The net result is, then, a 
damping of all Q-waves. The mechanism 
is identical with the interaction of the two 
Rayleigh waves in triple layer Couette flow 
(see for instance fig. toa in HOLMBOE 1953) 


/ 
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The long waves have a component field y, 
associated with the g-vorticities, such that it 
must be opposed by additional sliding vorticities 
of opposite sense to become stationary. The Ray- 
leigh field necessary for the forced stationary 
wave has the effect of slowing down the 
Q-waves at each level. The removal of this 
field will produce, therefore, a speeding up of 
the waves at all levels. The same type o 
reasoning applied above shows that the net 
result is an amplification of the wave, in agree- 
ment with the energy transfer theorem. 

The same result can be obtained analytically 
by calculating the instantaneous speed of the 
Q-waves at each level after the removal of 
the Rayleigh field. The transversal displace- 
ments, 

yp 
S (ii 


in the forced stationary wave represent the 
deformation of the Q-isopleths. Although the 
Q-waves initially propagate with different 
speeds C(z) at different levels, the fact that 
there is no tilt at the initial time will secure 
that the amplitudes of all the Q-waves are 
initially stationary, so we may use the neutral 
wave condition (3.5) at this moment, i.e. 


ae. 


U-C() 


From the two expressions for € we get 


Cie tee (7.1) 
ap 
When the values (6.1) and (6.3) are introduced 
for y and pr, we get 
k?-ı 


© UF eau) 


which is represented in fig. 4 for the values 
je =arand kee Ale 

It can be seen that C has the same sign as 
U when k >1 (short waves), and opposite 
sign when k< 1 (long waves). At great 
distances from the central level, C approaches 
asymptotically the value 


C,= (7.2) 
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Fig. 5. The a-wave in Th-flow. 


8.— The a-wave in Th-flow. The solutions 
we have analyzed so far lead only to a very 
special type of motion, namely one with initial 
vorticity isopleths 180° out of phase above and 
below the central level. Let us now examine 
the evolution of the field if at a given instant 
the vorticity isopleths are in phase throughout 
the field. In fig. 5 we have represented schemat- 
ically the problem we are considering. We 
have taken the values of € and & so that the 
layer z > 0 is identical with fig. 4. If we focus 
our attention on the Q-wave at any one level, 
we can immediately see that the propagation 
which would be caused by its own periodic 
vorticities acting alone is supported by the 
fields associated with all the vorticities in the 
same layer, but is opposed by the fields 
associated with the vorticities in the other 
layer. The motion as a whole has thus the 
characteristics of an a-wave in triple layer 
Couette flow as defined by HOLMBOE (1953). 
A comparison with fig. 4 shows at once that 
the a-wave has also at each level a relative 
motion towards the concave side of the profile. 
The speed of propagation will be, however, 
less than the corresponding speed of the b-wave. 


To compute the instantaneous speed of an 
a-wave, we have to know first the isolated 
effect of each of the two layers. We shall 
consider an auxiliary model in which the mean 
motion is 


U= -thz for z>0, 

Ue —z for z <0. 

The solution for the stationary wave is given by 
p=(k+U)e-* cos kx, for z>o, 


y =kekcos kx, for z<o, 
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Fig. 6. The wo-wave in Th-flow. 


where the constants have been so adjusted as 
to give bounded and continuous values of w 
at all levels. The continuity of y’ requires 


yp’ (0) = -1-R?=k? 
I 


The value k = = V2 gives, then, the stationary 


4) 


wave of the system. For all other values of k 
we shall have only forced stationary waves 
with the sliding vorticities at the central level 


Ay’ = (1 — 2k?) cos kx, 
corresponding to the Rayleigh field 
2k=1 

k 


2 


x“ 


e~kl=l cos kx. 


var 


When this field is subtracted from the total 
field, the remaining field, 


vom (=r Ue cos Riz 06.2 =0, 


(8.1) 


I 
Wy =— ek? cos kx, for z<o, 


2k 


is the instantaneous field associated with the 
4? 
D: 

It is obvious that the same instantaneous field 
of motion would be obtained by considering a 
Th-flow throughout, but with a field of ¢ 
which at that moment is identically zero in the 
lower layer and has the value of the b-wave 
in the upper layer. The instantaneous phase 
speed of the y-wave at each level is found to be 


u(k-2e- 1) 
2 


k-U 


vorticities q = — U’’C, where = - 


Cie) = (8.2) 


which is drawn in fig. 6. 
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At great distances from the central level, C 
approaches asymptotically the value 


Ey cine 
2 


sem nat rare 1:3) 


We have therefore been able to isolate the 
effect of each layer. A simple inspection of fig. 
4, fig. 5 and fig. 6 shows that 


Ya + Yo = 2%, 
or 
Wa = 2Yo — Wy = — Ue-*lzl cos kx, (8.4) 
where (6.4) has been introduced for y», and 
(8.1), (8.2) for y. Similarly 
Ur: 
vhs 


Therefore, the instantaneous phase speed of 
the a-wave at each level is, for all values of k, 
less than the corresponding speed of the mean 
flow. Consequently, all the a-waves wili devel- 
op a tilt in the same sense as the shear of the 
mean flow. 


9. — Comparison of Th-flow with skew- 
symmetric three-layer Couette flow with shear in 
the central layer. The Th-flow has a close re- 
semblance to skew-symmetric three layer 
Couette flow with no shear in the outer layer 
analyzed by HOLMB0E (1953, section 11). In fig. 
7 both types of flow are shown for comparison. 
The dispersion diagram for the Couette flow 
system has been reproduced from fig. 11 a in 
HoıMmBoE (1953). (Only the stationary b-wave 
appears in the diagram.) The spectrum is divided 
into two regions, one of short stable waves 
(L < L,) and the other of long unstable waves 
(L > L,). The phase speeds in this diagram are 
the phase speeds of the upper and lower 
component waves with reference to the fluid 
in the respective outer layers. The correspond- 
ing dispersion diagram for the Th-flow at the 
initial moment would represent the instan- 
taneous phase speeds of the wave—telative to 
the Auid—at great distances from the central 
level. To find these values we have subtracted 
from (7.2), (8.3) and (8.5) the speed of the 
mean flow at infinity (—1). 

The two dispersion diagrams show a striking 
similarity, especially for large values of the 
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Fig. 7. Dispersion diagrams for Th-flow and three layer Couette flow. 


wave lengths. This similarity suggests rather 
strongly a close analogy between the dynamic 
behavior of the two systems. It is reasonable 
to expect that the symmetric waves in Th-flow 
will behave similarly to such waves in Couette 
flow also for some time beyond the initial 
moment. The significant analytical difference 
between the two systems, however, is that in 
the Th-flow the waves will be subject to 
internal deformations which cannot be re- 
presented by fundamental modes. 


10. — Conditions for stationary waves in flow 
with arbitrary curved profiles. The technique 
which has been used for Th-flow can now be 
extended to more general profiles. In order to 
simplify the analysis, let the fluid be bounded 
by two parallel walls which intercept the z- 
axis at the points of ordinates o and h. The 
neutral waves in this system must be represented 
by a stream-function which (ji) is continuous, 
(ii) has continuous derivatives at all levels, (iii) 
vanishes at both boundaries. It will follow that 
the function U(z) must, in turn, satisfy very 
special conditions in order that such solutions 
can exist. 


ie) - 2rd LE 
Let us consider a function defined by 
, , Y 
a = (Uy — U'y) (10.1) 
U 
whose derivative is 
a! = ky? + (7 a) (10.2) 


From its definition (10.1) it follows that the 
zeros of x coincide with the zeros of y. This is so 
even if y and U vanish simultaneously at a 
level z = 2 , since in that case we should have 


Là 


lim «= - lim — y? = 


On the other hand, in view of (10.2), «’ is 
positive at all levels, and therefore « must have 
a discontinuity between each pair of zeros of 
y. This discontinuity can only be provided by 
the vanishing of U at some level where y #0. 
This is the well known Rayleigh’s theorem. 

The same property of the function « shows 
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Fig. 8. 


that U must have a zero between any two zeros 
of y. Therefore the number of zeros of y 
between the boundaries must be at most one 
unit less than the number of zeros of U. It 
follows that a profile with only one zero cannot 
have a stationary wave with a nodal plane. 

If the level U =o is not a nodal plane 
(y #0), U’’ must also vanish. To prove this 
it is enough to consider the function yy’ 
whose derivative, in view of (3.2), is given by 


LA YA © ih [/ 
op (Ba) tt 


If U=o with U’ being zero, the above 
expression and all the higher order derivatives 
of it become infinite. This implies that the 
function py’ is discontinuous, which violates 
the conditions for the stationary wave. The 
zeros of U must be, therefore, inflection points 
of the profile. 

In equation (3.4) the expression in paren- 
thesis cannot be positive for all z. For in that 
case y’’ would have the same sign as y at all 
levels, which is incompatible with the vanishing 
of y at both boundaries. It follows that the 
curvature and the speed of the basic flow cannot 
have the same sign at all levels. 

If we apply these considerations to the case 
of skew-symmetric flow (i.e. a low with a 

rofile which is skew-symmetric with reference 
to the level half way between the boundaries) 
we can see that all the properties we have 
derived here are in complete agreement with 
what was found by HOLMBOE (1953) for three- 
layer Couette flow systems. In fig. 8 b we have 
reproduced fig. 11 b of Hormpoer (1953). The 
Couette flow profiles of typeI can be considered 
as limiting cases of curved velocity profiles 
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having the property . = 0 at all levels. No 


stationary wave is possible in this case. The 
Couette flow profiles of types II and III can be 
considered as limiting cases of curved velocity 


11 


profiles having the property u” at all levels. 


Couette flow profiles of type III, with only 
one zero of U, have only one stationary wave. 
Couette flow profiles of type II, with three 
zeros of U, have two stationary waves, one 
of them with a nodal plane at the central level. 


11. — The sliding vorticities in the forced sta- 
tionary waves. In the treatment of Th-flow we 
were able to show that stationary wave (for 
k = 1 in Th-flow) represents a transition from 
short damped waves to long amplified waves. 
We can now generalize this result. 

If we consider again equation (3.2), the 
presence of an inflection point at the level 
where U = o implies the existence of two 
independent solutions, y, and y, which— 
together with their derivatives—are continuous 
and bounded for all values of z. The only 
admissible solution is, however, that which 
vanishes at both boundaries. If there is such a 
solution it means that the linear combination 


yp = Ay, + By, 


when suitable values are introduced for A and 
B, satisfies the conditions 


y (0) =y (h) =o, (11.1) 


where o and À are the ordinates of the bounda- 
ries. This type of solution exists, if at all, only 
for very special values of k which are called the 
eigenvalues of the problem. In this section we 
shall designate them by k,, and the correspond- 
ing solutions (eigenfunctions) by y, 1. e.: 


uw 

pa Rp ps = 0. (ar) 

For other values of k, any solution which 
vanishes at o will not cut again the z-axis at 
z = h. The Sturm theory for the differential 
equations of the form (3.2) shows how this 
intersection moves along z when we let k either 
increase or decrease from its eigenvalue k, 
(fig. 9). The eigenvalues k, are, therefore, 
“isolated”. By means of the Sturm’s theorem 
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it is easy to show that the symmetric flow of 
type II referred to at the end of the previous 
section has in fact two stationary waves. 

A modified form of the Sturm’s theorem 
will now be applied in order to find how the 
systems behave when the wavelength of the 
periodic motion is such that k#k,. As it was 
stated above, equation (3.2) has, in the absence 
of singularities, two independent solutions, y, 
and we, which are analytic for all z. These 
solutions can always be chosen in such a way 
that y, vanishes at z = o and yw, vanishes at 
z=h. We can get, therefore, a continuous 
solution (full line in left part of fig. 9 c) which 
has the values 
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and satisfies the boundary conditions, but has 
a jump in the derivative at the level z,. As 
we did in section 6, we can again consider this 
jump as produced by a Rayleigh field which 
is Laplacian everywhere except at the level 24. 
This field can always be represented by the 


expression 


pr=or, (11.3) 


where dr is a combination of expotentials. 
When wp is subtracted from the total field y, 
we get a function y, which satisfies the bound- 
ary conditions and has continuous derivatives 
at all levels, and it is defined by the relation 


P= Yo + pr. (11.4) 


The actual value of y, can be directly ob- 
tained when we replace (11.4) for y in (3.2). 
We can easily find that it must satisfy the non- 
homogeneous differential equation 


i Ces ee 
ip = (#+5) Wp = OT PR: (11.5) 


The solution of this equation satisfying (11.1) 
can always be found, by means of the Green’s 
function defined by the homogeneous equa- 
tion, provided that k + k.. 

We shall consider first a profile having only 
one stationary wave of wave number k, with 
a stream function of the periodic field satisfying 
(11.2) and the boundary conditions (11.1). For 
any other wave number k#k, we can always 
find a b-wave whose stream-funcion y, satisfies 
(11.5) and (11.1), ie. such that it becomes 
stationary by the addition of the Rayleigh 
field in (11.4). If we now multiply (11.2) by 
y, and (11.5) by y,, and subtract, we get 
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or, after integration between the boundaries 
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The integral on the left is always positive (due 
to our convention in the selection of the cell), 
while the integral on the right, for the types of 
profiles we are considering (U’/U<o) is 
always negative, so we can write 
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We see, then, that waves with wave numbers 
k>k, would require a Rayleigh field of the 
same sign to make them stationary, whereas 


for k<k, the required Rayleigh field must 
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have opposite sign. In other words, short b- 
waves need the help of a periodic field of the 
same sense in order to propagate with the 
relative phase speed required for the stationary 
condition, which in turn means that the b-field 
alone would propagate with less speed. The 
crest and trough lines, initially vertical, will 
develop tilt in the same sense as the shear of the 
mean flow. The same argument applied to long 
b-waves shows that they will develop tilt in 
Opposite sense to the mean flow. The character 
of the stationary wave as a transition between 
short damped waves and long amplified waves 
is thereby proved. 

If the system under consideration has two 
stationary waves, the same reasoning shows that, 
starting with very short waves and decreasing 
the value of k, every time we cross an eigenvalue 
of k the Rayleigh field which is necessary to 
keep the wave stationary changes its sign. The 
long waves (for k smaller than the smallest 
eigenvalue) and the short waves (for k larger 
than the largest eigenvalue) will behave in the 
same way, since both need the help of a Ray- 
leigh field of the same sign to become stationary 
(damped waves). The intermediate range will 
be of amplified waves, since they must be 
opposed by a Rayleigh field of different sign 
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in order to remain stationary. We have there- 
fore obtained a generalization of the results 
obtained by HOLMBOE (1953) concerning skew- 
symmetric three-layer Couette flow with oppo- 
site shear in the central and outer layers. 

The last result can be easily generalized to 
the case of profiles with any number of zeros. 
If the number of eigenvalues is odd, the shortest 
and longest waves will be damped waves. If 
the number of eigenvalues is even, the shortest 
waves are damped and the longest waves will 
be amplified. The intermediate ranges of the 
wavelengths are divided in regions which 
alternatively will correspond to damped and 
to amplified waves. 
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A Study of Oceanic Fronts’ 


By TOWNSEND CROMWELL and JOSEPH L. REID, Jr., Scripps Institution of Oceanography, 
University of California 


(Manuscript received 3 August 1955) 


Abstract 


Oceanic fronts and frontal layers are defined. Temperature data obtained at two fronts in 
the equatorial Pacific are presented and a particular circulation, considered basic to the existence 


of a front, is proposed. 


Vertically isothermal warm water overlies the stable frontal layer. This suggests that vertical 
mixing in the isothermal water maintains the frontal layer in the same way that wind-induced 


vertical mixing maintains sharp thermoclines. 


The data suggest that the depth of vertical mixing is different in the warm and cool water 
masses adjacent to one of the fronts. It is suggested that the change in vertical mixing regimes 
contributes to the horizontal gradients of properties across this front. 


Introductory Remarks 


Atmospheric “fronts” are studied daily by 
weather forecasters and the term is in wide- 
spread popular use. This paper is concerned 
with an analogous occurrence in the ocean. It 
is the phenomenon often called “current rip” 
or “tide rip” in English literature on the subject, 
although these terms have been applied to a 
number of different phenomena. Recent ob- 
servations in the central tropical Pacific revealed 
a subsurface structure near “current rips” 
similar to that at atmospheric fronts and led us 
to adopt meteorological terminology. 

The majority of investigations into the 
nature of oceanic fronts have been by the 
Japanese. Their interest centers around the fact 
that fish and whales are abundant along fronts 


(UpA, 1938, 1953, 1954). Upa (1938) has 


1 Contribution from the Scripps Institution of Ocea- 
nography, New Series No. 842. 

Work initiated by the Pacific Oceanic Fishery In- 
vestigations of the U.S. Fish and Wildlife Service and 
completed under the auspices of the Office of Naval 
Research and the California Cooperative Oceanic 
Fisheries Investigations. 


published a paper including a history of ob- 
servations of fronts and many descriptions of 
visual manifestations of the phenomenon. 

A surface front is often visible as a narrow 
band of agitated wavelets. This band may 
follow a sinuous path, and, if there is floating 
debris in the area, it concentrates along the 
front. There is frequently a sound, apparently 
produced by breaking wavelets, and which, on 
occasion, has been reported as a “roaring 
noise”. UDA (1938) reports a front audible at 
a distance of three miles. 


Definitions 


We know of no definition of an oceanic 
front. It seems necessary to start with one, but 
it will be inadequate and will be expanded by 
discussion and specific examples. 

A front will mean a band along the sea 
surface across which the density changes 
abruptly. Our measurements across fronts are 
of temperature by a recording thermograph. 
Throughout this paper, which deals with fronts 
in the open ocean where salinity variations are 
small, we will interpret abrupt temperature. 


Tellus VIII (1956), 1 


A STUDY OF OCEANIC FRONTS 95 


N 


] 
Hy 
[] 
que HER 


Fig er: 


SE HET 
Kr 6 OLA 
4 ER Sr ORO COMES 

où 440 5x SON 
Ÿ ENS annul 


Examples of surface temperature gradients in the central tropical Pacific. Ship speed during recording, 


about 9 knots so that each hourly interval (marked along the top) represents about 9 miles. Frontal system 
corresponding to trace a. appears in cross section in Fig. 3 and that corresponding to trace b., in Fig. 4. 


changes as indicative of abrupt changes in 
density. 

Perhaps the fuzziest part of the definition 
lies in the degree of abruptness of the density 
gradient required for a front to exist. Only 
qualitative remarks can be made about this 
requirement. 

Fig. 1 shows thermograph records which 
illustrate sea surface temperature gradients 
encountered in the central tropical Pacific. 
Generally the gradient of surface temperature 
in this region is of the order of degrees per 
ten or hundred miles as can be seen in the 
figure. At fronts, however, the surface tem- 
perature gradient is often of the order of degrees 
per hundredth or tenth miles. Examples of the 
gradients encountered at these fronts in the 
central Pacific appear on the thermograph traces 
in the figure. 

Besides density, temperature and salinity, 
additional properties may change abruptly 
across a front, including current velocity and 
concentrations of dissolved and suspended 
materials, sometimes with striking changes in 
water color and transparency. Only a change in 
density is required by the definition, however, 
for the existence of a front. 

From the front, a stable layer inclines down- 
ward, separating a warm water mass of low 
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density above, from a cooler one of greater 
density below. This transition zone will be 
called the frontal layer, a feature that isanalogous 
to that occurring in association with atmos- 
pheric fronts. There is an upper, warm boundary 
and a lower, cool boundary to the frontal layer. 
We wish to emphasize that the frontal layer 
extends in three dimensions and to refer to 
certain characteristics of the temperature 
gradient across it. In particular, an important 
characteristic is the clarity with which the warm 
and cool boundaries are defined by the BT 
traces. For example, on bathythermograph 
traces 61 and 62 of Fig. 2, the warm boundary 
is clearly defined by a sharp curvature maxi- 
mum on each trace, whereas the cool boundary 
is defined less clearly by minor curvature 
maxima. At some fronts the cool boundary is 
more sharply defined than it is here. 


Temperature distributions near fronts 


1. Bathythermograms were obtained close to 
a front in the equatorial Pacific on October 27, 
1952, by G. I. Murphy from the Pacific Oceanic 
Fishery Investigations’ hydrographic vessel 
SMITH (Fig. 2). The front was clearly visible 
as a band of agitated wavelets about 150 feet 
wide, oriented east-west. The two “cool” 
traces, 58 and 59, were obtained 9 miles and 
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Fig. 2. BT traces on the cool side (58 and 59) and 
on the warm side (61 and 62) of a front. BT trace no. 
60 is in the front. The warm boundary (w. b.) and 
cool boundary (c. b.) of the frontal layer are indicated. 
This frontal system appears in cross section in Fig. 3. 


150 feet south of the front, respectively. Trace 
60 was obtained at the northern edge of the 
zone of agitated wavelets. Traces 61 and 62 
show the thermal structure 3 and 5 miles 
north of the surface front. It is apparent that 
the horizontal temperature change is much 
greater across the front than it is for many 
miles in either the warm or cool water masses 
separated by the front. There is an abrupt 
density change across the front. 

We feel that the character of the vertical 
temperature gradient on the warm side of the 
front is significant. The warm water isvertically 
isothermal and the frontal layer begins abruptly 
at the bottom of this isothermal layer. This 
suggests that there is intense vertical mixing 
in the warm water above the frontal layer and 
also the further possibility that this mixing plays 
a part in maintaining the frontal layer in a 
manner that will be described in a later section. 

Enough BTs were obtained near this front 
to represent the thermal distribution in cross- 
section (Fig. 3). The warm and cool boundaries 
of the frontal layer are indicated using heavy 
solid lines. The thin lines are isotherms. Regions 
where water is nearly isothermal vertically are 


indicated by shading. The hypothesized basic 
water circulation, discussed below, is indicated 
by arrows. % 

The drawing in the vicinity of the surface 
front may need some explanation. BT trace 
No. 60 (Fig. 2) shows water of 72F at 288 feet 
at the northern edge of the frontal zone. Water 
of this same temperature occurs at 98 feet 
on BT trace No. 59, obtained ten minutes 
earlier and approximately 300 feet to the south. 
Thus, the 72F isotherm rises 190 feet over a 
horizontal distance of approximately 300 feet. 
The actual slope of this isotherm is at least 
190 ft./300 ft., or approximately 2 feet of rise 
per 3 feet horizontal distance. 

Because of the vertical exaggeration in the 
section (760: 1) the 72F isotherm appears to 
rise vertically beneath the surface front be- 
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Fig. 3. Vertical temperature distribution across front at 
1°12/N, 120°11/W, on October 27, 1952. Isotherms for 
each °F, thin lines. Frontal layer enclosed by heavy lines. 
Water nearly isothermal vertically, stippled. Hypoth- 
esized basic circulation indicated by arrows. The solid 
arrows indicate water involved in cell-like motion in 
this plane. Vertical exaggeration 760: 1. Surface thermo- 
graph trace across this front appears as Fig. 1, a. 
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Fig. 4. Vertical temperature distribution across front at 1°50’N, 172°00°W, on February s, 1950. Isotherms for 
each °F, thin lines, and for each .2F, thin dashed lines. Frontal layer enclosed by heavy lines (where clearly defined 
by BT observations) and heavy dashed lines (where diffusely defined). Water nearly isothermal vertically, stip- 
pled. Hypothesized basic circulation indicated by arrows. The solid arrows indicate water involved in cell-like 
motion in this plane. Vertical exaggeration 760:1. Surface thermograph trace across this front appears as Fig. 1, b. 


tween 288 and 98 feet. Successive isotherms of 
higher value, found at considerable depth on 
BT trace No. 60 in the frontal zone, are found 
closer to the surface or not at all on BT trace 
No. 59 in the cool water. Thus, the nearly 
vertical line extending upward from 290 feet 
is made up of several isotherms, some of which 
intersect the sea surface at the front. 

The deep penetration of warm water locally 
beneath the surface front is based on a single 
BT observation, No. 60 (compared to those 
nearby). It suggests the existence of a steeply 
inclined “frontal wall”. Closely spaced BTs, 
obtained recently by R. A. Schwartzlose from 
the Scripps vessel HORIZON, at a front off 
Lower California, revealed no such feature. 
More observations are needed to determine 
whether or not a frontal wall is a real feature 
of some oceanic fronts. 

2. Data were obtained near a front en- 
countered Feb. 5, 1950. A cross section showing 
the thermal structure at this front was drawn 
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by CroMWELL (1953) and has been redrawn 
by us in the light of additional information 
about fronts discussed in this paper and to show 
greater detail (Fig. 4). Typical BTs near this 
front are reproduced in Cromwell’s paper. 
There was no observation immediately at the 
position of the front. 

This frontal layer is more complex than the 
other. It appears to be “broken”, 30 miles north 
of the surface front. In this region the frontal 
layer is spread upward and downward so that 
it occupies a large part of the surface layer. 
Further, the warm and cool boundaries of the 
frontal layer through this region are not 
discernible. The frontal layer is defined again, 
farther to the north. 

As at the front in Fig. 3, the warm water 
adjacent to this front is nearly isothermal and 
the frontal layer begins abruptly at the bottom 
of this isothermal layer. Data obtained near 
two additional fronts, not presented here, 
revealed this same basic thermal structure. 
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The circulation near fronts 


There are, to our knowledge, no published 
measurements of subsurface currents near oce- 
anic fronts. No such measurements were made 
at the fronts shown in the cross sections in this 
paper. However, it seems to us that a reason- 
able speculation can be made as exhibited by 
arrows in the sections. 


Surface convergence along fronts : 


Many observers have noted that foam, kelp, 
and floating debris sometimes collect in con- 
spicuous bands. It seems almost axiomatic that 
there is a component of surface flow directed 
into these zones. Feebly swimming plankton 
observed in fronts (BEEBE, 1926; Upa, 1938) 
are probably concentrated there by the con- 
vergent surface flow. The physical data, in 
particular the abrupt temperature changes 
across fronts, are consistent with this idea. 
Surface waters of different temperatures, by 
virtue of their different origins, are brought 
into juxtaposition by the convergent flow. 


Sinking at fronts: 


Continuity demands that there be a flow 
to compensate for the horizontal convergence 
that occurs along fronts. This could be ac- 
complished by sinking beneath the surface 
front and such flow has been observed, strong 
enough to sink fishing nets at fronts off Japan 
(Upa, 1938). Recently we observed a similar 
occurrence. 

On November 30, 1953, the authors placed 
a crude drag near a front off San Diego harbor 
to observe its motion. The drag consisted of 
two mutually perpendicular 3 ft. by 3 ft. planes 
of 1/32 inch sheet iron, attached by 40 feet 
of line to a surface float (corks), pole and flag. 
Upon reaching the front, the pole, which 
previously had lain tilted, became suddenly 
vertical and sank slowly as though pulled from 
below. Although we remained in the area for 
an hour and watched for the flag, it was not 


sighted again. 
Divergence at depth: 


Given these two features of the flow, there 
must be horizontal divergence at depth. If 
there is a main thermocline which remains 
stationary, the divergence must occur some- 
where in the surface layer. 

These three features of flow, 1) convergence 
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Fig. 5. Effect of wind-induced turbulence on the verti- 

cal temperature gradient in the ocean, Curve 1: Before 

a strong wind. Curve 2: After a strong wind has blown 
for some time. 


along the front, 2) sinking and 3) divergence 
below, are shown in the cross sections. Water 
which sinks in the front and which must have 
a temperature somewhere between that of the 
warm and cool water masses probably sinks 
to intermediate depths in and below the frontal 
layer as shown by the dashed arrows. In Fig. 3 
a jet is indicated downward beneath the surface 
front, but, again, it should be mentioned that 
more data are needed to determine whether or 
not a deep penetration of warm water is a real 
feature beneath some oceanic fronts. In the 
vertically isothermal water masses, where sta- 
bility would not oppose them, cell-like circula- 
tions are drawn. We have no information 
about where water rises to the surface to com- 
pensate for that which sinks along the front. 


An hypothesis concerning the maintenance 
of frontal layers 


The maintenance of frontal systems involves 
the problem of the maintenance of the thin, 
stable, frontal layer in spite of diffusion across 
the layer which should act to destroy it. A 
process that may maintain the frontal layer is 
suggested by a different, but, presumably, 
related phenomenon, i.e., the maintenance of 
sharp thermoclines by wind-induced turbulence. 


Some effects of wind-induced turbulence on the 
vertical temperature gradient : 


Curve 1 of Fig. 5 pictures a typical vertical 
temperature gradient in the ocean. Curve 2 
Tellus VIII (1956), 1 


ASS UID Y, FOFNO CE ANTEC 


10 
u 

D8 JULY, 195] 
a 62°00'N, 33°00'W 
w6 

uw 

a 

24 

= 

=2 

= 


O 


16 17 
DATE 

à 0 
Cp 100 
pi 100 
= 200 
A 200 
= 300 
Es 300 
= 400 
= 400 

50 30 
: x 146 


RRO NALS 99 
19 20 
(GCT) 
Q 
100 
200 
300 
400 
Rome SO 
30 
AO 50 


Fig. 6. Selected BTs illustrating ‘‘confined mixing’’ in the surface layer. Wind and BT data obtained at weather 
patrol station by USCGC UNIMAK. Times of BT observations indicated in wind speed diagram by arrows. 


presents the gradient in the same region after 
a strong wind has blown for some time. This 
picture is virtually the same as figure 2 in 
H. ©. Misc. 15360 prepared by SCHULE, et al. 
(1952). 

The strong wind increases turbulence at the 
sea surface and, presumably, downwardthrough 
the surface layer. This turbulence alters the 
vertical temperature gradient in three respects 
pertinent to the present discussion: 1. Mixing, 
which is made more intense by the turbulence, 
reduces or eliminates the gradient throughout 
a layer adjacent to the sea surface. This is shown 
schematically in the figure between the surface 
and depth A. 2. The gradient is intensified in 
a layer immediately subjacent to the mixed 
layer, as illustrated between depths A and B 
in the figure. 3. The transition from the surface 
layer to the thermocline layer becomes more 
abrupt, i.e., the curvature at the “knee” of 
the trace is increased. One or more of these 
events have been described by various authors 
(SVERDRUP, et al., 1942, p. 497; Anonymous, 
1946; LAFOND, 1954), but to our knowledge 
no data illustrating their occurence have been 
published. 

Francıs and STOMMEL (1953) used fifteen 
BT series obtained at weather ship stations 
before and after (and in some cases during) 
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gales to show that the depth of the thermocline 
increases during a strong wind. These data, 
lent to us by the Woods Hole Oceano- 
graphic Institution, show clearly the events 
described in the preceding paragraph. Selected 
BTs from one of the series appear as Fig. 6. 

What occurs, as a result of the increasing 
wind, might be called “confined mixing” 
because the turbulent mixing induced at the sea 
surface is bounded below by a stable layer. The 
stability of the confining, stable layer is in- 
creased during the process. A concise physical 
explanation of how or why this phenomenon 
occurs, has not been given. 


The maintenance of frontal layers : 


The correlation that appears in the cross 
sections among the isothermal warm water 
mass, the subjacent stable frontal layer and the 
abrupt transition between the two suggests that 
the frontal layers are maintained by the same 
process that maintains sharp thermoclines, i.e., 
confined mixing. The question arises as to the 
source of the turbulence. 

The frontal layer is overlain by a mixed 
layer of neutral stability. Turbulence within 
the mixed layer is not damped by action of the 
Archimedean force. Turbulence created at the 
sea surface by the wind or by cooling and 


Ioo 


subsequent overturn can extend downward to 
the frontal layer. There is another source of 
turbulence which may be of basic importance 
in the maintenance of frontal layers. 

Circulations like those drawn in the warm 
water in each section have been suggested by 
Upa (1938) and have been reported to occur 
in Lake Mead (VETTER, 1953). It is necessary, 
now, to postulate that the “mixed” water 
overlying the frontal layer is not perfectly 
homogeneous. One would expect hetero- 
geneities of density to exist near the surface of 
this water mass by virtue of small scale varia- 
tions in the factors that control the heat and 
water transfer processes at the boundary; for 
example, small scale variations in wind speed. 
If the warm water mass is not perfectly homo- 
geneous, the hypothesized circulation in the 
warm water intermittently will bring parcels 
of low-density water into vertical juxtaposi- 
tion beneath relatively dense water. Instability, 
so produced, will provide a continuing source 
of energy for turbulent overturn. Stommel has 
suggested that closed circulations, required by 
the presence of lateral boundaries, may be of 
importance in producing isothermal layers 
bounded by thermoclines (POLLAK, 1954). A 
front seems to provide a “fluid boundary” 
which allows this process to occur in the open 
ocean. 


Changes in properties across fronts 


Convergent flow appears to be partially 
responsible for the change of properties that 
occurs across a surface front. Surface waters, 
of different properties by virtue of their 
different origins, would be brought into juxta- 
position by flows like those shown in the cross 
sections. Water motion hypothesized to occur 
near the front shown in Fig. 4 leads to an 
additional partial explanation for the change in 
properties that occurs across this front. Both 
the warm and cool water masses are nearly 
isothermal vertically. It is proposed that cell-like 
circulations in the warm and cool water masses 
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lead to turbulent överturn and mixing, as 
already described. This overturn maintains 
vertically mixed layers of differing thickness in 
the warm and cool water. 

Specifically, in Fig. 4, vertical mixing extends 
from the sea surface to a depth greater than 
500 feet in the cool water south of the front and 
to a depth only a little greater than 100 feet in 
the warm water just north of the front. This 
vertical mixing has virtually eliminated the 
vertical temperature gradient. Mixing in the 
warm water north of the front is confined by 
the stable frontal layer. South of the front, 
vertical turbulence has extended downward 
to the main thermocline and it seems reasonable 
that the surface temperature has been lowered 
by mixing with this deeper water. This process 
may contribute to the change in temperature 
observed horizontally across the front. 

The mixing that redistributes the tempera- 
ture gradient will also redistribute gradients 
of other properties and will, in general, contrib- 
ute to the changes in all properties horizontally 
across the front. For example, the resulting 
surface velocities in the cool and warm water 
masses separated by the front will, in general, 
be different, and thus the vertical mixing may 
contribute to the horizontal shear across such 
a front. Analogous reasoning suggests that there 
should be large shear across the frontal layer. 
Shear, then, is envisaged as possibly a result 
of the processes that maintain such a front. 
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The Role of the Gulf Stream in the Prediction of Iceberg 
Distribution in the North Atlantic” 


By LOUIS A. POST, USN Hydrographic Office, Washington, D.C. 


(Manuscript received February I, 1955) 


Abstract 


A twenty-six year correlation is made between the yearly frequency fluctuations of icebergs 
south of the 48th parallel and the preseasonal sea surface temperature anomalies at Key West 
in the Straits of Florida and those of the Labrador Current three years later. 

Good agreement among these variables is attributed to (r) the immediate effect of the relative 
strength of the Gulf Stream in barring the southward flow of bergs, and (2) the effect that 
varying strength of the Gulf Stream has upon the volume of warm water that escapes to the 
north, the effect in turn upon the rate of transport of icebergs, and the potential number 
eventually to be carried into the shipping lanes. 

An iceberg prediction table is offered, making possible a prediction for the current year by 
using the sea temperature anomaly at Key West for March of this year together with the 
anomaly of three years earlier. 

“The opinions and assertions contained herein are the private ones of the author and do not 


reflect the views of the Navy Department or the Naval Establishment.’” 


The drift of icebergs from their birthplace, 
the glaciers of Greenland, southward and 
into the heavily-traveled shipping lanes of the 
North Atlantic Ocean is a menace to the 
safety of navigation, long recognized and 
feared. Only since the “Titanic” sinking in 
1912 after collision with an iceberg has any 
system been in operation for guarding against 
this danger. The International Ice Patrol was 
established to warn passing vessels of the limits 
of danger from day to day throughout each 
iceberg season. Prescribed steamer tracks were 
laid out and their use recommended during 
different parts of each season depending upon 
the numbers of icebergs and their southward 
extent. To facilitate a comparison of the sea- 
sonal and yearly abundance of bergs in the 
most critical area, the number of bergs drifting 
south of latitude 48° N has been compiled to 
obtain monthly and yearly averages. 


The numbers of bergs found south of this 
latitude have varied markedly, from as few 
as O to So per year up to a maximum of 
1,351 in 1929. Although the average is a little 
over 400 for the season, the majority of years 
have totaled either far above of far below 
this value. Investigators have concluded that 
such a phenomenal variation could result from 
fluctuations in several controlling factors: 
wind, currents, pack ice, winter temperature, 
and the production rate of the icebergs them- 
selves, or their rate of calving from the glaciers. 


! Predictions made for the last three iceberg seasons 
have been verified: 


See — Newsweek Magazine, July 13, 1953, 
Science Newsletter, May 1, 1954, 
Science Newsletter, April 30, 1955. 


However, like all new theory worth serious considera- 


tion, many years may be required before complete 
verification. 
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It was apparent that if the proper balance 
between the causative variates and the iceberg 
distribution could be measured, a means for 
forecasting the severity of each season well in 
advance might be obtained. 

A great amount of research and many papers 
have resulted from efforts to determine this 
balance. The work of Smith and of Soule in 
this direction is well known. Smrrx (1928) 
developed a formula for iceberg prediction 
based on the atmospheric pressure distribution 
over Greenland. He has attributed discrepancies 
between calculated numbers of bergs and the 
extreme recorded values to abnormalities, such 
as “ice jams, variations in precipitation, winter 
storms, summer calms etc.” A similar approach 
is taken by SCHELL (1952); yet, since it is 
the extreme rather than the average in iceberg 
distribution which occurs most frequently, 
reliable forecasting techniques must prove 
their ability to correlate well with the extreme 
as well as with average conditions. Because 
no one cause has been shown to account for 
the differences in berg distribution for past 
and predicted years, the opinion is prevalent 
that no one cause is predominant. 

The writer differs with this opinion and 
suggests that the reason no better correlations 
have been found to exist is perhaps a tendency 
to magnify the complexity of related but 
secondary causes instead of locating a single 
cause which alone may be capable of producing 
at least the greater part of the fluctuation in 
iceberg distribution. It is the purpose of this 
paper to provide evidence that fluctuation in 
the Gulf Stream System is such a single cause 
and to derive a prediction technique for 
icebergs south of latitude 48°N from causative 
variates in this system. 


Introduction 


Present-day methods for the prediction of 
the yearly distribution of icebergs in the North 
Atlantic Ocean from atmospheric variables 
have proven neither conclusive in theory nor 
infallible in practice. The present study was 
begun in the belief that in the deeper, less 
transitory influence of ocean currents may be 
found a more reliable means for forecasting 
iceberg distribution, that is, one that is able 
singly to effect the extreme periodic fluctuation 
in this phenomenon. 
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Fig. 1. Shift of the Northern Edge of the Gulf Stream 
during the iceberg season. Typical monthly limits as 
derived from Surface Dynamic Current Charts of the 
International Ice Patrol, from Ice Patrol Bulletins 
(1932—1950). Northern Edge of Gulf Stream is the 
typical position of the minimum dynamic anomaly line 
within the high-speed core in this region. Depths are 
given in fathoms. 


It is further a purpose of this study to de- 
monstrate, and apply to a prediction technique, 
the relation existing between the ocean’s pulse 
in the Florida Straits where it may be felt most 
clearly, and its influence, delayed but preserved 
intact during months and even years until the 
same pulse will swell in certain measure its 
spreading arteries and veins. 

Because of information made available by 
the International Ice Patrol in the vicinity of 
the Grand Banks, and because this region is 
one in which the conjoining Gulf Stream and 
Labrador Current display distinguishable char- 
acteristics of relative heat and transport, it 
is a natural location in which to test empiri- 
cally the above hypothesis and to develop 
therefrom a method for forecasting the sever- 
ity of the iceberg season. 

Figure 1 illustrates the seasonal shift of 
the northern edge of the Gulf Stream with 
change in its “relative” velocity off the 
Grand Banks. The expression “relative” is 
used because it is not certain to what extent 
the Gulf Stream and the Labrador Current 
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are effective in controlling this shift. However, 
the northward displacement between March 
and June closely conforms to increasing 
velocity of the Gulf Stream from monthly 
drift observations (U. S. Navy HyDROGRAPHIC 
OFFICE, 1944). Arrows indicate the typical 
paths icebergs follow during different parts 
of the season. It is apparent that increasing 
“relative” strength of the Gulf Stream acts 
to reduce the southward progression of bergs, 
and Ice Patrol records frequently show that 
when the westward salient in June towards 
the Banks is well developed or occurs earlier 
in the season, the supply of bergs is effectively 
cut off. 

Now it would appear that if an entire season 
or critical month should show this “relative” 
strength of the Gulf Stream to be above 
average, a light iceberg season might result; 
but if it should be below the average, more 
bergs might be permitted to advance to a 
more southerly latitude. 

The writer selected Key West in the Straits 
of Florida as a point most suitable from its 
location and availability of data for the 
establishment of an empirical correlation 
between the yearly iceberg concentration 
south of latitude 48° N and contemporary 
changes in the strength of the Gulf Stream. 

Here, the major portion of the current is 
closely confined by banks of earth so that its 
fluctuation in velocity should provide a truly 
reliable measure of fluctuation in mass trans- 
port. Further, the Gulf Stream in this region is 
far enough “upstream” from the Grand 
Banks to allow sufficient leeway in preparing 
forecast information, yet is an integral part, 
not only of the North Atlantic Eddy, but 
also of the entire “feedback” circulation 
systems of the Arctic and of Baffin Bay. It is 
likely, further, that in this region the Gulf 
Stream will reflect in definite proportion to 
their ultimate effects fluctuations of far- 
reaching influence occurring from causes 
elsewhere in the system. 

Figure 2 displays the general area of the 
study and the cycle of events which forms the 
basis for the correlation. The results of this 
study indicate: (1) An increase in transport 
of the Gulf Stream measured at (A) produces 
changes of like sign at (B) in the critical region 
of iceberg advance. (2) The immediate effect 
at (B) is a strengthened current and increased 
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Fig. 2. Area of Study. Arrows with heavy solid line 
indicate warming current. Arrows with dashed line 
indicate cooling current. Transport times are estimates 
based on SCHOTT, 1942, United States Coastguard and 
U.S. Navy Hydrographic. Office, 1944. 

resistance to the invasion of icebergs into the 
shipping lanes to the south, favoring a lighter 
iceberg count (SOULE, 1949). (3) An increase 
in volume of the current at (B) after some 
time increases the volume of warm Irminger 
Current water meeting the cold East Greenland 
Current off Kap Farvel at (C). Later an 
increased volume of warmer than average 
water reaches (D) in the region of berg forma- 
tion, favoring reduction in ice cover, greater 
solar heating of the surface water and of the 
overlying atmosphere, and calving from the 
glaciers of a larger than usual number of 
bergs. (4) The influx of a greater volume of 
warm water into Baffin Bay must result in a 
correspondingly greater efflux on the Labrador 
side, accompanied by a withdrawal from the 
Bay of greater amounts of icebergs and sea 
ice stripped away by the current. 

While evidence has accumulated (SmirH, 
1940) that currents influence both the pro- 
duction and transport of pack ice and icebergs, 
in the present study the emphasis is placed 


Tellus VIII (1956), 1 


THE GULF STREAM IN THE PREDICTION OF ICEBERGS 


nos ICEBERGS 


NO OF ICEBERGS (MARCH - JUNE) 


105 


SEA LEVEL DIFFERENCE 
IN FEET 
BERMUDA — CHARLESTON 


DEPARTURE FROM MEAN 


À TEMP ANOMALIES 


' 
w 


+5 


(38, à ER 
ewe 


1 
+ 


L 
2 BUS 3% 37 3 39 
YEAR 


1926 27 2 2 M = 


40 « 


42 43 44 43 46 47 48 49 50 Si 52 


Fig. 3. Number of icebergs south of 48th parallel versus sea temperature anomalies off 
Key West (1926—1952). Solid line indicates number of icebergs. Broken line indicates 
mean sea temperature anomalies (March). Anomalies (1941—1946) are from Daytona Beach, 
Fla. because of absence of data at Key West during these years. Light dotted line indi- 
cates the departure from mean sea level difference in feet, Bermuda minus Charleston, 
SOULE (1940). Sea surface temperature data are from U. S. Coast & Geodetic Survey Pub. TW-I. 


on the effect of changes in the transport of 
currents in the same area in which the berg 
count is made upon the local distribution of 
icebergs. 

The study of the effect of the Gulf Stream 
upon iceberg distribution is divided below 
into two parts: the immediate effect, or the 
effect of the Gulf Stream in inhibiting the 
southward drift of bergs in the Grand Banks 
area; and the delayed effect of changes in 
the strength of the Gulf Stream through its 
extension into Baffin Bay upon the eventual 
strength of the Labrador Current and upon 
the numbers of icebergs it is able to carry 
southward. 


The Immediate Effect 


Figure 3 presents a twenty-six year correla- 
tion between the number of icebergs counted 
south of parallel 48° N during the iceberg 
season (March—June) and preseasonal sea 
surface temperature anomalies at Key West 
in the Straits of Florida. The temperature data 
are those of the U.S. Coast and Geodetic 
Survey’s monthly averages of daily observa- 
tions at tidal gauge stations. An anomaly of 
mean monthly temperature, based on the con- 
cepts given below, should therefore reflect 
proportional change in the strength of the 
current and provide an extensive record of 
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this change in the absence of any long-period 
direct measurements of currents. Since three 
months is the maximum time lag for current 
between the Straits of Florida and the Grand 
Banks, fluctuations occurring in the Straits 
earlier than January (three months prior to 
the opening of the iceberg season) may be 
considered ineffective. Actually the March 
anomalies agreed the best and these have been 
used after careful comparison with other 
months has made possible two conclusions: 
(1) that the changes reflected in the March 
anomalies persisted at Key West for three 
to six months, and (2) that the March anom- 
alies fitted the time lag in such a way as to 
arrive in the Grand Banks area at the most 
critical part of the iceberg season. The average 
percentage of bergs by months, based on so 
years of Ice Patrol data (1900—1950), is: 


.| a 

Oo} ) cs ea aes 
s[2| 5181882) % 8/5) 5] 8] Month 
Sele] a] al Sloat an} ol] ala 
1/4] 2| rz) 24) 35|27| 5| 2| rer. rl Percent 


The time lag required for fluctuations in 
the Florida Straits to reach the vicinity of 
Grand Banks and be effective the same season 
would thus appear to be between 1 and 3 
months. 
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A plausible interpretation of the correlation 
in figure 3 is that higher than normal tem- 
peratures of the water at Key West will be 
followed by a weaker than normal Gulf 
Stream offthe Grand Banks in order to permit 
a greater southward extent of icebergs and 
vice versa. This does not necessarily mean 
that high temperatures at Key West accompany 
a weaker current at that point. Two concepts 
may be advanced to explain the mechanisms 
involved in effecting the correlation found in 
figure 3. 

The first is that advanced by Iserın (1938) 
and (1940) according to which the edge of 
the current system contracts with increasing 
speed and expands with decreasing speed. 

Under the circumstances prevailing in 
the Straits of Florida “only a shallow band of 
coastal water separates the current from the 
shore’—and “the main thermocline at the 
western margin of the Florida Current is free 
to occupy whatever depth is required by the 
strength of the current”. Thus, a stronger 
current would tend to permit a lesser amount 
of warm water to invade the coastal water, 
whereas a weaker current would reverse this 
tendency (IsELIN, 1940, figures 26-B and B’). 

The other concept assumes that the strength 
is transmitted directly throughout the current 
system. This would mean that high tempera- 
tures at Key West accompany weaker current, 
and that this weakness, transmitted to the 
Grand Banks, is reflected in a larger berg 
count. 

Such indirect current measurements as are 
available tend to support both these concepts. 
The dotted line in figure 3, representing the 
sea level difference, Bermuda-Charleston (1926 
—1938) shows at least a qualitative inverse 
relationship between current and temperature. 

Thus it may be inferred that positive tem- 
perature anomalies in the Straits of Florida 
indicate a weaker current there and vice 
versa. The immediate effect of Gulf Stream 
water of higher speeds and lower tempera- 
tures (at Key West) in barring the south- 
ward flow of bergs via the Labrador Current 
is apparent. 


The Delayed Effect 


The immediate effect described above, 
pronounced over the period of record, appears 
to establish the responsiveness of the strength 
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Fig. 4. Time Lag — Labrador Current, West Green- 

land Current (minus 2 years) and Gulf Stream (minus 

3 yars) Solid line: Current Transport. Broken line: 

Current temperature. Light dashed line: Icebergs. La- 

brador Current and West Greenland Current are drawn 
from SOULE and Barnes (1941). 


of the Gulf Stream off the Grand Banks to 
changes one to three months earlier in the 
Florida Straits. However, as stated previously, 
this strength is a relative one. 

The question remains: can the immediate 
response, the primary influence of the Gulf 
Stream upon iceberg distribution, be further 
traced throughout the branches of the system 
until ultimately it makes itself felt again as a 
primary source of variation in the influence 
of the Labrador Current returning at the 
Grand Banks to its parent stream? 

The estimated time lag of three years 
between a change in the current in the Straits 
of Florida and consequent change in the 
Labrador Current via the Irminger Current 
and Baftin Bay is only approximate, but it 
seems likely that changes of several months 
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Fig. s. Correlation—Icebergs south of 48° N versus Gulf Stream Anomalies (°F.) 
of present year and present year minus 3. 


duration in the Florida Straits might be detect- 
able three years later if indeed they continue 
to be the predominant cause of change in the 
Labrador Current. 

It has been inferred that high temperatures 
at Key West precede a discharge of weak 
current to the north and conversely. The 
effect of a weakened Gulf Stream has been 
shown in figure 3 to favor a greater number of 
icebergs; correspondingly, the effect of a 
weakened Labrador Current would favor a 
decrease in this number. Therefore, tempera- 
tures at Key West must be in inverse relation 
to the iceberg population approximately three 
years later if they really reflect Gulf Stream 
fluctuations which influence the Labrador 
Current. 

Figure 4 supports this hypothesis by dem- 
onstrating the agreement between (1) Lab- 
rador Current volume and temperature anom- 
alies, measured at South Wolf Island, Lab- 
rador, (2) West Greenland Current volume 
and temperature anomalies at Kap Farvel two 
years earlier, and (3) Gulf Stream current and 
Key West temperature anomalies three years 
earlier, together with the resulting changes in 
iceberg distribution. 
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A conclusion previously drawn by SouLe 
(1940) is that the above “negative correlation 
between the departure from normal of the 
mean temperature and volume of flow of the 
West Greenland Current” implies that these 
fluctuations “are largely the result of fluctua- 
tions in the East Greenland Current rather 
than the Irminger Current’. The fact that 
figure 4 shows a positive correlation between 
Labrador Current speed and temperature 
appears to contradict this conclusion, because 
only by an increase in the net heat transfer 
via the Irminger Current could the Labrador 
Current receive an increase in both volume 
and temperature. Moreover, it appears ex- 
tremely unlikely that Gulf Stream fluctuations 
could have traversed the entire polar basin 
and reached the vicinity of Kap Farvel in 
the single year which separates these pheno- 
mena. 

‘A more suitable interpretation might be 
that although an initially stronger and colder 
Gulf Stream can result in a stronger and colder 
Irminger Current off Kap Farvel and even in a 
greater volume of colder than normal West 
Greenland Current, the total heat transfer 
into Baffin Bay will be greater thereby, and 
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Fig. 6. Relation between number of icebergs south of 48° N and 
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the result will be an increase in both the speed 
and temperature of the Labrador Current. 

The agreement between the observed Lab- 
rador Current and that inferred from figure 
4-C thus strongly suggests that anomalies of 
the former reflect, and were induced by, 
anomalies of the Gulf Stream System and that 
these are measurable, even three years in 
advance, in the Straits of Florida. A correlation 
coefficient of +. 72 was obtained between the 
transport of the Labrador Current and Gulf 
Stream three years earlier as given in figure 4. 

Figure 5 then demonstrates the correlation 
obtained between iceberg count and Gulf 
Stream temperature anomalies of the present 
year and three years earlier. The effect of 
positive anomalies (weak current) in the pres- 
ent year in influencing larger numbers of 
bergs in conjunction with negative anomalies 
(strong current) of three years earlier is 
apparent in figure s-A. The summation of 
the effect of anomalies of the present year 
and present year minus three upon iceberg 
distribution is given in figure 5-B. Correlation 
coefficients listed below show the linear relation 
between iceberg count and temperature anom- 
alies of the Gulf Stream for the present year 
alone and in combination with preceding years 
as indicated. 


Correlation Coefficients between iceberg count and 
Key West temperature anomalies for: 


Summation of Present Year and Present 


Present Year minus: 
Year 
2: |,2. ]»351,»4: Lomeli elie 
+.65 |+-.42|+-55|+-75|+-65|+-46|+-.35|+-49| +-54 


Successive Years without Present Year. 


Present Present Year minus: 
Year’ RE eee 
BEIRBZEUFUTLENNGZ: 
+.65 er -.76|--.53|+-43| +-51|--27|- 37 


It is evident that a more reliable correlation 
than that obtained for the present year alone 
results only by combining for each year of 
record the effect upon iceberg distribution 
of the Gulf Stream anomalies both of the 
present year and of three years earlier. These 
values for the years of record are correlated 
in figure 6, and probable predictable limits in 
iceberg distribution are indicated. These limits, 
as drawn, depart slightly from the linear 
relationship assumed in figure 5 in order to 
provide for the effect of theoretical extremes 
in the strength of the Gulf Stream and Labrador 
Current. Figure 7 correlates the numbers of 
icebergs observed (1926—1951) with numbers 
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Fig. 7. Comparison between observed (solid) and computed (dashed line) numbers 
of icebergs appearing south of 48° N, computed from figure 6. 


computed from figure 6 by using March 
temperature anomalies at Key West for the 
present year and present year minus three. 
A coefficient of + .94 was obtained between 
computed and observed numbers grouped as 
in figure 7 which allows a leeway of + 100 
bergs. 

The agreement supports the hypothesis that 
the Gulf Stream influence is the predominant 
one. The averages obtained in figure 6 are 
incorporated into an iceberg prediction table, 
table I. Monthly temperature means are ob- 
tained from “Surface Water Temperatures, 
Atlantic Coast”, U.S.C. & G.S. Pub. No. 
TW-1. Anomalies during the current year 
may be obtained directly from the Coast and 
Geodetic Survey at the end of each month. 

Iceberg predictions using table I should be 
interpreted with caution, despite the satis- 
factory correlation coefficients. For example, 
it will be noted that the observed and computed 
values in figure 7 are in disagreement for the 
year 1949. The dearth of icebergs that year is 
attributed (SOULE, 1949) to very abnormal 
winds that delayed the progress of the bergs 
until, with the advance of the warm season, 
they failed to survive the journey to the 48th 
parallel. 
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However, the present study has shown that 
such occurrences have been infrequent over 
an extended period and that the primary role 
in controlling the distribution of icebergs 
south of 48° N is indeed that of the Gulf 
Stream together with its branches and return- 
flow tributaries. 

Thus, paradoxically, the Gulf Stream plays 
the role of abettor to the threat of icebergs 
by its very effort to restrain them. For the 
greater its volume southward of the Grand 
Banks, though it impedes the progress of bergs 
temporarily, the greater must be the volume 
of water to reach Baffin Bay, and the greater 
ultimately must be the volume of the Labrador 
Current three years later. 


Summary 


1. The drift of icebergs south of latitude 
48° N is shown to be controlled seasonally by 
thé relative strengths of the Gulf Stream and 
the Labrador Current in that vicinity. 

2. The Straits of Florida is selected as 
ideally located for the measurement of anom- 
alies in the strength of the Gulf Stream. Key 
West temperature anomalies are shown to 
be related inversely to current speed and to 
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Table 1. Iceberg prediction table 


Number of Icebergs South of 48° N from Temp. (Present Year) and Temp. (Present Year—3) 
Mean for March 


Present Year—Gulf Stream Temperature (°F) Key West 


Mean 

| SRE EEE EE tn Ep Fe N 
5 : o | fe} | (6) | fo) | fo) | 100 100 100 350 350 350 600 
a 5 fo) | o | o | o | o | 100 | 100 | 350 | 350 | 350 | 500 | 600 
2 : o | o | o | o | 100 | 100 | 350 | 350 | 350 | 500 | 600 | 600 
3 : o | o | fo) | 100 | 100 | 200 | 350 | 350 | 450 | 600 | 600 | 600 
3 : o | fo) | 100 | 100 | 200 | 350 | 350 | 350 | 600 | 600 | 600 | 600 
= o | 100 | 100 | 100 | 350 | 350 | 450 | 600 | 600 | 700 | 800 | 800 
3 € fe) | 100 | 100 | 350 | 350 | 450 | 600 | 600 | 700 | 800 | 800 | 800 
E - 100 | 100 | 200 | 350 | 350 | 600 | 600 | 800 | 800 | 800 | 900 | >900 
jr : 100 | 100 | 350 | 350 | 600 | 600 | 800 | 800 | 800 | >900 | >900 | >900 
8 71.8 | 100 | 200 | 350 | 500 | 600 | 700 | 800 | 800 | >900 | >900 | >900 | >900 
8 70.8 | 100 | 350 | 350 | 600 | 600 | 800 | 800 | 900 | > 900 | >900 | >900 | >900 
£ 69.8 | 200 | 350 | 350 | 600 | 600 | 800 | 800 | >900 | >900 | >900 | >900 | >900 


For iceberg number = 100 correlation cofficient is +.94 (see figure 7) 


be synchronized with current fluctuation at 4. Current fluctuation in the Labrador 
the Grand Banks less than three months Current is found to be derived from the Gulf 
later. Stream and is shown to be proportional to 


3. A curve is drawn relating temperature fluctuations in the West Greenland Current 
anomalies at Key West to the iceberg count for two years earlier and to those of the Gulf 
a period of twenty-six years. It is found Stream at Key West three years earlier. 
that high temperature and weak current 5. Summation of the temperature anomalies 
coexist and effect a proportional weakening for the present year and three years earlier at 
in the “relative” Gulf Stream strength off Key West results in a closer correlation with 
Grand Banks together with increase in the iceberg distribution. The numbers of bergs 
southerly drift of icebergs one to three months resulting from both anomalies are incorporated 
later. in a prediction table. 
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Current Data on the Chemical Composition of Air and Precipitation V 


(For further information see Egnér, H., Eriksson., Tellus 7, 134—139, 195 5.) 
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CO,-values Sept.—Dec. 1955 in Scandinavia. 
(Cf. Fonselius, Koroleff: Tellus, 7, 2, pp. 258—265). . 
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ERTTERS-TO THE SEDITOR 
On »Forecasting with the Barotropic Model» 


Dear Sir, 

The recent article by Bouin (1955) contains among 
other results also a series of four 48-hour forecasts 
made by O. Haug and obtained, as it is said, by the 
use of Fjortoft’s “graphical” method. The actual 
procedure is, however, not explained. Reference is 
merely made to two articles by the author (FJORTOFT, 
1952, 1955). From information given me after the 
appearance of Bolin’s article it is clear that the 
method used by Haug is different from the method 
developed in the above mentioned articles. It is 
similar to the method proposed in the most recent 
article but differs in one important respect. 

There might be different numerical approaches to 
one and the same problem which in the limit ap- 
proach the true solution. Any such method, how- 
ever, suffers from certain deficiencies which neces- 
sarily lead to errors and the various causes to these 
errors are not very easy to separate. It is therefore 
useful to compare integration results obtained by 
different methods but it is also clear that a detailed 
discussion of the results is very necessary. I have pre- 
pared 48-hour forecasts using the method developed 
in my most recent article for some of the maps from 
which numerical forecasts with BESK have been 
done and were reported by Bolin. In view of table 2 
in Bolin’s article it seems to me important to ask 


Dear Sir, 

In TELLUS, Vol 7, BercrHorsson, Döös, et al 
(1955) report a series of numerical forecasts made on 
the Swedish computer BESK and compared with 
subjective forecasts made for the same dates by con- 
ventional methods. By way of verification, the 
authors list correlation coefficients between fore- 
cast and observed soomb height changes over a 
selected area and also subjective grades assigned to 
the prognoses by a board of five meteorologists. 

It is of interest to note that the correlation between 
the two sets of verification scores is only +0.35 e- 
valuated by O. Haug, Det Norske Meteorologiske In- 
stitutt). This may not be so surprising, since the first 
measure (the correlation) describes the success in 
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permission from the editor to get the results from 
these forecasts published, asking simultaneously the 
reader to excuse that the discussion is deferred to a 
later occasion. In the first line of the table below the 
correlation coefficient r between computed and ob- 
served 48-hour height changes is given together 
with the mean absolute error as obtained graphically. 
The second line shows the corresponding values 
obtained numerically with BESK. The verification 
was done for an area denoted as “large” in Bolin’s 
article and the values for the numerical forecasts are 
therefore somewhat different from those reported by 
Bolin. The table shows that there is no significant 
difference, on an average, in the results obtained 
graphically and numerically in these three cases. 
Further comparisons of this kind will be made. 


Sept. 
26—28 


Sept. 
27729 


Sept. 
28— 30 


1 | € V | € LE | € 


Fjortoft | 0.82 | 46m | 0.71 | 51m | 0.78 | 52m 


BESK | 0.83 | 44m | 0.53 | 60 m | 0.80 | 59m 
Ragnar Fjortoft 
Meteorological Institute of Norway 


predicting a soomb pattern after having started from 
a given pattern, whereas the subjective grade more 
nearly defines the closeness of agreement between 
the prognosis and the observed chart. 

The purpose of this note, however, is to point 
out the implications of a 2-way distribution of these 
verification scores, which is given below in Table x. 
The table classiffies the scores of the 24h, 48h, and 
72h numerical prognoses and the 24h and 48h scores 
of the conventional, subjective progs. The subjec- 
tive scores are from 5 to 1, with 5 standing for best 
agreement between predicted and observed. The spe- 
cification of class intervals conforms to the simple 
frequency rule: 


116 
SUBJECTIVE SCORES 
TN Re: 
Tr > +o.8o 18 15 [_4 0 o 
CORRELATION °7° — 0.80 [ro 19 17/5 ı 
COEFFICIENT 9.46 — 0.69 6|ı3_ 16 4 


8 
0.25 — 0.45 TOURS DID RL 
2% 22025 he at a 
Table 1. Contingency table, correlation coefficients 

vs subjective scores. 168 cases. 


that is, the intervals are chosen so as to make ap- 
proximately equal the frequency totals of Row 1 and 
Column 1, Row 2 and Column 2, and so on. 

It is apparent from inspecting Table 1 that the 
sample is skewed. The total of frequencies on either 
side of the upper-left-to-lower-right diagonal are 
approximately equal, which is an outcome of the 
way in which class intervals have been specified. 
Looking, however, at the extremes of the distri- 
bution (the lower left and upper right corners of the 
table as shown by the diagonal lines), it is seen that 
there are more prognoses for which the correlation 
was very low and the subjective score very high, 
than vice versa. In 4 cases, for example, the progs 
were given a top score yet they showed a correlation 
of height change, of less then +0.25. At the other 
extreme, on the other hand, there were no cases 
when the correlation was above + 0.80 and the fore- 
casts were graded below a 3. 
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Clearly, the low correlation between the sets of 
verifications (+0.35) is attributable to the dispersion 
of frequencies in Table 1. It is also suggested, how- 
ever, that this low correlation is due more to those 
cases in which the height change correlation was 
very low but the score was very high, than to those 
cases in which the height change correlation was 
very high and the score low. 

The implication here is of some interest, namely, 
that although the correlation between forecast and 
observed height changes.may be criticized as a meas- 
ure of prognostic accuracy, or prognostic value, it 
does tend to act as a lower limit as far as concerns 
measuring the similarity between the prog and the 
observed chart. Thus the similarity tends to be at 
least as good as is indicated by the correlation coefti- 
cient and in fact may be much better. 


William H. Best 
University of Stockholm 
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